Extracellular cytochrome nanowires appear to be ubiquitous in prokaryotes.

Electrically conductive appendages from the anaerobic bacterium Geobacter sulfurreducens, recently identified as extracellular cytochrome nanowires (ECNs), have received wide attention due to numerous potential applications. However, whether other organisms employ similar ECNs for electron transfer remains unknown. Here, using cryoelectron microscopy, we describe the atomic structures of two ECNs from two major orders of hyperthermophilic archaea present in deep-sea hydrothermal vents and terrestrial hot springs. Homologs of Archaeoglobus veneficus ECN are widespread among mesophilic methane-oxidizing Methanoperedenaceae, alkane-degrading Syntrophoarchaeales archaea, and in the recently described megaplasmids called Borgs. The ECN protein subunits lack similarities in their folds; however, they share a common heme arrangement, suggesting an evolutionarily optimized heme packing for efficient electron transfer. The detection of ECNs in archaea suggests that filaments containing closely stacked hemes may be a common and widespread mechanism for long-range electron transfer in both prokaryotic domains of life.

The transcriptional cofactor Tle3 reciprocally controls effector and central memory CD8+ T cell fates.

Antigen-experienced CD8+ T cells form effector and central memory T cells (TEM and TCM cells, respectively); however, the mechanism(s) controlling their lineage plasticity remains incompletely understood. Here we show that the transcription cofactor Tle3 critically regulates TEM and TCM cell fates and lineage stability through dynamic redistribution in antigen-responding CD8+ T cell genome. Genetic ablation of Tle3 promoted CD8+ TCM cell formation at the expense of CD8+ TEM cells. Lineage tracing showed that Tle3-deficient CD8+ TEM cells underwent accelerated conversion into CD8+ TCM cells while retaining robust recall capacity. Tle3 acted as a coactivator for Tbet to increase chromatin opening at CD8+ TEM cell-characteristic sites and to activate CD8+ TEM cell signature gene transcription, while engaging Runx3 and Tcf1 to limit CD8+ TCM cell-characteristic molecular features. Thus, Tle3 integrated functions of multiple transcription factors to guard lineage fidelity of CD8+ TEM cells, and manipulation of Tle3 activity could favor CD8+ TCM cell production.

Tcf1 preprograms the mobilization of glycolysis in central memory CD8+ T cells during recall responses.

The mechanisms underlying the heightened protection mediated by central memory CD8+ T (TCM) cells remain unclear. Here we show that the transcription factor Tcf1 was required in resting TCM cells to generate secondary effector CD8+ T cells and to clear pathogens during recall responses. Recall stimulation of CD8+ TCM cells caused extensive reprogramming of the transcriptome and chromatin accessibility, leading to rapid induction of glycolytic enzymes, cell cycle regulators and transcriptional regulators, including Id3. This cluster of genes did not require Tcf1 in resting CD8+ TCM cells, but depended on Tcf1 for optimal induction and chromatin opening in recall-stimulated CD8+ TCM cells. Tcf1 bound extensively to these recall-induced gene loci in resting CD8+ TCM cells and mediated chromatin interactions that positioned these genes in architectural proximity with poised enhancers. Thus, Tcf1 preprogramed a transcriptional program that supported the bioenergetic and proliferative needs of CD8+ TCM cells in case of a secondary challenge.

Inherited DNA methylation primes the establishment of accessible chromatin during genome activation.

For animals, epigenetic modifications can be globally or partially inherited from gametes after fertilization, and such information is required for proper transcriptional regulation, especially during the process of zygotic genome activation (ZGA). However, the mechanism underlying how the inherited epigenetic signatures affect transcriptional regulation during ZGA remains poorly understood. Here, we performed genome-wide profiling of chromatin accessibility during zebrafish ZGA, which is closely related to zygotic transcriptional regulation. We observed a clear trend toward a gradual increase in accessible chromatin during ZGA. Furthermore, accessible chromatin at the promoters displayed a sequential priority of emergence, and the locations of the accessible chromatin were precisely primed by the enrichment of unmethylated CpGs that were fully inherited from gametes. On the other hand, distal regions with high methylation levels that were inherited from the sperm facilitated the binding of DNA methylation-preferred transcription factors, such as Pou5f3 and Nanog, which contributed to the establishment of accessible chromatin at these loci. Our results demonstrate a model whereby inherited DNA methylation signatures from gametes prime the establishment of accessible chromatin during zebrafish ZGA through two distinct mechanisms.

Comparative analysis of metazoan chromatin organization.

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal 'arms', and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function.

