Microbiota Diurnal Rhythmicity Programs Host Transcriptome Oscillations.

The intestinal microbiota undergoes diurnal compositional and functional oscillations that affect metabolic homeostasis, but the mechanisms by which the rhythmic microbiota influences host circadian activity remain elusive. Using integrated multi-omics and imaging approaches, we demonstrate that the gut microbiota features oscillating biogeographical localization and metabolome patterns that determine the rhythmic exposure of the intestinal epithelium to different bacterial species and their metabolites over the course of a day. This diurnal microbial behavior drives, in turn, the global programming of the host circadian transcriptional, epigenetic, and metabolite oscillations. Surprisingly, disruption of homeostatic microbiome rhythmicity not only abrogates normal chromatin and transcriptional oscillations of the host, but also incites genome-wide de novo oscillations in both intestine and liver, thereby impacting diurnal fluctuations of host physiology and disease susceptibility. As such, the rhythmic biogeography and metabolome of the intestinal microbiota regulates the temporal organization and functional outcome of host transcriptional and epigenetic programs.

The Spectrum and Regulatory Landscape of Intestinal Innate Lymphoid Cells Are Shaped by the Microbiome.

Innate lymphoid cells (ILCs) are critical modulators of mucosal immunity, inflammation, and tissue homeostasis, but their full spectrum of cellular states and regulatory landscapes remains elusive. Here, we combine genome-wide RNA-seq, ChIP-seq, and ATAC-seq to compare the transcriptional and epigenetic identity of small intestinal ILCs, identifying thousands of distinct gene profiles and regulatory elements. Single-cell RNA-seq and flow and mass cytometry analyses reveal compartmentalization of cytokine expression and metabolic activity within the three classical ILC subtypes and highlight transcriptional states beyond the current canonical classification. In addition, using antibiotic intervention and germ-free mice, we characterize the effect of the microbiome on the ILC regulatory landscape and determine the response of ILCs to microbial colonization at the single-cell level. Together, our work characterizes the spectrum of transcriptional identities of small intestinal ILCs and describes how ILCs differentially integrate signals from the microbial microenvironment to generate phenotypic and functional plasticity.

Transcriptional Heterogeneity and Lineage Commitment in Myeloid Progenitors.

Within the bone marrow, stem cells differentiate and give rise to diverse blood cell types and functions. Currently, hematopoietic progenitors are defined using surface markers combined with functional assays that are not directly linked with in vivo differentiation potential or gene regulatory mechanisms. Here, we comprehensively map myeloid progenitor subpopulations by transcriptional sorting of single cells from the bone marrow. We describe multiple progenitor subgroups, showing unexpected transcriptional priming toward seven differentiation fates but no progenitors with a mixed state. Transcriptional differentiation is correlated with combinations of known and previously undefined transcription factors, suggesting that the process is tightly regulated. Histone maps and knockout assays are consistent with early transcriptional priming, while traditional transplantation experiments suggest that in vivo priming may still allow for plasticity given strong perturbations. These data establish a reference model and general framework for studying hematopoiesis at single-cell resolution.

Tissue-resident macrophage enhancer landscapes are shaped by the local microenvironment.

Macrophages are critical for innate immune defense and also control organ homeostasis in a tissue-specific manner. They provide a fitting model to study the impact of ontogeny and microenvironment on chromatin state and whether chromatin modifications contribute to macrophage identity. Here, we profile the dynamics of four histone modifications across seven tissue-resident macrophage populations. We identify 12,743 macrophage-specific enhancers and establish that tissue-resident macrophages have distinct enhancer landscapes beyond what can be explained by developmental origin. Combining our enhancer catalog with gene expression profiles and open chromatin regions, we show that a combination of tissue- and lineage-specific transcription factors form the regulatory networks controlling chromatin specification in tissue-resident macrophages. The environment is capable of shaping the chromatin landscape of transplanted bone marrow precursors, and even differentiated macrophages can be reprogrammed when transferred into a new microenvironment. These results provide a comprehensive view of macrophage regulatory landscape and highlight the importance of the microenvironment, along with pioneer factors in orchestrating identity and plasticity.

The role of the local environment and epigenetics in shaping macrophage identity and their effect on tissue homeostasis.

