Protocol to dissociate epithelia from non-pregnant and pregnant mouse cervical tissue for single-cell RNA-sequencing

A challenge in studying cervical epithelial cell biology at the single-cell level is that differentiated subtypes, in particular mucus-secreting goblet cells, are sensitive to disassociating enzymes making isolation of all epithelial subpopulations difficult. Here we present a protocol to dissociate epithelia from non-pregnant and pregnant mouse cervical tissue for single-cell RNA-sequencing. We describe steps for harvesting cervices, preparing cervical tissue, dissociation of cervical cells, and viability checks. We then detail library preparation, sequencing, and procedure for data analysis. For complete details on the use and execution of this protocol, please refer to Cooley et al. (2023).1.

Enhancer regulatory networks globally connect non-coding breast cancer loci to cancer genes

Genetic studies have associated thousands of enhancers with breast cancer. However, the vast majority have not been functionally characterized. Thus, it remains unclear how variant-associated enhancers contribute to cancer. Here, we perform single-cell CRISPRi screens of 3,512 regulatory elements associated with breast cancer to measure the impact of these regions on transcriptional phenotypes. Analysis of >500,000 single-cell transcriptomes in two breast cancer cell lines shows that perturbation of variant-associated enhancers disrupts breast cancer gene programs. We observe variant-associated enhancers that directly or indirectly regulate the expression of cancer genes. We also find one-to-multiple and multiple-to-one network motifs where enhancers indirectly regulate cancer genes. Notably, multiple variant-associated enhancers indirectly regulate TP53. Comparative studies illustrate sub-type specific functions between enhancers in ER+ and ER- cells. Finally, we developed the pySpade package to facilitate analysis of single-cell enhancer screens. Overall, we demonstrate that enhancers form regulatory networks that link cancer genes in the genome, providing a more comprehensive understanding of the contribution of enhancers to breast cancer development.

Large-scale production of human blastoids amenable to modeling blastocyst development and maternal-fetal cross talk

Recent advances in human blastoids have opened new avenues for modeling early human development and implantation. One limitation of our first protocol for human blastoid generation was relatively low efficiency. We now report an optimized protocol for the efficient generation of large quantities of high-fidelity human blastoids from naive pluripotent stem cells. This enabled proteomics analysis that identified phosphosite-specific signatures potentially involved in the derivation and/or maintenance of the signaling states in human blastoids. Additionally, we uncovered endometrial stromal effects in promoting trophoblast cell survival, proliferation, and syncytialization during co-culture with blastoids and blastocysts. Side-by-side single-cell RNA sequencing revealed similarities and differences in transcriptome profiles between pre-implantation blastoids and blastocysts, as well as post-implantation cultures, and uncovered a population resembling early migratory trophoblasts during co-culture with endometrial stromal cells. Our optimized protocol will facilitate broader use of human blastoids as an accessible, perturbable, scalable, and tractable model for human blastocysts.

Modeling post-implantation stages of human development into early organogenesis with stem-cell-derived peri-gastruloids

In vitro stem cell models that replicate human gastrulation have been generated, but they lack the essential extraembryonic cells needed for embryonic development, morphogenesis, and patterning. Here, we describe a robust and efficient method that prompts human extended pluripotent stem cells to self-organize into embryo-like structures, termed peri-gastruloids, which encompass both embryonic (epiblast) and extraembryonic (hypoblast) tissues. Although peri-gastruloids are not viable due to the exclusion of trophoblasts, they recapitulate critical stages of human peri-gastrulation development, such as forming amniotic and yolk sac cavities, developing bilaminar and trilaminar embryonic discs, specifying primordial germ cells, initiating gastrulation, and undergoing early neurulation and organogenesis. Single-cell RNA-sequencing unveiled transcriptomic similarities between advanced human peri-gastruloids and primary peri-gastrulation cell types found in humans and non-human primates. This peri-gastruloid platform allows for further exploration beyond gastrulation and may potentially aid in the development of human fetal tissues for use in regenerative medicine.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

CHD-associated enhancers shape human cardiomyocyte lineage commitment

Enhancers orchestrate gene expression programs that drive multicellular development and lineage commitment. Thus, genetic variants at enhancers are thought to contribute to developmental diseases by altering cell fate commitment. However, while many variant-containing enhancers have been identified, studies to endogenously test the impact of these enhancers on lineage commitment have been lacking. We perform a single-cell CRISPRi screen to assess the endogenous roles of 25 enhancers and putative cardiac target genes implicated in genetic studies of congenital heart defects (CHDs). We identify 16 enhancers whose repression leads to deficient differentiation of human cardiomyocytes (CMs). A focused CRISPRi validation screen shows that repression of TBX5 enhancers delays the transcriptional switch from mid- to late-stage CM states. Endogenous genetic deletions of two TBX5 enhancers phenocopy epigenetic perturbations. Together, these results identify critical enhancers of cardiac development and suggest that misregulation of these enhancers could contribute to cardiac defects in human patients.