Refined DNase-seq protocol and data analysis reveals intrinsic bias in transcription factor footprint identification.

Sequencing of DNase I hypersensitive sites (DNase-seq) is a powerful technique for identifying cis-regulatory elements across the genome. We studied the key experimental parameters to optimize performance of DNase-seq. Sequencing short fragments of 50-100 base pairs (bp) that accumulate in long internucleosome linker regions was more efficient for identifying transcription factor binding sites compared to sequencing longer fragments. We also assessed the potential of DNase-seq to predict transcription factor occupancy via generation of nucleotide-resolution transcription factor footprints. In modeling the sequence-specific DNase I cutting bias, we found a strong effect that varied over more than two orders of magnitude. This indicates that the nucleotide-resolution cleavage patterns at many transcription factor binding sites are derived from intrinsic DNase I cleavage bias rather than from specific protein-DNA interactions. In contrast, quantitative comparison of DNase I hypersensitivity between states can predict transcription factor occupancy associated with particular biological perturbations.

Dr.seq: a quality control and analysis pipeline for droplet sequencing.

MOTIVATION: Drop-seq has recently emerged as a powerful technology to analyze gene expression from thousands of individual cells simultaneously. Currently, Drop-seq technology requires refinement and quality control (QC) steps are critical for such data analysis. There is a strong need for a convenient and comprehensive approach to obtain dedicated QC and to determine the relationships between cells for ultra-high-dimensional datasets. RESULTS: We developed Dr.seq, a QC and analysis pipeline for Drop-seq data. By applying this pipeline, Dr.seq provides four groups of QC measurements for given Drop-seq data, including reads level, bulk-cell level, individual-cell level and cell-clustering level QC. We assessed Dr.seq on simulated and published Drop-seq data. Both assessments exhibit reliable results. Overall, Dr.seq is a comprehensive QC and analysis pipeline designed for Drop-seq data that is easily extended to other droplet-based data types. AVAILABILITY AND IMPLEMENTATION: Dr.seq is freely available at: http://www.tongji.edu.cn/∼zhanglab/drseq and https://bitbucket.org/tarela/drseq CONTACT: yzhang@tongji.edu.cn SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CRISPR-DO for genome-wide CRISPR design and optimization.

MOTIVATION: Despite the growing popularity in using CRISPR/Cas9 technology for genome editing and gene knockout, its performance still relies on well-designed single guide RNAs (sgRNA). In this study, we propose a web application for the Design and Optimization (CRISPR-DO) of guide sequences that target both coding and non-coding regions in spCas9 CRISPR system across human, mouse, zebrafish, fly and worm genomes. CRISPR-DO uses a computational sequence model to predict sgRNA efficiency, and employs a specificity scoring function to evaluate the potential of off-target effect. It also provides information on functional conservation of target sequences, as well as the overlaps with exons, putative regulatory sequences and single-nucleotide polymorphisms (SNPs). The web application has a user-friendly genome-browser interface to facilitate the selection of the best target DNA sequences for experimental design. AVAILABILITY AND IMPLEMENTATION: CRISPR-DO is available at http://cistrome.org/crispr/ CONTACT: qiliu@tongji.edu.cn or hanxu@jimmy.harvard.edu or xsliu@jimmy.harvard.eduSupplementary information: Supplementary data are available at Bioinformatics online.

CistromeMap: a knowledgebase and web server for ChIP-Seq and DNase-Seq studies in mouse and human.

SUMMARY: Transcription and chromatin regulators, and histone modifications play essential roles in gene expression regulation. We have created CistromeMap as a web server to provide a comprehensive knowledgebase of all of the publicly available ChIP-Seq and DNase-Seq data in mouse and human. We have also manually curated metadata to ensure annotation consistency, and developed a user-friendly display matrix for quick navigation and retrieval of data for specific factors, cells and papers. Finally, we provide users with summary statistics of ChIP-Seq and DNase-Seq studies.

Single cell RNA sequencing unravels mechanisms underlying senescence-like phenotypes of alveolar macrophages.