Macrophages provide a critical systemic network cells of the innate immune system. Emerging data suggest that in addition, they have important tissue-specific functions that range from clearance of surfactant from the lungs to neuronal pruning and establishment of gut homeostasis. The differentiation and tissue-specific activation of macrophages require precise regulation of gene expression, a process governed by epigenetic mechanisms such as DNA methylation, histone modification and chromatin structure. We argue that epigenetic regulation of macrophages is determined by lineage- and tissue-specific transcription factors controlled by the built-in programming of myeloid development in combination with signaling from the tissue environment. Perturbation of epigenetic mechanisms of tissue macrophage identity can affect normal macrophage tissue function and contribute to pathologies ranging from obesity and autoimmunity to neurodegenerative diseases.

Distinct Murine Mucosal Langerhans Cell Subsets Develop from Pre-dendritic Cells and Monocytes.

Langerhans cells (LCs) populate the mucosal epithelium, a major entry portal for pathogens, yet their ontogeny remains unclear. We found that, in contrast to skin LCs originating from self-renewing radioresistant embryonic precursors, oral mucosal LCs derive from circulating radiosensitive precursors. Mucosal LCs can be segregated into CD103(+)CD11b(lo) (CD103(+)) and CD11b(+)CD103(-) (CD11b(+)) subsets. We further demonstrated that similar to non-lymphoid dendritic cells (DCs), CD103(+) LCs originate from pre-DCs, whereas CD11b(+) LCs differentiate from both pre-DCs and monocytic precursors. Despite this ontogenetic discrepancy between skin and mucosal LCs, the transcriptomic signature and immunological function of oral LCs highly resemble those of skin LCs but not DCs. These findings, along with the epithelial position, morphology, and expression of the LC-associated phenotype strongly suggest that oral mucosal LCs are genuine LCs. Collectively, in a tissue-dependent manner, murine LCs differentiate from at least three distinct precursors (embryonic, pre-DC, and monocytic) in steady state.

A Novel MIP-1-Expressing Macrophage Subtype in BAL Fluid from Healthy Volunteers.

Tissue availability remains an important limitation of single-cell genomic technologies for investigating cellular heterogeneity in human health and disease. BAL represents a minimally invasive approach to assessing an individual's lung cellular environment for diagnosis and research. However, the lack of high-quality, healthy lung reference data is a major obstacle to using single-cell approaches to study a plethora of lung diseases. Here, we performed single-cell RNA sequencing on over 40,000 cells isolated from the BAL of four healthy volunteers. Of the six cell types or lineages we identified, macrophages were consistently the most numerous across individuals. Our analysis confirmed the expression of marker genes defining cell types despite background signals because of the ambient RNA found in many single-cell studies. We assessed the variability of gene expression across macrophages and defined a distinct subpopulation of cells expressing a set of genes associated with Macrophage Inflammatory Protein 1 (MIP-1). RNA in situ hybridization and reanalysis of published lung single-cell data validated the presence of this macrophage subpopulation. Thus, our study characterizes lung macrophage heterogeneity in healthy individuals and provides a valuable resource for future studies to understand the lung environment in health and disease.

A negative feedback loop of transcription factors specifies alternative dendritic cell chromatin States.

During hematopoiesis, cells originating from the same stem cell reservoir differentiate into distinct cell types. The mechanisms enabling common progenitors to differentiate into alternative cell fates are not fully understood. Here, we identify cell-fate-determining transcription factors (TFs) governing dendritic cell (DC) development by annotating the enhancer landscapes of the DC lineage. Combining these analyses with detailed overexpression, knockdown, and ChIP-Seq studies, we show that Irf8 functions as a plasmacytoid DC epigenetic and fate-determining TF, regulating massive, cell-specific chromatin changes in thousands of pDC enhancers. Importantly, Irf8 forms a negative feedback loop with Cebpb, a monocyte-derived DC epigenetic fate-determining TF. We show that using this circuit logic, a pulse of TF expression can stably define epigenetic and transcriptional states, regardless of the microenvironment. More broadly, our study proposes a general paradigm that allows closely related cells with a similar set of signal-dependent factors to generate differential and persistent enhancer landscapes.

Activation of the 15-lipoxygenase pathway in aspirin-exacerbated respiratory disease.