Dynamic states of cervical epithelia during pregnancy and epithelial barrier disruption

The cervical epithelium undergoes changes in proliferation, differentiation, and function that are critical to ensure fertility and maintain pregnancy. Here, we identify cervical epithelial subtypes in non-pregnant, pregnant, and in labor mice using single-cell transcriptome and spatial analysis. We identify heterogeneous subpopulations of epithelia displaying spatial and temporal specificity. Notably in pregnancy, two goblet cell subtypes are present in the most luminal layers with one goblet population expanding earlier in pregnancy than the other goblet population. The goblet populations express novel protective factors and distinct mucosal networks. Single-cell analysis in a model of cervical epithelial barrier disruption indicates untimely basal cell proliferation precedes the expansion of goblet cells with diminished mucosal integrity. These data demonstrate how the cervical epithelium undergoes continuous remodeling to maintain dynamic states of homeostasis in pregnancy and labor, and provide a framework to understand perturbations in epithelial health that increase the risk of premature birth.

Global chromatin landscapes identify candidate noncoding modifiers of cardiac rhythm

Comprehensive cis-regulatory landscapes are essential for accurate enhancer prediction and disease variant mapping. Although cis-regulatory element (CRE) resources exist for most tissues and organs, many rare - yet functionally important - cell types remain overlooked. Despite representing only a small fraction of the heart's cellular biomass, the cardiac conduction system (CCS) unfailingly coordinates every life-sustaining heartbeat. To globally profile the mouse CCS cis-regulatory landscape, we genetically tagged CCS component-specific nuclei for comprehensive assay for transposase-accessible chromatin-sequencing (ATAC-Seq) analysis. Thus, we established a global CCS-enriched CRE database, referred to as CCS-ATAC, as a key resource for studying CCS-wide and component-specific regulatory functions. Using transcription factor (TF) motifs to construct CCS component-specific gene regulatory networks (GRNs), we identified and independently confirmed several specific TF sub-networks. Highlighting the functional importance of CCS-ATAC, we also validated numerous CCS-enriched enhancer elements and suggested gene targets based on CCS single-cell RNA-Seq data. Furthermore, we leveraged CCS-ATAC to improve annotation of existing human variants related to cardiac rhythm and nominated a potential enhancer-target pair that was dysregulated by a specific SNP. Collectively, our results established a CCS-regulatory compendium, identified novel CCS enhancer elements, and illuminated potential functional associations between human genomic variants and CCS component-specific CREs.

Computational identification of clonal cells in single-cell CRISPR screens

BACKGROUND

Single-cell CRISPR screens are powerful tools to understand genome function by linking genetic perturbations to transcriptome-wide phenotypes. However, since few cells can be affordably sequenced in these screens, biased sampling of cells could affect data interpretation. One potential source of biased sampling is clonal cell expansion.

RESULTS

Here, we identify clonal cells in single cell screens using multiplexed sgRNAs as barcodes. We find that the cells in each clone share transcriptional similarities and bear segmental copy number changes. These analyses suggest that clones are genetically distinct. Finally, we show that the transcriptional similarities of clonally expanded cells contribute to false positives in single-cell CRISPR screens.

CONCLUSIONS

Experimental conditions that reduce clonal expansion or computational filtering of clonal cells will improve the reliability of single-cell CRISPR screens.

Single-Cell Genomics: Catalyst for Cell Fate Engineering

As we near a complete catalog of mammalian cell types, the capability to engineer specific cell types on demand would transform biomedical research and regenerative medicine. However, the current pace of discovering new cell types far outstrips our ability to engineer them. One attractive strategy for cellular engineering is direct reprogramming, where induction of specific transcription factor (TF) cocktails orchestrates cell state transitions. Here, we review the foundational studies of TF-mediated reprogramming in the context of a general framework for cell fate engineering, which consists of: discovering new reprogramming cocktails, assessing engineered cells, and revealing molecular mechanisms. Traditional bulk reprogramming methods established a strong foundation for TF-mediated reprogramming, but were limited by their small scale and difficulty resolving cellular heterogeneity. Recently, single-cell technologies have overcome these challenges to rapidly accelerate progress in cell fate engineering. In the next decade, we anticipate that these tools will enable unprecedented control of cell state.

