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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Pluripotency reprogramming by competent and incompetent POU factors uncovers temporal dependency for Oct4 and Sox2.

Oct4, along with Sox2 and Klf4 (SK), can induce pluripotency but structurally similar factors like Oct6 cannot. To decode why Oct4 has this unique ability, we compare Oct4-binding, accessibility patterns and transcriptional waves with Oct6 and an Oct4 mutant defective in the dimerization with Sox2 (Oct4defSox2). We find that initial silencing of the somatic program proceeds indistinguishably with or without Oct4. Oct6 mitigates the mesenchymal-to-epithelial transition and derails reprogramming. These effects are a consequence of differences in genome-wide binding, as the early binding profile of Oct4defSox2 resembles Oct4, whilst Oct6 does not bind pluripotency enhancers. Nevertheless, in the Oct6-SK condition many otherwise Oct4-bound locations become accessible but chromatin opening is compromised when Oct4defSox2 occupies these sites. We find that Sox2 predominantly facilitates chromatin opening, whilst Oct4 serves an accessory role. Formation of Oct4/Sox2 heterodimers is essential for pluripotency establishment; however, reliance on Oct4/Sox2 heterodimers declines during pluripotency maintenance.

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.

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.

Variation in DNA-Damage Responses to an Inhalational Carcinogen (1,3-Butadiene) in Relation to Strain-Specific Differences in Chromatin Accessibility and Gene Transcription Profiles in C57BL/6J and CAST/EiJ Mice.

BACKGROUND: The damaging effects of exposure to environmental toxicants differentially affect genetically distinct individuals, but the mechanisms contributing to these differences are poorly understood. Genetic variation affects the establishment of the gene regulatory landscape and thus gene expression, and we hypothesized that this contributes to the observed heterogeneity in individual responses to exogenous cellular insults. OBJECTIVES: We performed an in vivo study of how genetic variation and chromatin organization may dictate susceptibility to DNA damage, and influence the cellular response to such damage, caused by an environmental toxicant. MATERIALS AND METHODS: We measured DNA damage, messenger RNA (mRNA) and microRNA (miRNA) expression, and genome-wide chromatin accessibility in lung tissue from two genetically divergent inbred mouse strains, C57BL/6J and CAST/EiJ, both in unexposed mice and in mice exposed to a model DNA-damaging chemical, 1,3-butadiene. RESULTS: Our results showed that unexposed CAST/EiJ and C57BL/6J mice have very different chromatin organization and transcription profiles in the lung. Importantly, in unexposed CAST/EiJ mice, which acquired relatively less 1,3-butadiene-induced DNA damage, we observed increased transcription and a more accessible chromatin landscape around genes involved in detoxification pathways. Upon chemical exposure, chromatin was significantly remodeled in the lung of C57BL/6J mice, a strain that acquired higher levels of 1,3-butadiene-induced DNA damage, around the same genes, ultimately resembling the molecular profile of CAST/EiJ. CONCLUSIONS: These results suggest that strain-specific changes in chromatin and transcription in response to chemical exposure lead to a "compensation" for underlying genetic-driven interindividual differences in the baseline chromatin and transcriptional state. This work represents an example of how chemical and environmental exposures can be evaluated to better understand gene-by-environment interactions, and it demonstrates the important role of chromatin response in transcriptomic changes and, potentially, in deleterious effects of exposure. https://doi.org/10.1289/EHP1937.

Simultaneous measurement of chromatin accessibility, DNA methylation, and nucleosome phasing in single cells.

Gaining insights into the regulatory mechanisms that underlie the transcriptional variation observed between individual cells necessitates the development of methods that measure chromatin organization in single cells. Here I adapted Nucleosome Occupancy and Methylome-sequencing (NOMe-seq) to measure chromatin accessibility and endogenous DNA methylation in single cells (scNOMe-seq). scNOMe-seq recovered characteristic accessibility and DNA methylation patterns at DNase hypersensitive sites (DHSs). An advantage of scNOMe-seq is that sequencing reads are sampled independently of the accessibility measurement. scNOMe-seq therefore controlled for fragment loss, which enabled direct estimation of the fraction of accessible DHSs within individual cells. In addition, scNOMe-seq provided high resolution of chromatin accessibility within individual loci which was exploited to detect footprints of CTCF binding events and to estimate the average nucleosome phasing distances in single cells. scNOMe-seq is therefore well-suited to characterize the chromatin organization of single cells in heterogeneous cellular mixtures.

PIF1 family DNA helicases suppress R-loop mediated genome instability at tRNA genes.

Saccharomyces cerevisiae encodes two Pif1 family DNA helicases, Pif1 and Rrm3. Rrm3 promotes DNA replication past stable protein complexes at tRNA genes (tDNAs). We identify a new role for the Pif1 helicase: promotion of replication and suppression of DNA damage at tDNAs. Pif1 binds multiple tDNAs, and this binding is higher in rrm3Δ cells. Accumulation of replication intermediates and DNA damage at tDNAs is higher in pif1Δ rrm3Δ than in rrm3Δ cells. DNA damage at tDNAs in the absence of these helicases is suppressed by destabilizing R-loops while Pif1 and Rrm3 binding to tDNAs is increased upon R-loop stabilization. We propose that Rrm3 and Pif1 promote genome stability at tDNAs by displacing the stable multi-protein transcription complex and by removing R-loops. Thus, we identify tDNAs as a new source of R-loop-mediated DNA damage. Given their large number and high transcription rate, tDNAs may be a potent source of genome instability.