A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility.

We applied a combinatorial indexing assay, sci-ATAC-seq, to profile genome-wide chromatin accessibility in ∼100,000 single cells from 13 adult mouse tissues. We identify 85 distinct patterns of chromatin accessibility, most of which can be assigned to cell types, and ∼400,000 differentially accessible elements. We use these data to link regulatory elements to their target genes, to define the transcription factor grammar specifying each cell type, and to discover in vivo correlates of heterogeneity in accessibility within cell types. We develop a technique for mapping single cell gene expression data to single-cell chromatin accessibility data, facilitating the comparison of atlases. By intersecting mouse chromatin accessibility with human genome-wide association summary statistics, we identify cell-type-specific enrichments of the heritability signal for hundreds of complex traits. These data define the in vivo landscape of the regulatory genome for common mammalian cell types at single-cell resolution.

Stable integrant-specific differences in bimodal HIV-1 expression patterns revealed by high-throughput analysis.

HIV-1 gene expression is regulated by host and viral factors that interact with viral motifs and is influenced by proviral integration sites. Here, expression variation among integrants was followed for hundreds of individual proviral clones within polyclonal populations throughout successive rounds of virus and cultured cell replication, with limited findings using CD4+ cells from donor blood consistent with observations in immortalized cells. Tracking clonal behavior by proviral "zip codes" indicated that mutational inactivation during reverse transcription was rare, while clonal expansion and proviral expression states varied widely. By sorting for provirus expression using a GFP reporter in the nef open reading frame, distinct clone-specific variation in on/off proportions were observed that spanned three orders of magnitude. Tracking GFP phenotypes over time revealed that as cells divided, their progeny alternated between HIV transcriptional activity and non-activity. Despite these phenotypic oscillations, the overall GFP+ population within each clone was remarkably stable, with clones maintaining clone-specific equilibrium mixtures of GFP+ and GFP- cells. Integration sites were analyzed for correlations between genomic features and the epigenetic phenomena described here. Integrants inserted in the sense orientation of genes were more frequently found to be GFP negative than those in the antisense orientation, and clones with high GFP+ proportions were more distal to repressive H3K9me3 peaks than low GFP+ clones. Clones with low frequencies of GFP positivity appeared to expand more rapidly than clones for which most cells were GFP+, even though the tested proviruses were Vpr-. Thus, much of the increase in the GFP- population in these polyclonal pools over time reflected differential clonal expansion. Together, these results underscore the temporal and quantitative variability in HIV-1 gene expression among proviral clones that are conferred in the absence of metabolic or cell-type dependent variability, and shed light on cell-intrinsic layers of regulation that affect HIV-1 population dynamics.

Predicting gene expression in the human malaria parasite Plasmodium falciparum using histone modification, nucleosome positioning, and 3D localization features.

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors.

Pulmonary osteoclast-like cells in silica induced pulmonary fibrosis.

The pathophysiology of silicosis is poorly understood, limiting development of therapies for those who have been exposed to the respirable particle. We explored mechanisms of silica-induced pulmonary fibrosis in human lung samples collected from patients with occupational exposure to silica and in a longitudinal mouse model of silicosis using multiple modalities including whole-lung single-cell RNA sequencing and histological, biochemical, and physiologic assessments. In addition to pulmonary inflammation and fibrosis, intratracheal silica challenge induced osteoclast-like differentiation of alveolar macrophages and recruited monocytes, driven by induction of the osteoclastogenic cytokine, receptor activator of nuclear factor κΒ ligand (RANKL) in pulmonary lymphocytes, and alveolar type II cells. Anti-RANKL monoclonal antibody treatment suppressed silica-induced osteoclast-like differentiation in the lung and attenuated pulmonary fibrosis. We conclude that silica induces differentiation of pulmonary osteoclast-like cells leading to progressive lung injury, likely due to sustained elaboration of bone-resorbing proteases and hydrochloric acid. Interrupting osteoclast-like differentiation may therefore constitute a promising avenue for moderating lung damage in silicosis.

Single-cell analysis of chromatin and expression reveals age- and sex-associated alterations in the human heart.

