Regulation of transcription through acetylation of H3K122 on the lateral surface of the histone octamer.

Histone modifications are key regulators of chromatin function. However, little is known to what extent histone modifications can directly impact on chromatin. Here, we address how a modification within the globular domain of histones regulates chromatin function. We demonstrate that H3K122ac can be sufficient to stimulate transcription and that mutation of H3K122 impairs transcriptional activation, which we attribute to a direct effect of H3K122ac on histone-DNA binding. In line with this, we find that H3K122ac defines genome-wide genetic elements and chromatin features associated with active transcription. Furthermore, H3K122ac is catalyzed by the coactivators p300/CBP and can be induced by nuclear hormone receptor signaling. Collectively, this suggests that transcriptional regulators elicit their effects not only via signaling to histone tails but also via direct structural perturbation of nucleosomes by directing acetylation to their lateral surface.

Hedgehog-FGF signaling axis patterns anterior mesoderm during gastrulation.

The mechanisms used by embryos to pattern tissues across their axes has fascinated developmental biologists since the founding of embryology. Here, using single-cell technology, we interrogate complex patterning defects and define a Hedgehog (Hh)-fibroblast growth factor (FGF) signaling axis required for anterior mesoderm lineage development during gastrulation. Single-cell transcriptome analysis of Hh-deficient mesoderm revealed selective deficits in anterior mesoderm populations, culminating in defects to anterior embryonic structures, including the pharyngeal arches, heart, and anterior somites. Transcriptional profiling of Hh-deficient mesoderm during gastrulation revealed disruptions to both transcriptional patterning of the mesoderm and FGF signaling for mesoderm migration. Mesoderm-specific Fgf4/Fgf8 double-mutants recapitulated anterior mesoderm defects and Hh-dependent GLI transcription factors modulated enhancers at FGF gene loci. Cellular migration defects during gastrulation induced by Hh pathway antagonism were mitigated by the addition of FGF4 protein. These findings implicate a multicomponent signaling hierarchy activated by Hh ligands from the embryonic node and executed by FGF signals in nascent mesoderm to control anterior mesoderm patterning.

Transcription factor Ebf1 regulates differentiation stage-specific signaling, proliferation, and survival of B cells.

The transcription factor Ebf1 is an important determinant of early B lymphopoiesis. To gain insight into the functions of Ebf1 at distinct stages of differentiation, we conditionally inactivated Ebf1. We found that Ebf1 is required for the proliferation, survival, and signaling of pro-B cells and peripheral B-cell subsets, including B1 cells and marginal zone B cells. The proliferation defect of Ebf1-deficient pro-B cells and the impaired expression of multiple cell cycle regulators are overcome by transformation with v-Abl. The survival defect of transformed Ebf1(fl/fl) pro-B cells can be rescued by the forced expression of the Ebf1 targets c-Myb or Bcl-x(L). In mature B cells, Ebf1 deficiency interferes with signaling via the B-cell-activating factor receptor (BAFF-R)- and B-cell receptor (BCR)-dependent Akt pathways. Moreover, Ebf1 is required for germinal center formation and class switch recombination. Genome-wide analyses of Ebf1-mediated gene expression and chromatin binding indicate that Ebf1 regulates both common and distinct sets of genes in early and late stage B cells. By regulating important components of transcription factor and signaling networks, Ebf1 appears to be involved in the coordination of cell proliferation, survival, and differentiation at multiple stages of B lymphopoiesis.

Early B cell factor 1 regulates B cell gene networks by activation, repression, and transcription- independent poising of chromatin.

The transcription factor early B cell factor-1 (Ebf1) is a key determinant of B lineage specification and differentiation. To gain insight into the molecular basis of Ebf1 function in early-stage B cells, we combined a genome-wide ChIP sequencing analysis with gain- and loss-of-function transcriptome analyses. Among 565 genes that are occupied and transcriptionally regulated by Ebf1, we identified large sets involved in (pre)-B cell receptor and Akt signaling, cell adhesion, and migration. Interestingly, a third of previously described Pax5 targets was found to be occupied by Ebf1. In addition to Ebf1-activated and -repressed genes, we identified targets at which Ebf1 induces chromatin changes that poise the genes for expression at subsequent stages of differentiation. Poised chromatin states on specific targets could also be established by Ebf1 expression in T cells but not in NIH 3T3 cells, suggesting that Ebf1 acts as a "pioneer" factor in a hematopoietic chromatin context.