BACKGROUND: Aspirin-exacerbated respiratory disease (AERD) is characterized by asthma, chronic rhinosinusitis with nasal polyps (CRSwNP), and an intolerance of medications that inhibit cyclooxygenase-1. Patients with AERD have more severe upper and lower respiratory tract disease than do aspirin-tolerant patients with CRSwNP. A dysregulation in arachidonic acid metabolism is thought to contribute to the enhanced sinonasal inflammation in AERD. OBJECTIVE: Our aim was to utilize an unbiased approach investigating arachidonic acid metabolic pathways in AERD. METHODS: Single-cell RNA sequencing (10× Genomics, Pleasanton, Calif) was utilized to compare the transcriptional profile of nasal polyp (NP) cells from patients with AERD and patients with CRSwNP and map differences in the expression of select genes among identified cell types. Findings were confirmed by traditional real-time PCR. Lipid mediators in sinonasal tissue were measured by mass spectrometry. Localization of various proteins within NPs was assessed by immunofluorescence. RESULTS: The gene encoding for 15-lipooxygenase (15-LO), ALOX15, was significantly elevated in NPs of patients with AERD compared to NPs of patients with CRSwNP (P < .05) or controls (P < .001). ALOX15 was predominantly expressed by epithelial cells. Expression levels significantly correlated with radiographic sinus disease severity (r = 0.56; P < .001) and were associated with asthma. The level of 15-oxo-eicosatetraenoic acid (15-Oxo-ETE), a downstream product of 15-LO, was significantly elevated in NPs from patients with CRSwNP (27.93 pg/mg of tissue) and NPs from patients with AERD (61.03 pg/mg of tissue) compared to inferior turbinate tissue from controls (7.17 pg/mg of tissue [P < .001]). Hydroxyprostaglandin dehydrogenase, an enzyme required for 15-Oxo-ETE synthesis, was predominantly expressed in mast cells and localized near 15-LO+ epithelium in NPs from patients with AERD. CONCLUSIONS: Epithelial and mast cell interactions, leading to the synthesis of 15-Oxo-ETE, may contribute to the dysregulation of arachidonic acid metabolism via the 15-LO pathway and to the enhanced sinonasal disease severity observed in AERD.

Microglia Adopt Longitudinal Transcriptional Changes After Traumatic Brain Injury.

BACKGROUND: Traumatic brain injury (TBI) is an under-recognized public health threat. Even mild brain injuries can lead to long-term neurologic impairment. Microglia play a fundamental role in the development and progression of this ensuing neurologic impairment. Despite this, a microglia-specific injury signature has yet to be identified. We hypothesized that TBI would lead to long-term changes in the transcriptional profile of microglial pathways associated with the development of subsequent neurologic impairment. MATERIALS AND METHODS: Male C57BL/6 mice underwent TBI via a controlled cortical impact and were followed longitudinally. FACSorted microglia from TBI mice were subjected to Quantiseq 3'-biased RNA sequencing at 7, 30, and 90 d after TBI. K-means clustering on 396 differentially expressed genes was performed, and gene ontology enrichment analysis was used to determine corresponding enriched processes. RESULTS: Differentially expressed genes in microglia exhibited four main patterns of expression over the course of TBI. In particular, we identified four gene clusters which corresponded to the host defense response, synaptic plasticity, lipid remodeling, and membrane polarization. CONCLUSIONS: Transcriptional profiling within individual populations of microglia after TBI remains a critical unmet research need within the field of TBI. This focused study identified several physiologic processes within microglia that may be associated with development of long-term neurologic impairment after TBI. These data demonstrate the capability of longitudinal transcriptional profiling to uncover potential cell-specific targets for the treatment of TBI.

Single-Cell Transcriptomic Analysis of Human Lung Provides Insights into the Pathobiology of Pulmonary Fibrosis.