Alveolar macrophages (AMs) are resident innate immune cells that play vital roles in maintaining lung physiological functions. However, the effects of aging on their dynamics, heterogeneity, and transcriptional profiles remain to be fully elucidated. Through single cell RNA sequencing (scRNA-seq), we identified CBFß as an indispensable transcription factor that ensures AM self-renewal. Intriguingly, despite transcriptome similarities of proliferating cells, AMs from aged mice exhibited reduced embryonic stem cell-like features. Aged AMs also displayed compromised DNA repair abilities, potentially leading to obstructed cell cycle progression and an elevation of senescence markers. Consistently, AMs from aged mice exhibited impaired self-renewal ability and reduced sensitivity to GM-CSF. Decreased CBFß was observed in the cytosol of AMs from aged mice. Similar senescence-like phenotypes were also found in human AMs. Taken together, these findings suggest that AMs in aged hosts demonstrate senescence-like phenotypes, potentially facilitated by the abrogated CBF ß activity.

Proximal immune-epithelial progenitor interactions drive chronic tissue sequelae post COVID-19.

The long-term physiological consequences of SARS-CoV-2, termed Post-Acute Sequelae of COVID-19 (PASC), are rapidly evolving into a major public health concern. The underlying cellular and molecular etiology remain poorly defined but growing evidence links PASC to abnormal immune responses and/or poor organ recovery post-infection. Yet, the precise mechanisms driving non-resolving inflammation and impaired tissue repair in the context of PASC remain unclear. With insights from three independent clinical cohorts of PASC patients with abnormal lung function and/or viral infection-mediated pulmonary fibrosis, we established a clinically relevant mouse model of post-viral lung sequelae to investigate the pathophysiology of respiratory PASC. By employing a combination of spatial transcriptomics and imaging, we identified dysregulated proximal interactions between immune cells and epithelial progenitors unique to the fibroproliferation in respiratory PASC but not acute COVID-19 or idiopathic pulmonary fibrosis (IPF). Specifically, we found a central role for lung-resident CD8+ T cell-macrophage interactions in maintaining Krt8hi transitional and ectopic Krt5+ basal cell progenitors, thus impairing alveolar regeneration and driving fibrotic sequelae after acute viral pneumonia. Mechanistically, CD8+ T cell derived IFN-γ and TNF stimulated lung macrophages to chronically release IL-1ß, resulting in the abnormal accumulation of dysplastic epithelial progenitors and fibrosis. Notably, therapeutic neutralization of IFN-γ and TNF, or IL-1ß after the resolution of acute infection resulted in markedly improved alveolar regeneration and restoration of pulmonary function. Together, our findings implicate a dysregulated immune-epithelial progenitor niche in driving respiratory PASC. Moreover, in contrast to other approaches requiring early intervention, we highlight therapeutic strategies to rescue fibrotic disease in the aftermath of respiratory viral infections, addressing the current unmet need in the clinical management of PASC and post-viral disease.

Proximal immune-epithelial progenitor interactions drive chronic tissue sequelae post COVID-19.

The long-term physiological consequences of SARS-CoV-2, termed Post-Acute Sequelae of COVID-19 (PASC), are rapidly evolving into a major public health concern. The underlying cellular and molecular etiology remain poorly defined but growing evidence links PASC to abnormal immune responses and/or poor organ recovery post-infection. Yet, the precise mechanisms driving non-resolving inflammation and impaired tissue repair in the context of PASC remain unclear. With insights from three independent clinical cohorts of PASC patients with abnormal lung function and/or viral infection-mediated pulmonary fibrosis, we established a clinically relevant mouse model of post-viral lung sequelae to investigate the pathophysiology of respiratory PASC. By employing a combination of spatial transcriptomics and imaging, we identified dysregulated proximal interactions between immune cells and epithelial progenitors unique to the fibroproliferation in respiratory PASC but not acute COVID-19 or idiopathic pulmonary fibrosis (IPF). Specifically, we found a central role for lung-resident CD8+ T cell-macrophage interactions in maintaining Krt8hi transitional and ectopic Krt5+ basal cell progenitors, thus impairing alveolar regeneration and driving fibrotic sequelae after acute viral pneumonia. Mechanistically, CD8+ T cell derived IFN-γ and TNF stimulated lung macrophages to chronically release IL-1ß, resulting in the abnormal accumulation of dysplastic epithelial progenitors and fibrosis. Notably, therapeutic neutralization of IFN-γ and TNF, or IL-1ß after the resolution of acute infection resulted in markedly improved alveolar regeneration and restoration of pulmonary function. Together, our findings implicate a dysregulated immune-epithelial progenitor niche in driving respiratory PASC. Moreover, in contrast to other approaches requiring early intervention, we highlight therapeutic strategies to rescue fibrotic disease in the aftermath of respiratory viral infections, addressing the current unmet need in the clinical management of PASC and post-viral disease.

Intrinsic bias estimation for improved analysis of bulk and single-cell chromatin accessibility profiles using SELMA.