Blastocyst-like structures generated from human pluripotent stem cells

Limited access to embryos has hampered the study of human embryogenesis and disorders that occur during early pregnancy. Human pluripotent stem cells provide an alternative means to study human development in a dish^1-7. Recent advances in partial embryo models derived from human pluripotent stem cells have enabled human development to be examined at early post-implantation stages^8-14. However, models of the pre-implantation human blastocyst are lacking. Starting from naive human pluripotent stem cells, here we developed an effective three-dimensional culture strategy with successive lineage differentiation and self-organization to generate blastocyst-like structures in vitro. These structures-which we term 'human blastoids'-resemble human blastocysts in terms of their morphology, size, cell number, and composition and allocation of different cell lineages. Single-cell RNA-sequencing analyses also reveal the transcriptomic similarity of blastoids to blastocysts. Human blastoids are amenable to embryonic and extra-embryonic stem cell derivation and can further develop into peri-implantation embryo-like structures in vitro. Using chemical perturbations, we show that specific isozymes of protein kinase C have a critical function in the formation of the blastoid cavity. Human blastoids provide a readily accessible, scalable, versatile and perturbable alternative to blastocysts for studying early human development, understanding early pregnancy loss and gaining insights into early developmental defects.

Identification and characterization of cellular heterogeneity within the developing renal interstitium

Kidney formation requires the coordinated growth of multiple cell types including the collecting ducts, nephrons, vasculature and interstitium. There is a long-held belief that interactions between progenitors of the collecting ducts and nephrons are primarily responsible for kidney development. However, over the last several years, it has become increasingly clear that multiple aspects of kidney development require signaling from the interstitium. How the interstitium orchestrates these various roles is poorly understood. Here, we show that during development the interstitium is a highly heterogeneous patterned population of cells that occupies distinct positions correlated to the adjacent parenchyma. Our analysis indicates that the heterogeneity is not a mere reflection of different stages in a linear developmental trajectory but instead represents several novel differentiated cell states. Further, we find that β-catenin has a cell autonomous role in the development of a medullary subset of the interstitium and that this non-autonomously affects the development of the adjacent epithelia. These findings suggest the intriguing possibility that the different interstitial subtypes may create microenvironments that play unique roles in development of the adjacent epithelia and endothelia.

Rational Reprogramming of Cellular States by Combinatorial Perturbation

Ectopic expression of transcription factors (TFs) can reprogram cell state. However, because of the large combinatorial space of possible TF cocktails, it remains difficult to identify TFs that reprogram specific cell types. Here, we develop Reprogram-Seq to experimentally screen thousands of TF cocktails for reprogramming performance. Reprogram-Seq leverages organ-specific cell-atlas data with single-cell perturbation and computational analysis to predict, evaluate, and optimize TF combinations that reprogram a cell type of interest. Focusing on the cardiac system, we perform Reprogram-Seq on MEFs using an undirected library of 48 cardiac factors and, separately, a directed library of 10 epicardial-related TFs. We identify a combination of three TFs, which efficiently reprogram MEFs to epicardial-like cells that are transcriptionally, molecularly, morphologically, and functionally similar to primary epicardial cells. Reprogram-Seq holds promise to accelerate the generation of specific cell types for regenerative medicine.

A Cellular Anatomy of the Normal Adult Human Prostate and Prostatic Urethra

A comprehensive cellular anatomy of normal human prostate is essential for solving the cellular origins of benign prostatic hyperplasia and prostate cancer. The tools used to analyze the contribution of individual cell types are not robust. We provide a cellular atlas of the young adult human prostate and prostatic urethra using an iterative process of single-cell RNA sequencing (scRNA-seq) and flow cytometry on ~98,000 cells taken from different anatomical regions. Immunohistochemistry with newly derived cell type-specific markers revealed the distribution of each epithelial and stromal cell type on whole mounts, revising our understanding of zonal anatomy. Based on discovered cell surface markers, flow cytometry antibody panels were designed to improve the purification of each cell type, with each gate confirmed by scRNA-seq. The molecular classification, anatomical distribution, and purification tools for each cell type in the human prostate create a powerful resource for experimental design in human prostate disease.

Using single-cell RNA sequencing, immunofluorescence, and flow cytometry, Henry et al. create a cellular anatomy of the normal human prostate and provide the tools to identify, isolate, and localize every cell type. They identify two additional epithelial cell types enriched in the prostatic urethra and proximal prostatic ducts.