Sex differences and age-related changes in the human heart at the tissue, cell, and molecular level have been well-documented and many may be relevant for cardiovascular disease. However, how molecular programs within individual cell types vary across individuals by age and sex remains poorly characterized. To better understand this variation, we performed single-nucleus combinatorial indexing (sci) ATAC- and RNA-Seq in human heart samples from nine donors. We identify hundreds of differentially expressed genes by age and sex and find epigenetic signatures of variation in ATAC-Seq data in this discovery cohort. We then scale up our single-cell RNA-Seq analysis by combining our data with five recently published single nucleus RNA-Seq datasets of healthy adult hearts. We find variation such as metabolic alterations by sex and immune changes by age in differential expression tests, as well as alterations in abundance of cardiomyocytes by sex and neurons with age. In addition, we compare our adult-derived ATAC-Seq profiles to analogous fetal cell types to identify putative developmental-stage-specific regulatory factors. Finally, we train predictive models of cell-type-specific RNA expression levels utilizing ATAC-Seq profiles to link distal regulatory sequences to promoters, quantifying the predictive value of a simple TF-to-expression regulatory grammar and identifying cell-type-specific TFs. Our analysis represents the largest single-cell analysis of cardiac variation by age and sex to date and provides a resource for further study of healthy cardiac variation and transcriptional regulation at single-cell resolution.

Pulmonary osteoclast-like cells in silica induced pulmonary fibrosis.

The pathophysiology of silicosis is poorly understood, limiting development of therapies for those who have been exposed to the respirable particle. We explored the mechanisms of silica-induced pulmonary fibrosis in a mouse model using multiple modalities including whole-lung single-nucleus RNA sequencing. These analyses revealed that in addition to pulmonary inflammation and fibrosis, intratracheal silica challenge induced osteoclast-like differentiation of alveolar macrophages and recruited monocytes, driven by induction of the osteoclastogenic cytokine, receptor activator of nuclear factor-κB ligand (RANKL) in pulmonary lymphocytes and alveolar type II cells. Furthermore, anti-RANKL monoclonal antibody treatment suppressed silica-induced osteoclast-like differentiation in the lung and attenuated silica-induced pulmonary fibrosis. We conclude that silica induces osteoclast-like differentiation of distinct recruited and tissue resident monocyte populations, leading to progressive lung injury, likely due to sustained elaboration of bone resorbing proteases and hydrochloric acid. Interrupting osteoclast-like differentiation may therefore constitute a promising avenue for moderating lung damage in silicosis.

Airway epithelial interferon response to SARS-CoV-2 is inferior to rhinovirus and heterologous rhinovirus infection suppresses SARS-CoV-2 replication.

Common alphacoronaviruses and human rhinoviruses (HRV) induce type I and III interferon (IFN) responses important to limiting viral replication in the airway epithelium. In contrast, highly pathogenic betacoronaviruses including SARS-CoV-2 may evade or antagonize RNA-induced IFN I/III responses. In airway epithelial cells (AECs) from children and older adults we compared IFN I/III responses to SARS-CoV-2 and HRV-16, and assessed whether pre-infection with HRV-16, or pretreatment with recombinant IFN-ß or IFN-λ, modified SARS-CoV-2 replication. Bronchial AECs from children (ages 6-18 years) and older adults (ages 60-75 years) were differentiated ex vivo to generate organotypic cultures. In a biosafety level 3 (BSL-3) facility, cultures were infected with SARS-CoV-2 or HRV-16, and RNA and protein was harvested from cell lysates 96 h. following infection and supernatant was collected 48 and 96 h. following infection. In additional experiments cultures were pre-infected with HRV-16, or pre-treated with recombinant IFN-ß1 or IFN-λ2 before SARS-CoV-2 infection. In a subset of experiments a range of infectious concentrations of HRV-16, SARS-CoV-2 WA-01, SARS-CoV-2 Delta variant, and SARS-CoV-2 Omicron variant were studied. Despite significant between-donor heterogeneity SARS-CoV-2 replicated 100 times more efficiently than HRV-16. IFNB1, INFL2, and CXCL10 gene expression and protein production following HRV-16 infection was significantly greater than following SARS-CoV-2. IFN gene expression and protein production were inversely correlated with SARS-CoV-2 replication. Treatment of cultures with recombinant IFNß1 or IFNλ2, or pre-infection of cultures with HRV-16, markedly reduced SARS-CoV-2 replication. In addition to marked between-donor heterogeneity in IFN responses and viral replication, SARS-CoV-2 (WA-01, Delta, and Omicron variants) elicits a less robust IFN response in primary AEC cultures than does rhinovirus, and heterologous rhinovirus infection, or treatment with recombinant IFN-ß1 or IFN-λ2, reduces SARS-CoV-2 replication, although to a lesser degree for the Delta and Omicron variants.