Tissue- and strain-specific effects of a genotoxic carcinogen 1,3-butadiene on chromatin and transcription.

Epigenetic effects of environmental chemicals are under intense investigation to fill existing knowledge gaps between environmental/occupational exposures and adverse health outcomes. Chromatin accessibility is one prominent mechanism of epigenetic control of transcription, and understanding of the chemical effects on both could inform the causal role of epigenetic alterations in disease mechanisms. In this study, we hypothesized that baseline variability in chromatin organization and transcription profiles among various tissues and mouse strains influence the outcome of exposure to the DNA damaging chemical 1,3-butadiene. To test this hypothesis, we evaluated DNA damage along with comprehensive quantification of RNA transcripts (RNA-seq), identification of accessible chromatin (ATAC-seq), and characterization of regions with histone modifications associated with active transcription (ChIP-seq for acetylation at histone 3 lysine 27, H3K27ac). We collected these data in the lung, liver, and kidney of mice from two genetically divergent strains, C57BL/6J and CAST/EiJ, that were exposed to clean air or to 1,3-butadiene (~600 ppm) for 2 weeks. We found that tissue effects dominate differences in both gene expression and chromatin states, followed by strain effects. At baseline, xenobiotic metabolism was consistently more active in CAST/EiJ, while immune system pathways were more active in C57BL/6J across tissues. Surprisingly, even though all three tissues in both strains harbored butadiene-induced DNA damage, little transcriptional effect of butadiene was observed in liver and kidney. Toxicologically relevant effects of butadiene in the lung were on the pathways of xenobiotic metabolism and inflammation. We also found that variability in chromatin accessibility across individuals (i.e., strains) only partially explains the variability in transcription. This study showed that variation in the basal states of epigenome and transcriptome may be useful indicators for individuals or tissues susceptible to genotoxic environmental chemicals.

Bile acid fitness determinants of a Bacteroides fragilis isolate from a human pouchitis patient.

IMPORTANCE: The Gram-negative bacterium Bacteroides fragilis is a common member of the human gut microbiota that colonizes multiple host niches and can influence human physiology through a variety of mechanisms. Identification of genes that enable B. fragilis to grow across a range of host environments has been impeded in part by the relatively limited genetic tractability of this species. We have developed a high-throughput genetic resource for a B. fragilis strain isolated from a UC pouchitis patient. Bile acids limit microbial growth and are altered in abundance in UC pouches, where B. fragilis often blooms. Using this resource, we uncovered pathways and processes that impact B. fragilis fitness in bile and that may contribute to population expansions during bouts of gut inflammation.

Simultaneous Measurement of DNA Methylation and Nucleosome Occupancy in Single Cells Using scNOMe-Seq.

Single-cell Nucleosome Occupancy and Methylome sequencing (scNOMe-seq) is a multimodal assay that simultaneously measures endogenous DNA methylation and nucleosome occupancy (i.e., chromatin accessibility) in single cells. scNOMe-seq combines the activity of a GpC Methyltransferase, an enzyme which methylates cytosines in GpC dinucleotides, with bisulfite conversion, whereby unmethylated cytosines are converted into thymines. Because GpC Methyltransferase acts only on cytosines present in non-nucleosomal regions of the genome, the subsequent bisulfite conversion step not only detects the endogenous DNA methylation, but also reveals the genome-wide pattern of chromatin accessibility. Implementing this technology at the single-cell level helps to capture the dynamics governing methylation and accessibility vary across individual cells and cell types. Here, we provide a scalable plate-based protocol for preparing scNOMe-seq libraries from single nucleus suspensions.

GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership.

Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.

GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership.

Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.

Single-cell genomics improves the discovery of risk variants and genes of atrial fibrillation.

Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 46 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits.

Bile acid fitness determinants of a Bacteroides fragilis isolate from a human pouchitis patient.