Global Analysis of Enhancer Targets Reveals Convergent Enhancer-Driven Regulatory Modules

Single-cell screens enable high-throughput functional assessment of enhancers in their endogenous genomic context. However, the design of current studies limits their application to identifying the primary gene targets of enhancers. Here, we improve the experimental and computational parameters of single-cell enhancer screens to identify the secondary gene targets of enhancers. Our analysis of >500 putative enhancers in K562 cells reveals an interwoven enhancer-driven gene regulatory network. We find that enhancers from distinct genomic loci converge to modulate the expression of common sub-modules, including the α- and β-globin loci, by directly regulating transcription factors. Our analysis suggests that several genetic variants associated with myeloid blood cell traits alter the activity of a distal enhancer of MYB (~140 kb away), with downstream consequences on hemoglobin genes expression and cell state. These data have implications for the understanding of enhancer-associated traits and emphasize the flexibility of controlling transcriptional systems by modifying enhancer activity.

Xie et al. apply improved strategies for single-cell screens to identify an enhancer-driven transcriptional regulatory network in K562 cells. They demonstrate that the same group of genes can be indirectly regulated by enhancers from distinct genomic loci. These data have implications for the understanding of enhancer-associated traits.

SCINA: A Semi-Supervised Subtyping Algorithm of Single Cells and Bulk Samples

Advances in single-cell RNA sequencing (scRNA-Seq) have allowed for comprehensive analyses of single cell data. However, current analyses of scRNA-Seq data usually start from unsupervised clustering or visualization. These methods ignore prior knowledge of transcriptomes and the probable structures of the data. Moreover, cell identification heavily relies on subjective and possibly inaccurate human inspection afterwards. To address these analytical challenges, we developed SCINA (Semi-supervised Category Identification and Assignment), a semi-supervised model that exploits previously established gene signatures using an expectation-maximization (EM) algorithm. SCINA is applicable to scRNA-Seq and flow cytometry/CyTOF data, as well as other data of similar format. We applied SCINA to a wide range of datasets, and showed its accuracy, stability and efficiency, which exceeded most popular unsupervised approaches. SCINA discovered an intermediate stage of oligodendrocytes from mouse brain scRNA-Seq data. SCINA also detected immune cell population changes in cytometry data in a genetically-engineered mouse model. Furthermore, SCINA performed well with bulk gene expression data. Specifically, we identified a new kidney tumor clade with similarity to FH-deficient tumors (FHD), which we refer to as FHD-like tumors (FHDL). Overall, SCINA provides both methodological advances and biological insights from perspectives different from traditional analytical methods.

Using Gjd3-CreEGFP mice to examine atrioventricular node morphology and composition

The atrioventricular node (AVN) coordinates the timing of atrial and ventricular contraction to optimize cardiac performance. To study this critical function using mouse genetics, however, new reagents are needed that allow AVN-specific manipulation. Here we describe a novel Gjd3-CreEGFP mouse line that successfully recombines floxed alleles within the AVN beginning at E12.5. These mice have been engineered to express CreEGFP under the control of endogenous Gjd3 regulatory elements without perturbing native protein expression. Detailed histological analysis of Gjd3-CreEGFP mice reveals specific labeling of AVN cardiomyocytes and a subset of cardiac endothelial cells. Importantly, we show that Gjd3-CreEGFP mice have preserved cardiac mechanical and electrical function. In one application of our newly described mouse line, we provide a three-dimensional (3D) view of the AVN using tissue clearing combined with confocal microscopy. With this 3D model as a reference, we identify specific AVN sub-structures based on marker staining characteristics. In addition, we use our Gjd3-CreEGFP mice to guide microdissection of the AVN and construction of a single-cell atlas. Thus, our results establish a new transgenic tool for AVN-specific recombination, provide an updated model of AVN morphology, and describe a roadmap for exploring AVN cellular heterogeneity.

Identification of functionally distinct fibro-inflammatory and adipogenic stromal subpopulations in visceral adipose tissue of adult mice