Airway epithelial interferon response to SARS-CoV-2 is inferior to rhinovirus and heterologous rhinovirus infection suppresses SARS-CoV-2 replication.

INTRODUCTION: Common alphacoronaviruses and human rhinoviruses (HRV) induce type I and III interferon (IFN) responses important to limiting viral replication in the airway epithelium. In contrast, highly pathogenic betacoronaviruses including SARS-CoV-2 may evade or antagonize RNA-induced IFN I/III responses. METHODS: In airway epithelial cells (AECs) from children and older adults we compared IFN I/III responses to SARS-CoV-2 and HRV-16, and assessed whether pre-infection with HRV-16, or pretreatment with recombinant IFN-ß or IFN-λ, modified SARS-CoV-2 replication. Bronchial AECs from children (ages 6-18 yrs.) and older adults (ages 60-75 yrs.) were differentiated ex vivo to generate organotypic cultures. In a biosafety level 3 (BSL-3) facility, cultures were infected with SARS-CoV-2 or HRV-16, and RNA and protein was harvested from cell lysates 96 hrs. following infection and supernatant was collected 48 and 96 hrs. following infection. In additional experiments cultures were pre-infected with HRV-16, or pre-treated with recombinant IFN-ß1 or IFN-λ2 before SARS-CoV-2 infection. RESULTS: Despite significant between-donor heterogeneity SARS-CoV-2 replicated 100 times more efficiently than HRV-16. IFNB1, INFL2, and CXCL10 gene expression and protein production following HRV-16 infection was significantly greater than following SARS-CoV-2. IFN gene expression and protein production were inversely correlated with SARS-CoV-2 replication. Treatment of cultures with recombinant IFNß1 or IFNλ2, or pre-infection of cultures with HRV-16, markedly reduced SARS-CoV-2 replication. DISCUSSION: In addition to marked between-donor heterogeneity in IFN responses and viral replication, SARS-CoV-2 elicits a less robust IFN response in primary AEC cultures than does rhinovirus, and heterologous rhinovirus infection, or treatment with recombinant IFN-ß1 or IFN-λ2, markedly reduces SARS-CoV-2 replication.

An In Vitro Assay to Detect tRNA-Isopentenyl Transferase Activity.

N6-isopentenyladenosine RNA modifications are functionally diverse and highly conserved among prokaryotes and eukaryotes. One of the most highly conserved N6-isopentenyladenosine modifications occurs at the A37 position in a subset of tRNAs. This modification improves translation efficiency and fidelity by increasing the affinity of the tRNA for the ribosome. Mutation of enzymes responsible for this modification in eukaryotes are associated with several disease states, including mitochondrial dysfunction and cancer. Therefore, understanding the substrate specificity and biochemical activities of these enzymes is important for understanding of normal and pathologic eukaryotic biology. A diverse array of methods has been employed to characterize i6A modifications. Herein is described a direct approach for the detection of isopentenylation by Mod5. This method utilizes incubation of RNAs with a recombinant isopentenyl transferase, followed by RNase T1 digestion, and 1-dimensional gel electrophoresis analysis to detect i6A modifications. In addition, the potential adaptability of this protocol to characterize other RNA-modifying enzymes is discussed.

Aggregation of Mod5 is affected by tRNA binding with implications for tRNA gene-mediated silencing.

Mod5 is a multifunctional protein that modifies a subset of tRNAs in the cytoplasm and is also required for an RNA-mediated form of transcriptional silencing. Previous in vivo studies have shown that the nuclear silencing function of Mod5 does not require that the causative tRNA gene encode a Mod5 substrate, although Mod5 is still required. However, previous data have not directly tested whether Mod5 can directly bind substrate and nonsubstrate RNAs. We herein demonstrate that Mod5 directly binds to both substrate and nonsubstrate RNAs, including a highly structured, non-tRNA sequence (5S-rRNA), consistent with previous in vivo data. Furthermore, we show that some RNAs drastically change the aggregation behavior of Mod5 with implications for tRNA gene-mediated silencing.