Bacteroides fragilis comprises 1-5% of the gut microbiota in healthy humans but can expand to >50% of the population in ulcerative colitis (UC) patients experiencing inflammation. The mechanisms underlying such microbial blooms are poorly understood, but the gut of UC patients has physicochemical features that differ from healthy patients and likely impact microbial physiology. For example, levels of the secondary bile acid deoxycholate (DC) are highly reduced in the ileoanal J-pouch of UC colectomy patients. We isolated a B. fragilis strain from a UC patient with pouch inflammation (i.e. pouchitis) and developed it as a genetic model system to identify genes and pathways that are regulated by DC and that impact B. fragilis fitness in DC and crude bile. Treatment of B. fragilis with a physiologically relevant concentration of DC reduced cell growth and remodeled transcription of one-quarter of the genome. DC strongly induced expression of chaperones and select transcriptional regulators and efflux systems and downregulated protein synthesis genes. Using a barcoded collection of ≈50,000 unique insertional mutants, we further defined B. fragilis genes that contribute to fitness in media containing DC or crude bile. Genes impacting cell envelope functions including cardiolipin synthesis, cell surface glycosylation, and systems implicated in sodium-dependent bioenergetics were major bile acid fitness factors. As expected, there was limited overlap between transcriptionally regulated genes and genes that impacted fitness in bile when disrupted. Our study provides a genome-scale view of a B. fragilis bile response and genetic determinants of its fitness in DC and crude bile.

Ribonucleotide reductase regulatory subunit M2 drives glioblastoma TMZ resistance through modulation of dNTP production.

During therapy, adaptations driven by cellular plasticity are partly responsible for driving the inevitable recurrence of glioblastoma (GBM). To investigate plasticity-induced adaptation during standard-of-care chemotherapy temozolomide (TMZ), we performed in vivo single-cell RNA sequencing in patient-derived xenograft (PDX) tumors of GBM before, during, and after therapy. Comparing single-cell transcriptomic patterns identified distinct cellular populations present during TMZ therapy. Of interest was the increased expression of ribonucleotide reductase regulatory subunit M2 (RRM2), which we found to regulate dGTP and dCTP production vital for DNA damage response during TMZ therapy. Furthermore, multidimensional modeling of spatially resolved transcriptomic and metabolomic analysis in patients' tissues revealed strong correlations between RRM2 and dGTP. This supports our data that RRM2 regulates the demand for specific dNTPs during therapy. In addition, treatment with the RRM2 inhibitor 3-AP (Triapine) enhances the efficacy of TMZ therapy in PDX models. We present a previously unidentified understanding of chemoresistance through critical RRM2-mediated nucleotide production.

Multiomic analysis reveals cellular and epigenetic plasticity in intestinal pouches of ulcerative colitis patients.

Background & Aims: Total proctocolectomy with ileal pouch anal anastomosis (IPAA) is the standard of care for patients with severe treatment resistant ulcerative colitis (UC). Despite improvements in patient outcomes, about 50% of patients will develop inflammation of the pouch within 1-2 years following surgery. Establishment of UC pouches is associated with profound histological changes of the mucosa. A detailed characterization of these changes on a cellular and molecular level is crucial for an improved understanding of pouch physiology and diseases management. Methods: We generated cell-type-resolved transcriptional and epigenetic atlases of UC pouches using scRNA-seq and scATAC-seq data from paired biopsy samples from the ileal pouch and ileal segment above the pouch (pre-pouch) of UC-IPAA patients (n=6, female=2) without symptoms. We also collected data from paired biopsies of the terminal ileum (TI) and ascending colon (AC) from healthy controls (n=6, female=3). Results: We identified novel populations of colon-like absorptive and secretory epithelial cells, constituting a significant proportion of the epithelial cell fraction in the pouch but not in matched pre-pouch samples. Pouch-specific enterocytes expressed colon-specific genes, including CEACAM5, CA2. However, in contrast to normal colonic epithelium, these cells also expressed a range of inflammatory and secretory genes, similar to previously detected gene expression signatures in IBD patients. Comparison to longitudinal bulk RNA-seq data from UC pouches demonstrated that colon-like epithelial cells are present early after pouch functionalization and independently of subsequent pouchitis. Finally, single cell chromatin accessibility revealed activation colonic transcriptional regulators, including CDX1, NFIA, and EHF. Conclusion: UC pouches are characterized by partial colonic metaplasia of the epithelium. These data constitute a resource of transcriptomic and epigenetic signatures of cell populations in the pouch and provide an anchor for understanding the underlying molecular mechanisms of pouchitis.