White adipose tissue (WAT) remodeling is dictated by coordinated interactions between adipocytes and resident stromal-vascular cells; however, the functional heterogeneity of adipose stromal cells has remained unresolved. We combined single-cell RNA-sequencing and FACS to identify and isolate functionally distinct subpopulations of PDGFRβ+ stromal cells within visceral WAT of adult mice. LY6C- CD9- PDGFRβ+ cells represent highly adipogenic visceral adipocyte precursor cells (‘APCs’), whereas LY6C+ PDGFRβ+ cells represent fibro-inflammatory progenitors (‘FIPs’). FIPs lack adipogenic capacity, display pro-fibrogenic/pro-inflammatory phenotypes, and can exert an anti-adipogenic effect on APCs. The pro-inflammatory phenotype of PDGFRβ+ cells is regulated, at least in part, by NR4A nuclear receptors. These data highlight the functional heterogeneity of visceral WAT perivascular cells, and provide insight into potential cell-cell interactions impacting adipogenesis and inflammation. These improved strategies to isolate FIPs and APCs from visceral WAT will facilitate the study of physiological WAT remodeling and mechanisms leading to metabolic dysfunction.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Fat tissue, also known as white adipose tissue, specializes in storing excess calories. Much of this storage happens under the skin, but fat tissue can also build up inside the abdomen and surround organs, where it is known as ‘visceral’ fat. When visceral fat tissue is unhealthy, it may help diseases such as diabetes and heart disease to develop.

Unhealthy fat tissue contains enlarged fat cells, which may die from overwork. The stress this places on the surrounding tissue activates the immune system, causing inflammation and the build-up of collagen fibers around the cells (a condition known as fibrosis). Not all people develop this type of unhealthy fat tissue, but we do not yet understand why.

In many tissues, blood vessels serve as a home for several types of adult stem cells that help to rejuvenate the tissue following damage. To identify these cells, Hepler et al. analyzed the genes used by more than 3,000 cells living around the blood vessels in the visceral fat of adult mice. Recent work had already revealed that stem cells called adipocyte precursor cells live in this region. Hepler et al. now reveal the presence of a second group of cells, termed fibro-inflammatory progenitor cells (or FIPs for short).

To investigate the roles of each cell type in more detail, Hepler et al. developed a new technique to isolate the adipocyte precursor cells from other cell types. When grown in the right conditions in petri dishes, the adipocyte precursor cells were able to form new fat cells. They could also make new fat cells when transplanted into mice that lacked fat tissue. By contrast, the FIPs can suppress the activity of adipocyte precursor cells and activate immune cells. They may also help fibrosis to develop.

It is not yet clear whether FIPs are present in human fat tissue. But, if they are, understanding them in greater detail may suggest new ways to treat diabetes and heart disease in obese people.

Frequent sgRNA-barcode recombination in single-cell perturbation assays

Simultaneously detecting CRISPR-based perturbations and induced transcriptional changes in the same cell is a powerful approach to unraveling genome function. Several lentiviral approaches have been developed, some of which rely on the detection of distally located genetic barcodes as an indirect proxy of sgRNA identity. Since barcodes are often several kilobases from their corresponding sgRNAs, viral recombination-mediated swapping of barcodes and sgRNAs is feasible. Using a self-circularization-based sgRNA-barcode library preparation protocol, we estimate the recombination rate to be ~50% and we trace this phenomenon to the pooled viral packaging step. Recombination is random, and decreases the signal-to-noise ratio of the assay. Our results suggest that alternative approaches can increase the throughput and sensitivity of single-cell perturbation assays.

Tet-Assisted Bisulfite Sequencing (TAB-seq)

5-Hydroxymethylcytosine (5hmC) is a modified form of cytosine, which has recently been found in mammalian cells and tissues. 5hmC is derived from 5-methylcytosine (5mC) by Ten-eleven translocation (TET) family protein-mediated oxidation and may regulate gene expression. Numerous affinity-based profiling methods have been developed to help understand the exact function of 5hmC in the genome. However, these methods have a relatively low resolution (~100 bp) without quantitative information of the modification percentage on each site. Here we demonstrated the detailed procedure of Tet-Assistant Bisulfite Sequencing (TAB-Seq), which can detect 5hmC at single-base resolution and quantify its abundance at each site. In this protocol, the genomic DNA is first treated with βGT and recombinant mTet1 consecutively to convert 5hmC to 5gmC and 5mC to 5caC, respectively. The treated genomic DNA can be directly applied to bisulfite treatment to detect 5hmC on specific loci or applied to whole-genome bisulfite sequencing as needed.

In Situ Capture of Chromatin Interactions by Biotinylated dCas9

Cis-regulatory elements (CREs) are commonly recognized by correlative chromatin features, yet the molecular composition of the vast majority of CREs in chromatin remains unknown. Here, we describe a CRISPR affinity purification in situ of regulatory elements (CAPTURE) approach to unbiasedly identify locus-specific chromatin-regulating protein complexes and long-range DNA interactions. Using an in vivo biotinylated nuclease-deficient Cas9 protein and sequence-specific guide RNAs, we show high-resolution and selective isolation of chromatin interactions at a single-copy genomic locus. Purification of human telomeres using CAPTURE identifies known and new telomeric factors. In situ capture of individual constituents of the enhancer cluster controlling human β-globin genes establishes evidence for composition-based hierarchical organization. Furthermore, unbiased analysis of chromatin interactions at disease-associated cis-elements and developmentally regulated super-enhancers reveals spatial features that causally control gene transcription. Thus, comprehensive and unbiased analysis of locus-specific regulatory composition provides mechanistic insight into genome structure and function in development and disease.