Genome-wide profiling of chromatin accessibility by DNase-seq or ATAC-seq has been widely used to identify regulatory DNA elements and transcription factor binding sites. However, enzymatic DNA cleavage exhibits intrinsic sequence biases that confound chromatin accessibility profiling data analysis. Existing computational tools are limited in their ability to account for such intrinsic biases and not designed for analyzing single-cell data. Here, we present Simplex Encoded Linear Model for Accessible Chromatin (SELMA), a computational method for systematic estimation of intrinsic cleavage biases from genomic chromatin accessibility profiling data. We demonstrate that SELMA yields accurate and robust bias estimation from both bulk and single-cell DNase-seq and ATAC-seq data. SELMA can utilize internal mitochondrial DNA data to improve bias estimation. We show that transcription factor binding inference from DNase footprints can be improved by incorporating estimated biases using SELMA. Furthermore, we show strong effects of intrinsic biases in single-cell ATAC-seq data, and develop the first single-cell ATAC-seq intrinsic bias correction model to improve cell clustering. SELMA can enhance the performance of existing bioinformatics tools and improve the analysis of both bulk and single-cell chromatin accessibility sequencing data.

Single-cell chromatin profiling of the primitive gut tube reveals regulatory dynamics underlying lineage fate decisions.

Development of the gastrointestinal system occurs after gut tube closure, guided by spatial and temporal control of gene expression. However, it remains unclear what forces regulate these spatiotemporal gene expression patterns. Here we perform single-cell chromatin profiling of the primitive gut tube to reveal organ-specific chromatin patterns that reflect the anatomical patterns of distinct organs. We generate a comprehensive map of epigenomic changes throughout gut development, demonstrating that dynamic chromatin accessibility patterns associate with lineage-specific transcription factor binding events to regulate organ-specific gene expression. Additionally, we show that loss of Sox2 and Cdx2, foregut and hindgut lineage-specific transcription factors, respectively, leads to fate shifts in epigenomic patterns, linking transcription factor binding, chromatin accessibility, and lineage fate decisions in gut development. Notably, abnormal expression of Sox2 in the pancreas and intestine impairs lineage fate decisions in both development and adult homeostasis. Together, our findings define the chromatin and transcriptional mechanisms of organ identity and lineage plasticity in development and adult homeostasis.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk.

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements. Here, we applied single-nucleus assay for transposase-accessible chromatin with sequencing to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell-type-specific elements and transcription factors, and prioritized functional CAD risk variants. We identified elements in smooth muscle cell transition states (for example, fibromyocytes) and functional variants predicted to alter smooth muscle cell- and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single-nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.

Ectopic Tcf1 expression instills a stem-like program in exhausted CD8+ T cells to enhance viral and tumor immunity.

Exhausted CD8+ T (Tex) cells are dysfunctional due to persistent antigen exposure in chronic viral infection and tumor contexts. A stem cell-like Tex (Tex-stem) subset can self-renew and differentiate into terminally exhausted (Tex-term) cells. Here, we show that ectopic Tcf1 expression potently promoted the generation of Tex-stem cells in both a chronic viral infection and preclinical tumor models. Tcf1 overexpression diminished coinhibitory receptor expression and enhanced polycytokine-producing capacity while retaining a heightened responses to checkpoint blockade, leading to enhanced viral and tumor control. Mechanistically, ectopically expressed Tcf1 exploited existing and novel chromatin accessible sites as transcriptional enhancers or repressors and modulated the transcriptome by enforcing pre-existing expression patterns in Tex-stem cells, such as enhanced suppression of Blimp1 and Bim and acquisition of new downstream genes, including Mx1, Tox2, and Runx3. These findings reveal a pronounced impact of ectopic Tcf1 expression on Tex functional restoration and highlight the therapeutic potential of harnessing Tcf1-enforced transcriptional programs.

ncHMR detector: a computational framework to systematically reveal non-classical functions of histone modification regulators.

Recently, several non-classical functions of histone modification regulators (HMRs), independent of their known histone modification substrates and products, have been reported to be essential for specific cellular processes. However, there is no framework designed for identifying such functions systematically. Here, we develop ncHMR detector, the first computational framework to predict non-classical functions and cofactors of a given HMR, based on ChIP-seq data mining. We apply ncHMR detector in ChIP-seq data-rich cell types and predict non-classical functions of HMRs. Finally, we experimentally reveal that the predicted non-classical function of CBX7 is biologically significant for the maintenance of pluripotency.