A Roadmap for the Human Gut Cell Atlas.

The number of studies investigating the human gastrointestinal tract using various single-cell profiling methods has increased substantially in the past few years. Although this increase provides a unique opportunity for the generation of the first comprehensive Human Gut Cell Atlas (HGCA), there remains a range of major challenges ahead. Above all, the ultimate success will largely depend on a structured and coordinated approach that aligns global efforts undertaken by a large number of research groups. In this Roadmap, we discuss a comprehensive forward-thinking direction for the generation of the HGCA on behalf of the Gut Biological Network of the Human Cell Atlas. Based on the consensus opinion of experts from across the globe, we outline the main requirements for the first complete HGCA by summarizing existing data sets and highlighting anatomical regions and/or tissues with limited coverage. We provide recommendations for future studies and discuss key methodologies and the importance of integrating the healthy gut atlas with related diseases and gut organoids. Importantly, we critically overview the computational tools available and provide recommendations to overcome key challenges.

Hedgehog signaling activates a mammalian heterochronic gene regulatory network controlling differentiation timing across lineages.

Many developmental signaling pathways have been implicated in lineage-specific differentiation; however, mechanisms that explicitly control differentiation timing remain poorly defined in mammals. We report that murine Hedgehog signaling is a heterochronic pathway that determines the timing of progenitor differentiation. Hedgehog activity was necessary to prevent premature differentiation of second heart field (SHF) cardiac progenitors in mouse embryos, and the Hedgehog transcription factor GLI1 was sufficient to delay differentiation of cardiac progenitors in vitro. GLI1 directly activated a de novo progenitor-specific network in vitro, akin to that of SHF progenitors in vivo, which prevented the onset of the cardiac differentiation program. A Hedgehog signaling-dependent active-to-repressive GLI transition functioned as a differentiation timer, restricting the progenitor network to the SHF. GLI1 expression was associated with progenitor status across germ layers, and it delayed the differentiation of neural progenitors in vitro, suggesting a broad role for Hedgehog signaling as a heterochronic pathway.

Divergence in alternative polyadenylation contributes to gene regulatory differences between humans and chimpanzees.

While comparative functional genomic studies have shown that inter-species differences in gene expression can be explained by corresponding inter-species differences in genetic and epigenetic regulatory mechanisms, co-transcriptional mechanisms, such as alternative polyadenylation (APA), have received little attention. We characterized APA in lymphoblastoid cell lines from six humans and six chimpanzees by identifying and estimating the usage for 44,432 polyadenylation sites (PAS) in 9518 genes. Although APA is largely conserved, 1705 genes showed significantly different PAS usage (FDR 0.05) between species. Genes with divergent APA also tend to be differentially expressed, are enriched among genes showing differences in protein translation, and can explain a subset of observed inter-species protein expression differences that do not differ at the transcript level. Finally, we found that genes with a dominant PAS, which is used more often than other PAS, are particularly enriched for differentially expressed genes.

Alternative polyadenylation mediates genetic regulation of gene expression.

Little is known about co-transcriptional or post-transcriptional regulatory mechanisms linking noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3' Seq to study the impact of genetic variation on alternative polyadenylation (APA) in the nuclear and total mRNA fractions of 52 HapMap Yoruba human lymphoblastoid cell lines. We mapped 602 APA quantitative trait loci (apaQTLs) at 10% FDR, of which 152 were nuclear specific. Effect sizes at intronic apaQTLs are negatively correlated with eQTL effect sizes. These observations suggest genetic variants can decrease mRNA expression levels by increasing usage of intronic PAS. We also identified 24 apaQTLs associated with protein levels, but not mRNA expression. Finally, we found that 19% of apaQTLs can be associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression, protein expression, and disease risk, and reveals uncharted modes of genetic regulation.

Systematic Comparison of High-throughput Single-Cell and Single-Nucleus Transcriptomes during Cardiomyocyte Differentiation.

