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.

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.

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.

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.