Multiplexed Engineering and Analysis of Combinatorial Enhancer Activity in Single Cells

The study of enhancers has been hampered by the scarcity of methods to systematically quantify their endogenous activity. We develop Mosaic-seq to systematically perturb enhancers and measure their endogenous activities at single-cell resolution. Mosaic-seq uses a CRISPR barcoding system to jointly measure a cell's transcriptome and its sgRNA modulators, thus quantifying the effects of dCas9-KRAB-mediated enhancer repression in single cells. Applying Mosaic-seq to 71 constituent enhancers from 15 super-enhancers, our analysis of 51,448 sgRNA-induced transcriptomes finds that only a small number of constituents are major effectors of target gene expression. Binding of p300 and RNAPII are key features of these constituents. We determine two key parameters of enhancer activity in single cells: their penetrance in a population and their contribution to expression in these cells. Through combinatorial interrogation, we find that simultaneous repression of multiple weak constituents can alter super-enhancer activity in a manner greatly exceeding repression of individual constituents.

An LXR-Cholesterol Axis Creates a Metabolic Co-Dependency for Brain Cancers

Small-molecule inhibitors targeting growth factor receptors have failed to show efficacy for brain cancers, potentially due to their inability to achieve sufficient drug levels in the CNS. Targeting non-oncogene tumor co-dependencies provides an alternative approach, particularly if drugs with high brain penetration can be identified. Here we demonstrate that the highly lethal brain cancer glioblastoma (GBM) is remarkably dependent on cholesterol for survival, rendering these tumors sensitive to Liver X receptor (LXR) agonist-dependent cell death. We show that LXR-623, a clinically viable, highly brain-penetrant LXRα-partial/LXRβ-full agonist selectively kills GBM cells in an LXRβ- and cholesterol-dependent fashion, causing tumor regression and prolonged survival in mouse models. Thus, a metabolic co-dependency provides a pharmacological means to kill growth factor-activated cancers in the CNS.

5mC oxidation by Tet2 modulates enhancer activity and timing of transcriptome reprogramming during differentiation

In mammals, cytosine methylation (5mC) is widely distributed throughout the genome but is notably depleted from active promoters and enhancers. While the role of DNA methylation in promoter silencing has been well documented, the function of this epigenetic mark at enhancers remains unclear. Recent experiments have demonstrated that enhancers are enriched for 5-hydroxymethylcytosine (5hmC), an oxidization product of the Tet family of 5mC dioxygenases and an intermediate of DNA demethylation. These results support the involvement of Tet proteins in the regulation of dynamic DNA methylation at enhancers. By mapping DNA methylation and hydroxymethylation at base resolution, we find that deletion of Tet2 causes extensive loss of 5hmC at enhancers, accompanied by enhancer hypermethylation, reduction of enhancer activity, and delayed gene induction in the early steps of differentiation. Our results reveal that DNA demethylation modulates enhancer activity, and its disruption influences the timing of transcriptome reprogramming during cellular differentiation.

N6-methyladenosine-dependent regulation of messenger RNA stability

N(6)-methyladenosine (m(6)A) is the most prevalent internal (non-cap) modification present in the messenger RNA of all higher eukaryotes. Although essential to cell viability and development, the exact role of m(6)A modification remains to be determined. The recent discovery of two m(6)A demethylases in mammalian cells highlighted the importance of m(6)A in basic biological functions and disease. Here we show that m(6)A is selectively recognized by the human YTH domain family 2 (YTHDF2) 'reader' protein to regulate mRNA degradation. We identified over 3,000 cellular RNA targets of YTHDF2, most of which are mRNAs, but which also include non-coding RNAs, with a conserved core motif of G(m(6)A)C. We further establish the role of YTHDF2 in RNA metabolism, showing that binding of YTHDF2 results in the localization of bound mRNA from the translatable pool to mRNA decay sites, such as processing bodies. The carboxy-terminal domain of YTHDF2 selectively binds to m(6)A-containing mRNA, whereas the amino-terminal domain is responsible for the localization of the YTHDF2-mRNA complex to cellular RNA decay sites. Our results indicate that the dynamic m(6)A modification is recognized by selectively binding proteins to affect the translation status and lifetime of mRNA.