A comprehensive reference map of all cell types in the human body is necessary for improving our understanding of fundamental biological processes and in diagnosing and treating disease. High-throughput single-cell RNA sequencing techniques have emerged as powerful tools to identify and characterize cell types in complex and heterogeneous tissues. However, extracting intact cells from tissues and organs is often technically challenging or impossible, for example in heart or brain tissue. Single-nucleus RNA sequencing provides an alternative way to obtain transcriptome profiles of such tissues. To systematically assess the differences between high-throughput single-cell and single-nuclei RNA-seq approaches, we compared Drop-seq and DroNc-seq, two microfluidic-based 3' RNA capture technologies that profile total cellular and nuclear RNA, respectively, during a time course experiment of human induced pluripotent stem cells (iPSCs) differentiating into cardiomyocytes. Clustering of time-series transcriptomes from Drop-seq and DroNc-seq revealed six distinct cell types, five of which were found in both techniques. Furthermore, single-cell trajectories reconstructed from both techniques reproduced expected differentiation dynamics. We then applied DroNc-seq to postmortem heart tissue to test its performance on heterogeneous human tissue samples. Our data confirm that DroNc-seq yields similar results to Drop-seq on matched samples and can be successfully used to generate reference maps for the human cell atlas.

Gut flora, Toll-like receptors and nuclear receptors: a tripartite communication that tunes innate immunity in large intestine.

Separating the large intestine from gut flora is a robust layer of epithelial cells. This barrier is armed with an array of recognizing receptors that collectively set the host innate response. Here, we use nuclear receptors (NRs) and Toll-like receptors (TLRs), suggested to act as second messengers in the communication between microorganisms and epithelial cells, as probes to assess the impact of gut flora on innate immunity in germ-free (GF) mice. Using quantitative real-time polymerase chain reaction analyses, we show that 37/49 NRs are expressed in colonic cells of GF mice. Of these, 5 can be modulated by resident flora: LXRalpha, RORgamma and CAR show reduced expression and Nur77 and GCNF display elevated expression in conventionally raised mice compared with GF. Moreover, increased expression levels of TLR-2 and TLR-5 are observed in specific pathogen-free (SPF) mice compared with GF mice, and CAR expression is connected to the TLR-2 signalling pathway. Infections of GF or SPF mice with Yersinia pseudotuberculosis, show that GF intestinal epithelial cells fail to respond, except for CAR, which is downregulated. In contrast, SPF epithelial cells show a downregulation of all the NRs except CAR, which appears to be unaffected. Our findings indicate that gut flora contributes to the development of an intact barrier function.

Two Pif1 Family DNA Helicases Cooperate in Centromere Replication and Segregation in Saccharomyces cerevisiae.

Pif1 family helicases are found in virtually all eukaryotes. Saccharomyces cerevisiae (Sc) encodes two Pif1 family helicases, ScPif1 and Rrm3 ScPif1 is multifunctional, required not only for maintenance of mitochondrial DNA but also for multiple distinct nuclear functions. Rrm3 moves with the replication fork and promotes movement of the fork through ∼1400 hard-to-replicate sites, including centromeres. Here we show that ScPif1, like Rrm3, bound robustly to yeast centromeres but only if the centromere was active. While Rrm3 binding to centromeres occurred in early to mid S phase, about the same time as centromere replication, ScPif1 binding occurred later in the cell cycle when replication of most centromeres is complete. However, the timing of Rrm3 and ScPif1 centromere binding was altered by the absence of the other helicase, such that Rrm3 centromere binding occurred later in pif1-m2 cells and ScPif1 centromere binding occurred earlier in rrm3Δ cells. As shown previously, the modest pausing of replication forks at centromeres seen in wild-type cells was increased in the absence of Rrm3 While a lack of ScPif1 did not result in increased fork pausing at centromeres, pausing was even higher in rrm3Δ pif1Δ cells than in rrm3Δ cells. Likewise, centromere function as monitored by the loss rate of a centromere plasmid was increased in rrm3Δ but not pif1Δ cells, and was even higher in rrm3Δ pif1Δ cells than in rrm3Δ cells. Thus, ScPif1 promotes centromere replication and segregation, but only in the absence of Rrm3 These data also hint at a potential post-S phase function for ScPif1 at centromeres. These studies add to the growing list of ScPif1 functions that promote chromosome stability.