Rationale: The contributions of diverse cell populations in the human lung to pulmonary fibrosis pathogenesis are poorly understood. Single-cell RNA sequencing can reveal changes within individual cell populations during pulmonary fibrosis that are important for disease pathogenesis. Objectives: To determine whether single-cell RNA sequencing can reveal disease-related heterogeneity within alveolar macrophages, epithelial cells, or other cell types in lung tissue from subjects with pulmonary fibrosis compared with control subjects. Methods: We performed single-cell RNA sequencing on lung tissue obtained from eight transplant donors and eight recipients with pulmonary fibrosis and on one bronchoscopic cryobiospy sample from a patient with idiopathic pulmonary fibrosis. We validated these data using in situ RNA hybridization, immunohistochemistry, and bulk RNA-sequencing on flow-sorted cells from 22 additional subjects. Measurements and Main Results: We identified a distinct, novel population of profibrotic alveolar macrophages exclusively in patients with fibrosis. Within epithelial cells, the expression of genes involved in Wnt secretion and response was restricted to nonoverlapping cells. We identified rare cell populations including airway stem cells and senescent cells emerging during pulmonary fibrosis. We developed a web-based tool to explore these data. Conclusions: We generated a single-cell atlas of pulmonary fibrosis. Using this atlas, we demonstrated heterogeneity within alveolar macrophages and epithelial cells from subjects with pulmonary fibrosis. These results support the feasibility of discovery-based approaches using next-generation sequencing technologies to identify signaling pathways for targeting in the development of personalized therapies for patients with pulmonary fibrosis.

Lipocalin-2 is a pathogenic determinant and biomarker of neuropsychiatric lupus.

Neuropsychiatric manifestations in lupus (NPSLE) affect ∼20-40% of patients. In the central nervous system, lipocalin-2 (LCN2) can promote injury through mechanisms directly linked to NPSLE, including brain barrier disruption, neurotoxicity, and glial activation. Since LCN2 is elevated in lupus and has been implicated in neuroinflammation, we investigated whether LCN2 is required for the pathogenesis of NPSLE. Here, we investigated the effects of LCN2 deficiency on the development of neurobehavioral deficits in the B6.Sle1.Sle3 (Sle1,3) mouse lupus model. Sle1,3 mice exhibited depression-like behavior and impaired spatial and recognition memory, and these deficits were attenuated in Sle1,3-LCN2KO mice. Whole-brain flow cytometry showed a significant increase in brain infiltrating leukocytes in Sle1,3 mice that was not reduced by LCN2 deficiency. RNA sequencing on sorted microglia revealed that several genes differentially expressed between B6 and Sle1,3 mice were regulated by LCN2, and that these genes are key mediators of the neuroinflammatory cascade. Importantly, LCN2 is upregulated in the cerebrospinal fluid of NPSLE patients across 2 different ethnicities. Our findings establish the Sle1,3 strain as an NPSLE model, demonstrate that LCN2 is a major regulator of the detrimental neuroimmune response in NPSLE, and identify CSF LCN2 as a novel biomarker for NPSLE.

A Beginner's Guide to Analysis of RNA Sequencing Data.

Since the first publications coining the term RNA-seq (RNA sequencing) appeared in 2008, the number of publications containing RNA-seq data has grown exponentially, hitting an all-time high of 2,808 publications in 2016 (PubMed). With this wealth of RNA-seq data being generated, it is a challenge to extract maximal meaning from these datasets, and without the appropriate skills and background, there is risk of misinterpretation of these data. However, a general understanding of the principles underlying each step of RNA-seq data analysis allows investigators without a background in programming and bioinformatics to critically analyze their own datasets as well as published data. Our goals in the present review are to break down the steps of a typical RNA-seq analysis and to highlight the pitfalls and checkpoints along the way that are vital for bench scientists and biomedical researchers performing experiments that use RNA-seq.

The role of chromatin dynamics in immune cell development.

In immune cells, as in all mammalian cells, nuclear DNA is wrapped around histones to form nucleosomes. The positioning and modifications of nucleosomes throughout the genome defines the chromatin state of the cell and has a large impact on gene regulation. Chromatin state is dynamic throughout immune cell development and activation. High-throughput open chromatin assays, such as DNase-seq, can be used to find regulatory element across the genome and, when combined with histone modifications, can specify their function. During hematopoiesis, distal regulatory elements, known as enhancers, are established by pioneer factors that alter chromatin state. Some of these enhancers are lost, some are gained, and some are maintained as a memory of the cell's developmental origin. The enhancer landscape is unique to the cell lineage-with different enhancers regulating the same promoter-and determines the mechanism of cell type-specific activation after exposure to stimuli. Histone modification and promoter architecture govern the diverse responses to stimulation. Furthermore, chromatin dynamics may explain the high plasticity of certain tissue-resident immune cell types. Future epigenomic research will depend on the development of more efficient experiments and better methods to associate enhancers with genes. The ultimate goal of mapping genome-wide chromatin state throughout the hematopoietic tree will help illuminate the mechanisms behind immune cell development and function.