Epigenetic memory at embryonic enhancers identified in DNA methylation maps from adult mouse tissues

Mammalian development requires cytosine methylation, a heritable epigenetic mark of cellular memory believed to maintain a cell’s unique gene expression pattern. However, it remains unclear how dynamic DNA methylation relates to cell-type specific gene expression and animal development. Here, by mapping base resolution methylomes in 17 adult mouse tissues at shallow coverage, we identify 302,864 tissue-specific differentially methylated regions (tsDMRs) and estimate that >6.7% of the mouse genome is variably methylated. Supporting a prominent role for DNA methylation in gene regulation, most tsDMRs occur at distal cis-regulatory elements. Surprisingly, some tsDMRs mark enhancers dormant in adult tissues but active in embryonic development. These “vestigial” enhancers are hypomethylated and lack active histone modifications in adult tissue, but nevertheless exhibit activity during embryonic development. Our results provide new insights into the role of DNA methylation at tissue-specific enhancers and suggest that epigenetic memory of embryonic development may be retained in adult tissues.

Abstract 1111: Whole genome bisulfite sequencing reveals tissue-specific DNA methylation in a normal mouse

Epigenetic alteration is a hallmark of cancer. In particular, compared to normal cells, the DNA methylomes of cancer cells often exhibit two contrasting phenomena: 1) hypomethylation of large regions of the genome, notably at regions of late replicating DNA; and 2) hypermethylation of focal sites, frequently at CpG-island containing promoters. While this epigenetic variation between cancer cells and normal counterparts has been well documented, the variation between different normal cell types is still unknown. In an effort to better understand the epigenetic variation of a normal individual, we profiled DNA methylation using whole genome bisulfite sequencing in 16 tissues isolated from an individual mouse. We observed a unique distribution of CpG methylation for each tissue, which cluster based on cell lineage and type. Most somatic tissues exhibit a small but consistent decrease in methylation at late-replicating regions of the genome, suggesting the beginnings of global hypomethylation. Our global analysis identified only one-eighth of the genome as tissue-specific differentially methylated regions. Remarkably, the vast majority of these regions exhibit hallmarks of cis-regulatory activity, including: small size (<500 bp), DNase I hypersensitivity, enhancer chromatin modifications (H3K4me1/H3K27ac), and evolutionary conservation. Our results reveal that most somatic cells are primed for global hypomethylation in regions of the genome that are typically hypomethylated in cancer cells, supporting previous observations that hypomethylation may be an indicator of cell age. In addition, our results in normal tissue suggest that hypomethylation in cancer cells could alter transcriptional regulation through disruption of cis-regulatory activity.

Citation Format: Gary C. Hon, Yin Shen, David F. McCleary, Lee Edsall, Samantha Kuan, Bing Ren. Whole genome bisulfite sequencing reveals tissue-specific DNA methylation in a normal mouse. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1111. doi:10.1158/1538-7445.AM2013-1111

Tet-assisted bisulfite sequencing of 5-hydroxymethylcytosine

A complete understanding of the potential function of 5-hydroxymethylcytosine (5-hmC), a DNA cytosine modification in mammalian cells, requires an accurate single-base resolution sequencing method. Here we describe a modified bisulfite-sequencing method, Tet-assisted bisulfite sequencing (TAB-seq), which can identify 5-hmC at single-base resolution, as well as determine its abundance at each modification site. This protocol involves β-glucosyltransferase (β-GT)-mediated protection of 5-hmC (glucosylation) and recombinant mouse Tet1(mTet1)-mediated oxidation of 5-methylcytosine (5-mC) to 5-carboxylcytosine (5-caC). After the subsequent bisulfite treatment and PCR amplification, both cytosine and 5-caC (derived from 5-mC) are converted to thymine (T), whereas 5-hmC reads as C. The treated genomic DNA is suitable for both whole-genome and locus-specific sequencing. The entire procedure (which does not include data analysis) can be completed in 14 d for whole-genome sequencing or 7 d for locus-specific sequencing.

Base-resolution analysis of 5-hydroxymethylcytosine in the mammalian genome

The study of 5-hydroxylmethylcytosines (5hmC) has been hampered by the lack of a method to map it at single-base resolution on a genome-wide scale. Affinity purification-based methods cannot precisely locate 5hmC nor accurately determine its relative abundance at each modified site. We here present a genome-wide approach, Tet-assisted bisulfite sequencing (TAB-Seq), that when combined with traditional bisulfite sequencing can be used for mapping 5hmC at base resolution and quantifying the relative abundance of 5hmC as well as 5mC. Application of this method to embryonic stem cells not only confirms widespread distribution of 5hmC in the mammalian genome but also reveals sequence bias and strand asymmetry at 5hmC sites. We observe high levels of 5hmC and reciprocally low levels of 5mC near but not on transcription factor-binding sites. Additionally, the relative abundance of 5hmC varies significantly among distinct functional sequence elements, suggesting different mechanisms for 5hmC deposition and maintenance.

CBX3 regulates efficient RNA processing genome-wide

CBX5, CBX1, and CBX3 (HP1α, β, and γ, respectively) play an evolutionarily conserved role in the formation and maintenance of heterochromatin. In addition, CBX5, CBX1, and CBX3 may also participate in transcriptional regulation of genes. Recently, CBX3 binding to the bodies of a subset of genes has been observed in human and murine cells. However, the generality of this phenomenon and the role CBX3 may play in this context are unknown. Genome-wide localization analysis reveals CBX3 binding at genic regions, which strongly correlates with gene activity across multiple cell types. Depletion of CBX3 resulted in down-regulation of a subset of target genes. Loss of CBX3 binding leads to a more dramatic accumulation of unspliced nascent transcripts. In addition, we observed defective recruitment of splicing factors, including SNRNP70, to CBX3 target genes. Collectively, our data suggest a role for CBX3 in aiding in efficient cotranscriptional RNA processing.

Next-generation genomics: an integrative approach

Integrating results from diverse experiments is an essential process in our effort to understand the logic of complex systems, such as development, homeostasis and responses to the environment. With the advent of high-throughput methods--including genome-wide association (GWA) studies, chromatin immunoprecipitation followed by sequencing (ChIP-seq) and RNA sequencing (RNA-seq)--acquisition of genome-scale data has never been easier. Epigenomics, transcriptomics, proteomics and genomics each provide an insightful, and yet one-dimensional, view of genome function; integrative analysis promises a unified, global view. However, the large amount of information and diverse technology platforms pose multiple challenges for data access and processing. This Review discusses emerging issues and strategies related to data integration in the era of next-generation genomics.

Predictive chromatin signatures in the mammalian genome

The DNA sequence of an organism is a blueprint of life: it harbors not only the information about proteins and other molecules produced in each cell, but also instructions on when and where such molecules are made. Chromatin, the structure of histone and DNA that has co-evolved with eukaryotic genome, also contains information that indicates the function and activity of the underlying DNA sequences. Such information exists in the form of covalent modifications to the histone proteins that comprise the nucleosome. Thanks to the development of high throughput technologies such as DNA microarrays and next generation DNA sequencing, we have begun to associate the various combinations of chromatin modification patterns with functional sequences in the human genome. Here, we review the rapid progress from descriptive observations of histone modification profiles to highly predictive models enabling use of chromatin signatures to enumerate novel functional sequences in mammalian genomes that have escaped previous detection.

Comprehensive curation and analysis of global interaction networks in Saccharomyces cerevisiae

BACKGROUND

The study of complex biological networks and prediction of gene function has been enabled by high-throughput (HTP) methods for detection of genetic and protein interactions. Sparse coverage in HTP datasets may, however, distort network properties and confound predictions. Although a vast number of well substantiated interactions are recorded in the scientific literature, these data have not yet been distilled into networks that enable system-level inference.

RESULTS

We describe here a comprehensive database of genetic and protein interactions, and associated experimental evidence, for the budding yeast Saccharomyces cerevisiae, as manually curated from over 31,793 abstracts and online publications. This literature-curated (LC) dataset contains 33,311 interactions, on the order of all extant HTP datasets combined. Surprisingly, HTP protein-interaction datasets currently achieve only around 14% coverage of the interactions in the literature. The LC network nevertheless shares attributes with HTP networks, including scale-free connectivity and correlations between interactions, abundance, localization, and expression. We find that essential genes or proteins are enriched for interactions with other essential genes or proteins, suggesting that the global network may be functionally unified. This interconnectivity is supported by a substantial overlap of protein and genetic interactions in the LC dataset. We show that the LC dataset considerably improves the predictive power of network-analysis approaches. The full LC dataset is available at the BioGRID (http://www.thebiogrid.org) and SGD (http://www.yeastgenome.org/) databases.

CONCLUSION

Comprehensive datasets of biological interactions derived from the primary literature provide critical benchmarks for HTP methods, augment functional prediction, and reveal system-level attributes of biological networks.