DNase-seq predicts regions of rotational nucleosome stability across diverse human cell types.

DNase-seq is primarily used to identify nucleosome-depleted DNase I hypersensitive (DHS) sites genome-wide that correspond to active regulatory elements. However, ≈ 40 yr ago it was demonstrated that DNase I also digests with a ≈ 10-bp periodicity around nucleosomes matching the exposure of the DNA minor groove as it wraps around histones. Here, we use DNase-seq data from 49 samples representing diverse cell types to reveal this digestion pattern at individual loci and predict genomic locations where nucleosome rotational positioning, the orientation of DNA with respect to the histone surface, is stably maintained. We call these regions DNase I annotated regions of nucleosome stability (DARNS). Compared to MNase-seq experiments, we show DARNS correspond well to annotated nucleosomes. Interestingly, many DARNS are positioned over only one side of annotated nucleosomes, suggesting that the periodic digestion pattern attenuates over the nucleosome dyad. DARNS reproduce the arrangement of nucleosomes around transcription start sites and are depleted at ubiquitous DHS sites. We also generated DARNS from multiple lymphoblast cell line (LCL) samples. We found that LCL DARNS were enriched at DHS sites present in most of the original 49 samples but absent in LCLs, while multi-cell-type DARNS were enriched at LCL-specific DHS sites. This indicates that variably open DHS sites are often occupied by rotationally stable nucleosomes in cell types where the DHS site is closed. DARNS provide additional information about precise DNA orientation within individual nucleosomes not available from other nucleosome positioning assays and contribute to understanding the role of chromatin in gene regulation.

Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity.

The human body contains thousands of unique cell types, each with specialized functions. Cell identity is governed in large part by gene transcription programs, which are determined by regulatory elements encoded in DNA. To identify regulatory elements active in seven cell lines representative of diverse human cell types, we used DNase-seq and FAIRE-seq (Formaldehyde Assisted Isolation of Regulatory Elements) to map "open chromatin." Over 870,000 DNaseI or FAIRE sites, which correspond tightly to nucleosome-depleted regions, were identified across the seven cell lines, covering nearly 9% of the genome. The combination of DNaseI and FAIRE is more effective than either assay alone in identifying likely regulatory elements, as judged by coincidence with transcription factor binding locations determined in the same cells. Open chromatin common to all seven cell types tended to be at or near transcription start sites and to be coincident with CTCF binding sites, while open chromatin sites found in only one cell type were typically located away from transcription start sites and contained DNA motifs recognized by regulators of cell-type identity. We show that open chromatin regions bound by CTCF are potent insulators. We identified clusters of open regulatory elements (COREs) that were physically near each other and whose appearance was coordinated among one or more cell types. Gene expression and RNA Pol II binding data support the hypothesis that COREs control gene activity required for the maintenance of cell-type identity. This publicly available atlas of regulatory elements may prove valuable in identifying noncoding DNA sequence variants that are causally linked to human disease.

Transcription initiation patterns indicate divergent strategies for gene regulation at the chromatin level.

The application of deep sequencing to map 5' capped transcripts has confirmed the existence of at least two distinct promoter classes in metazoans: "focused" promoters with transcription start sites (TSSs) that occur in a narrowly defined genomic span and "dispersed" promoters with TSSs that are spread over a larger window. Previous studies have explored the presence of genomic features, such as CpG islands and sequence motifs, in these promoter classes, but virtually no studies have directly investigated the relationship with chromatin features. Here, we show that promoter classes are significantly differentiated by nucleosome organization and chromatin structure. Dispersed promoters display higher associations with well-positioned nucleosomes downstream of the TSS and a more clearly defined nucleosome free region upstream, while focused promoters have a less organized nucleosome structure, yet higher presence of RNA polymerase II. These differences extend to histone variants (H2A.Z) and marks (H3K4 methylation), as well as insulator binding (such as CTCF), independent of the expression levels of affected genes. Notably, differences are conserved across mammals and flies, and they provide for a clearer separation of promoter architectures than the presence and absence of CpG islands or the occurrence of stalled RNA polymerase. Computational models support the stronger contribution of chromatin features to the definition of dispersed promoters compared to focused start sites. Our results show that promoter classes defined from 5' capped transcripts not only reflect differences in the initiation process at the core promoter but also are indicative of divergent transcriptional programs established within gene-proximal nucleosome organization.