Search

Here we present advancements in single-cell combinatorial indexed Assay for Transposase Accessible Chromatin (sciATAC) to measure chromatin accessibility that leverage nanowell chips to achieve atlas-scale cell throughput (>105cells) at low cost. The platform leverages the core of the sciATAC workflow where multiple indexed tagmentation reactions are performed, followed by pooling and distribution to a second set of reaction wells for polymerase chain reaction (PCR)-based indexing. In this work, we instead leverage a chip containing 5184 nanowells at the PCR stage of indexing, enabling a 52-fold improvement in scale and reduction in per-cell preparation costs. We detail three variants that balance cell throughput and depth of coverage, and apply these methods to banked mouse brain tissue, producing maps of cell types as well as neuronal subtypes that include integration with existing single-cell Assay for Transposase Accessible Chromatin (scATAC) and scRNA-seq data sets. Our optimized workflow achieves a high fraction of reads that fall within called peaks (>80%) and low cell doublet rates. The high cell coverage technique produces high unique reads per cell, while retaining high enrichment for open chromatin regions, enabling the assessment of >70,000 unique accessible loci on average for each cell profiled. When compared to current methods in the field, our technique provides similar or superior per-cell information with very low levels of cell-to-cell cross talk, and achieves this at a cost point much lower than existing assays.

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of theHomo sapiensandMus musculusgenomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and availableviathe ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL;https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environmentsviaCromwell. Access to the pipelines and dataviathe cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.Database URL:https://www.encodeproject.org/

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. Here we perform dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit the stem and progenitor stage. Collectively, these data provide insights into the temporally regulated synergy of the cis- and the trans-regulatory components underlying hematopoietic lineage commitment and differentiation., Lineage differentiation and commitment is driven by transcription regulators and chromatin changes. Here the authors report daily profiling of chromatin accessibility and transcriptome changes during human erythropoiesis, relating these changes to lineage potential between erythropoiesis and megakaryopoieis.

Targeted sequencing remains a valuable technique for clinical and research applications. However, many existing technologies suffer from pervasive guanine-cytosine (GC) sequence content bias, high input DNA requirements, and high cost for custom panels. We have developed Cas12a-Capture, a low-cost and highly scalable method for targeted sequencing. The method utilizes preprogrammed guide RNAs to direct CRISPR-Cas12a cleavage of double-stranded DNA in vitro and then takes advantage of the resulting four to five nucleotide overhangs for selective ligation with a custom sequencing adapter. Addition of a second sequencing adapter and enrichment for ligation products generates a targeted sequence library. We first performed a pilot experiment with 7176 guides targeting 3.5 Mb of DNA. Using these data, we modeled the sequence determinants of Cas12a-Capture efficiency, then designed an optimized set of 11,438 guides targeting 3.0 Mb. The optimized guide set achieves an average 64-fold enrichment of targeted regions with minimal GC bias. Cas12a-Capture variant calls had strong concordance with Illumina Platinum Genome calls, especially for single nucleotide variants, which could be improved by applying basic variant quality heuristics. We believe Cas12a-Capture has a wide variety of potential clinical and research applications and is amendable for selective enrichment for any double-stranded DNA template or genome.

Seasonal influenza results in 3 to 5 million cases of severe disease and 250,000 to 500,000 deaths annually. Macrophages have been implicated in both the resolution and progression of the disease, but the drivers of these outcomes are poorly understood. We probed mouse lung transcriptomic datasets using the Digital Cell Quantifier algorithm to predict immune cell subsets that correlated with mild or severe influenza A virus (IAV) infection outcomes. We identified a novel lung macrophage population that transcriptionally resembled small serosal cavity macrophages and correlated with mild disease. Until now, the study of serosal macrophage translocation in the context of infections has been neglected. Here, we show that pleural macrophages (PMs) migrate from the pleural cavity to the lung after infection with pH1N1 A/California/04/2009 IAV. We found that the depletion of PMs increased morbidity and pulmonary inflammation. There were increased proinflammatory cytokines in the pleural cavity and an influx of neutrophils within the lung. Our results show PMs are recruited to the lung during IAV infection and contribute to recovery from influenza. This study expands our knowledge of PM plasticity and provides a new source of lung macrophages independent of monocyte recruitment and local proliferation.GRAPHICAL ABSTRACT

High-throughput single-cell epigenomic assays can resolve cell type heterogeneity in complex tissues, however, spatial orientation is lost. Here, we present single-cell combinatorial indexing on Microbiopsies Assigned to Positions for the Assay for Transposase Accessible Chromatin, or sciMAP-ATAC, as a method for highly scalable, spatially resolved, single-cell profiling of chromatin states. sciMAP-ATAC produces data of equivalent quality to non-spatial sci-ATAC and retains the positional information of each cell within a 214 micron cubic region, with up to hundreds of tracked positions in a single experiment. We apply sciMAP-ATAC to assess cortical lamination in the adult mouse primary somatosensory cortex and in the human primary visual cortex, where we produce spatial trajectories and integrate our data with non-spatial single-nucleus RNA and other chromatin accessibility single-cell datasets. Finally, we characterize the spatially progressive nature of cerebral ischemic infarction in the mouse brain using a model of transient middle cerebral artery occlusion., Spatial orientation of cells in an interconnected network is lost in high-throughput single-cell epigenomic assays. Here the authors present sciMAP-ATAC to produce spatially resolved single-cell ATAC-seq data.

Despite the known heritable nature of autism spectrum disorder (ASD), studies have primarily identified risk genes with de novo variants (DNVs). To capture the full spectrum of ASD genetic risk, we performed a two-stage analysis of rare de novo and inherited coding variants in 42,607 ASD cases, including 35,130 new cases recruited online by SPARK. In the first stage, we analyzed 19,843 cases with one or both biological parents and found that known ASD or neurodevelopmental disorder (NDD) risk genes explain nearly 70% of the genetic burden conferred by DNVs. In contrast, less than 20% of genetic risk conferred by rare inherited loss-of-function (LoF) variants are explained by known ASD/NDD genes. We selected 404 genes based on the first stage of analysis and performed a meta-analysis with an additional 22,764 cases and 236,000 population controls. We identified 60 genes with exome-wide significance (p < 2.5e-6), including five new risk genes (NAV3, ITSN1, MARK2, SCAF1, and HNRNPUL2). The association of NAV3 with ASD risk is entirely driven by rare inherited LoFs variants, with an average relative risk of 4, consistent with moderate effect. ASD individuals with LoF variants in the four moderate risk genes (NAV3, ITSN1, SCAF1, and HNRNPUL2, n = 95) have less cognitive impairment compared to 129 ASD individuals with LoF variants in well-established, highly penetrant ASD risk genes (CHD8, SCN2A, ADNP, FOXP1, SHANK3) (59% vs. 88%, p= 1.9e-06). These findings will guide future gene discovery efforts and suggest that much larger numbers of ASD cases and controls are needed to identify additional genes that confer moderate risk of ASD through rare, inherited variants.

To capture the full spectrum of genetic risk for autism, we performed a two-stage analysis of rare de novo and inherited coding variants in 42,607 autism cases, including 35,130 new cases recruited online by SPARK. We identified 60 genes with exome-wide significance (P 2.5 × 10

Targeted sequencing remains a valuable technique for clinical and research applications. However, many existing technologies suffer from pervasive GC sequence content bias, high input DNA requirements, and high cost for custom panels. We have developed Cas12a-Capture, a low-cost and highly scalable method for targeted sequencing. The method utilizes preprogramed guide RNAs to direct CRISPR-Cas12a cleavage of double stranded DNAin vitroand then takes advantage of the resulting four to five nucleotide overhangs for selective ligation with a custom sequencing adapter. Addition of a second sequencing adapter and enrichment for ligation products generates a targeted sequence library. We first performed a pilot experiment with 7,176 guides targeting 3.5 megabases of DNA. Using these data, we modeled the sequence determinants of Cas12a-Capture efficiency, then designed an optimized set of 11,438 guides targeting 3.0 megabases. The optimized guide set achieves an average 64-fold enrichment of targeted regions with minimal GC bias. Cas12a-Capture variant calls had strong concordance with Illumina Platinum Genome calls, especially for SNVs, which could be improved by applying basic variant quality heuristics. We believe Cas12a-Capture has a wide variety of potential clinical and research applications and is amendable for selective enrichment for any double stranded DNA template or genome.

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. We performed dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit pluripotency. Collectively, these data provide novel insights into the global regulatory landscape during hematopoiesis.

Triple negative breast cancers (TNBC) constitute one-sixth of invasive female breast cancer cases and are the most likely to develop resistance to treatment via genetic and/or epigenetic adaptation into drug tolerant persister (DTP) states. We applied single-cell ATAC-seq and RNA-seq to characterize the dynamic regulatory and transcriptional landscape in five basal-like TNBC cell lines in response to the MEK inhibitor Trametinib. We observed surprisingly few shared changes between lines, indicating substantial heterogeneity in the emergence of DTP states. However, we identified a shift toward a common state based on the novel observation of the preferential loss of cell line-specific regulatory elements and gene expression. Integration of the two modalities enabled a granular dissection of dynamic regulatory mechanisms, which revealed highly context-dependent roles of regulatory elements. This work highlights the heterogeneity of response, yet suggests homogenization occurs in the form of the preferential loss of epigenetic configurations unique to each BCCL.

Natural hosts of simian immunodeficiency virus (SIV) avoid AIDS despite lifelong infection. Here, we examined how this outcome is achieved by comparing a natural SIV host, African green monkey (AGM) to an AIDS susceptible species, rhesus macaque (RM). To asses gene expression profiles from acutely SIV infected AGMs and RMs, we developed a systems biology approach termed Conserved Gene Signature Analysis (CGSA), which compared RNA sequencing data from rectal AGM and RM tissues to various other species. We found that AGMs rapidly activate, and then maintain, evolutionarily conserved regenerative wound healing mechanisms in mucosal tissue. The wound healing protein fibronectin shows distinct tissue distribution and abundance kinetics in AGMs. Furthermore, AGM monocytes exhibit an embryonic development and repair/regeneration signature featuring TGF-β and concomitant reduced expression of inflammatory genes compared to RMs. This regenerative wound healing process likely preserves mucosal integrity and prevents inflammatory insults that underlie immune exhaustion in RMs., Here, the authors compare gene expression signatures in rectal tissues of African green monkeys (AGMs) and rhesus macaques (RMs) acutely infected with simian immunodeficiency virus and find that AGMs rapidly activate and maintain evolutionarily conserved regenerative wound healing mechanisms.

Ebola virus (EBOV) infections are characterized by a pronounced lymphopenia that is highly correlative with fatalities. However, the mechanisms leading to T-cell depletion remain largely unknown. Here, we demonstrate that both viral mRNAs and antigens are detectable in CD4+ T cells despite the absence of productive infection. A protein phosphatase 1 inhibitor, 1E7-03, and siRNA-mediated suppression of viral antigens were used to demonstrate de novo synthesis of viral RNAs and antigens in CD4+ T cells, respectively. Cell-to-cell fusion of permissive Huh7 cells with non-permissive Jurkat T cells impaired productive EBOV infection suggesting the presence of a cellular restriction factor. We determined that viral transcription is partially impaired in the fusion T cells. Lastly, we demonstrate that exposure of T cells to EBOV resulted in autophagy through activation of ER-stress related pathways. These data indicate that exposure of T cells to EBOV results in an abortive infection, which likely contributes to the lymphopenia observed during EBOV infections., Author summary Lymphopenia is a common characteristic of the disease caused by EBOV. We determined that despite the apparent lack of productive infection, EBOV is capable of entering T cells and producing both viral RNAs and proteins. Furthermore, we demonstrate that EBOV causes an abortive infection in T cells due to the presence of a cellular restriction factor. The abortive infection was associated with cell death following ER-stress induced autophagy. Collectively, these findings suggest that abortive infection in T cells is likely to contribute to lymphopenia during Ebola virus disease, which is uniformly linked with the severity of the disease. All EBOV vaccine candidates utilize GP as the sole antigen inducing a protective antibody response and in some clinical trials were shown to induce adverse side effects. The present study suggests that these effects can be associated with GP, which may lead to abortive infection of the vaccine construct in T cells contributing to the inflammatory response to the vaccines.

Most bacteria exist in aggregated, three-dimensional structures called biofilms. Although biofilms play important ecological roles in natural and engineered settings, they can also pose societal problems, for example, when they grow in plumbing systems or on medical implants. Understanding the processes that promote the growth and disassembly of biofilms could lead to better strategies to manage these structures. We had previously shown that Helicobacter pylori bacteria are repulsed by high concentrations of a self-produced molecule, AI-2, and that H. pylori mutants deficient in AI-2 sensing form larger and more homogeneously spaced biofilms. Here, we used computer simulations of biofilm formation to show that local H. pylori behavior of repulsion from high AI-2 could explain the overall architecture of H. pylori biofilms. Our findings demonstrate that it is possible to change global biofilm organization by manipulating local cell behaviors, which suggests that simple strategies targeting cells at local scales could be useful for controlling biofilms in industrial and medical settings., Bacteria are often found living in aggregated multicellular communities known as biofilms. Biofilms are three-dimensional structures that confer distinct physical and biological properties to the collective of cells living within them. We used agent-based modeling to explore whether local cellular interactions were sufficient to give rise to global structural features of biofilms. Specifically, we asked whether chemorepulsion from a self-produced quorum-sensing molecule, autoinducer-2 (AI-2), was sufficient to recapitulate biofilm growth and cellular organization observed for biofilms of Helicobacter pylori, a common bacterial resident of human stomachs. To carry out this modeling, we modified an existing platform, Individual-based Dynamics of Microbial Communities Simulator (iDynoMiCS), to incorporate three-dimensional chemotaxis, planktonic cells that could join or leave the biofilm structure, and cellular production of AI-2. We simulated biofilm growth of previously characterized H. pylori strains with various AI-2 production and sensing capacities. Using biologically plausible parameters, we were able to recapitulate both the variation in biofilm mass and cellular distributions observed with these strains. Specifically, the strains that were competent to chemotax away from AI-2 produced smaller and more heterogeneously spaced biofilms, whereas the AI-2 chemotaxis-defective strains produced larger and more homogeneously spaced biofilms. The model also provided new insights into the cellular demographics contributing to the biofilm patterning of each strain. Our analysis supports the idea that cellular interactions at small spatial and temporal scales are sufficient to give rise to larger-scale emergent properties of biofilms. IMPORTANCE Most bacteria exist in aggregated, three-dimensional structures called biofilms. Although biofilms play important ecological roles in natural and engineered settings, they can also pose societal problems, for example, when they grow in plumbing systems or on medical implants. Understanding the processes that promote the growth and disassembly of biofilms could lead to better strategies to manage these structures. We had previously shown that Helicobacter pylori bacteria are repulsed by high concentrations of a self-produced molecule, AI-2, and that H. pylori mutants deficient in AI-2 sensing form larger and more homogeneously spaced biofilms. Here, we used computer simulations of biofilm formation to show that local H. pylori behavior of repulsion from high AI-2 could explain the overall architecture of H. pylori biofilms. Our findings demonstrate that it is possible to change global biofilm organization by manipulating local cell behaviors, which suggests that simple strategies targeting cells at local scales could be useful for controlling biofilms in industrial and medical settings.

There are no approved therapeutics for the treatment of dengue disease despite the global prevalence of dengue virus (DENV) and its mosquito vectors. DENV infections can lead to vascular complications, hemorrhage, and shock due to the ability of DENV to infect a variety of immune and nonimmune cell populations. Increasingly, studies have implicated the host response as a major contributor to severe disease. Inflammatory products of various cell types, including responding T cells, mast cells (MCs), and infected monocytes, can contribute to immune pathology. In this study, we show that the host response to DENV infection in immunocompetent mice recapitulates transcriptional changes that have been described in human studies. We found that DENV infection strongly induced metabolic dysregulation, complement signaling, and inflammation. DENV also affected the immune cell content of the spleen and liver, enhancing NK, NKT, and CD8 + T cell activation. The MC-stabilizing drug ketotifen reversed many of these responses without suppressing memory T cell formation and induced additional changes in the transcriptome and immune cell composition of the spleen, consistent with reduced inflammation. This study provides a global transcriptional map of immune activation in DENV target organs of an immunocompetent host and supports the further development of targeted immunomodulatory strategies to treat DENV disease. IMPORTANCE Dengue virus (DENV), which causes febrile illness, is transmitted by mosquito vectors throughout tropical and subtropical regions of the world. Symptoms of DENV infection involve damage to blood vessels and, in rare cases, hemorrhage and shock. Currently, there are no targeted therapies to treat DENV infection, but it is thought that drugs that target the host immune response may be effective in limiting symptoms that result from excessive inflammation. In this study, we measured the host transcriptional response to infection in multiple DENV target organs using a mouse model of disease. We found that DENV infection induced metabolic dysregulation and inflammatory responses and affected the immune cell content of the spleen and liver. The use of the mast cell stabilization drug ketotifen reversed many of these responses and induced additional changes in the transcriptome and immune cell repertoire that contribute to decreased dengue disease.

Ebola virus (EBOV) and Reston virus (RESTV) are members of the Ebolavirus genus which greatly differ in their pathogenicity. While EBOV causes a severe disease in humans characterized by a dysregulated inflammatory response and elevated cytokine and chemokine production, there are no reported disease-associated human cases of RESTV infection, suggesting that RESTV is nonpathogenic for humans. The underlying mechanisms determining the pathogenicity of different ebolavirus species are not yet known. In this study, we dissected the host response to EBOV and RESTV infection in primary human monocyte-derived macrophages (MDMs). As expected, EBOV infection led to a profound proinflammatory response, including strong induction of type I and type III interferons (IFNs). In contrast, RESTV-infected macrophages remained surprisingly silent. Early activation of IFN regulatory factor 3 (IRF3) and NF-κB was observed in EBOV-infected, but not in RESTV-infected, MDMs. In concordance with previous results, MDMs treated with inactivated EBOV and Ebola virus-like particles (VLPs) induced NF-κB activation mediated by Toll-like receptor 4 (TLR4) in a glycoprotein (GP)-dependent manner. This was not the case in cells exposed to live RESTV, inactivated RESTV, or VLPs containing RESTV GP, indicating that RESTV GP does not trigger TLR4 signaling. Our results suggest that the lack of immune activation in RESTV-infected MDMs contributes to lower pathogenicity by preventing the cytokine storm observed in EBOV infection. We further demonstrate that inhibition of TLR4 signaling abolishes EBOV GP-mediated NF-κB activation. This finding indicates that limiting the excessive TLR4-mediated proinflammatory response in EBOV infection should be considered as a potential supportive treatment option for EBOV disease. IMPORTANCE Emerging infectious diseases are a major public health concern, as exemplified by the recent devastating Ebola virus (EBOV) outbreak. Different ebolavirus species are associated with widely varying pathogenicity in humans, ranging from asymptomatic infections for Reston virus (RESTV) to severe disease with fatal outcomes for EBOV. In this comparative study of EBOV- and RESTV-infected human macrophages, we identified key differences in host cell responses. Consistent with previous data, EBOV infection is associated with a proinflammatory signature triggered by the surface glycoprotein (GP), which can be inhibited by blocking TLR4 signaling. In contrast, infection with RESTV failed to stimulate a strong host response in infected macrophages due to the inability of RESTV GP to stimulate TLR4. We propose that disparate proinflammatory host signatures contribute to the differences in pathogenicity reported for ebolavirus species and suggest that proinflammatory pathways represent an intriguing target for the development of novel therapeutics.

Fatal outcomes of Ebola virus (EBOV) infections are typically preceded by a ‘sepsis-like’ syndrome and lymphopenia despite T cells being resistant to Ebola infection. The mechanisms that lead to T lymphocytes death remain largely unknown; however, the degree of lymphopenia is highly correlative with fatalities. Here we investigated whether the addition of EBOV or its envelope glycoprotein (GP) to isolated primary human CD4+ T cells induced cell death. We observed a significant decrease in cell viability in a GP-dependent manner, which is suggestive of a direct role of GP in T cell death. Using immunoprecipitation assays and flow cytometry, we demonstrate that EBOV directly binds to CD4+ T cells through interaction of GP with TLR4. Transcriptome analysis revealed that the addition of EBOV to CD4+ T cells results in the significant upregulation of pathways associated with interferon signaling, pattern recognition receptors and intracellular activation of NFκB signaling pathway. Both transcriptome analysis and specific inhibitors allowed identification of apoptosis and necrosis as mechanisms associated with the observed T cell death following exposure to EBOV. The addition of the TLR4 inhibitor CLI-095 significantly reduced CD4+ T cell death induced by GP. EBOV stimulation of primary CD4+ T cells resulted in a significant increase in secreted TNFα; inhibition of TNFα-mediated signaling events significantly reduced T cell death while inhibitors of both necrosis and apoptosis similarly reduced EBOV-induced T cell death. Lastly, we show that stimulation with EBOV or GP augments monocyte maturation as determined by an overall increase in expression levels of markers of differentiation. Subsequently, the increased rates of cellular differentiation resulted in higher rates of infection further contributing to T cell death. These results demonstrate that GP directly subverts the host’s immune response by increasing the susceptibility of monocytes to EBOV infection and triggering lymphopenia through direct and indirect mechanisms., Author summary The latest outbreak of Ebola virus (EBOV) in West Africa resulted in more than 28,000 human infections including more than 11,000 deaths thus highlighting the necessity for the development of countermeasures. Monocytes and dendritic cells are among the primary targets of EBOV infection; infection of these critical antigen presenting cells contributes to the immune deficiency observed in Ebola virus disease (EVD). In contrast, lymphocytes are resistant to EBOV infection; however, in fatal EVD, pronounced lymphopenia is uniformly observed. Here we report that T lymphocyte cell death in the absence of detectable infection was observed in an EBOV glycoprotein (GP)-dependent manner. Using transcriptome analysis of EBOV-stimulated CD4+ T cells we show upregulation of both toll-like receptor 4 (TLR4) and cell death associated pathways. Furthermore, we demonstrate that EBOV increases susceptibility of monocytes to infection by promoting cellular differentiation. Both EBOV-induced monocyte differentiation and cell death of T lymphocytes result from a direct interaction between GP and TLR4. Blocking of TLR4 signaling significantly reduced both EBOV-induced T cell death and infection of monocytes. These data contribute to understanding of the ‘immune paralysis’ during EBOV infections and provide evidence for the development of targeted therapies for the treatment of EVD.

Natural hosts of simian immunodeficiency virus (SIV) avoid AIDS despite lifelong infection. Here, we examined how this outcome is achieved by comparing a natural SIV host, African green monkey (AGM) to an AIDS susceptible species, rhesus macaque (RM). To asses gene expression profiles from acutely SIV infected AGMs and RMs, we developed a systems biology approach termed Conserved Gene Signature Analysis (CGSA), which compared RNA sequencing data from rectal AGM and RM tissues to various other species. We found that AGMs rapidly activate, and then maintain, evolutionarily conserved regenerative wound healing mechanisms in mucosal tissue. The wound healing protein fibronectin shows distinct tissue distribution and abundance kinetics in AGMs. Furthermore, AGM monocytes exhibit an embryonic development and repair/regeneration signature featuring TGF-β and concomitant reduced expression of inflammatory genes compared to RMs. This regenerative wound healing process likely preserves mucosal integrity and prevents inflammatory insults that underlie immune exhaustion in RMs.

The non-human primate reference transcriptome resource (NHPRTR, available online at http://nhprtr.org/) aims to generate comprehensive RNA-seq data from a wide variety of non-human primates (NHPs), from lemurs to hominids. In the 2012 Phase I of the NHPRTR project, 19 billion fragments or 3.8 terabases of transcriptome sequences were collected from pools of ∼ 20 tissues in 15 species and subspecies. Here we describe a major expansion of NHPRTR by adding 10.1 billion fragments of tissue-specific RNA-seq data. For this effort, we selected 11 of the original 15 NHP species and subspecies and constructed total RNA libraries for the same ∼ 15 tissues in each. The sequence quality is such that 88% of the reads align to human reference sequences, allowing us to compute the full list of expression abundance across all tissues for each species, using the reads mapped to human genes. This update also includes improved transcript annotations derived from RNA-seq data for rhesus and cynomolgus macaques, two of the most commonly used NHP models and additional RNA-seq data compiled from related projects. Together, these comprehensive reference transcriptomes from multiple primates serve as a valuable community resource for genome annotation, gene dynamics and comparative functional analysis.

The severity of influenza virus infection is largely determined by the interplay between the virus and the host. Understanding the pathological drivers of influenza A viruses serves to rapidly characterize the host response to newly emerging influenza viruses and derive prognostic signatures of disease severity. We implemented a differential network approach to uncover novel distinctions between the global transcriptional host response to highly pathogenic and low pathogenic H1N1, H7N9, H7N7, and H5N1 influenza A infections in BALB/c mice. We identified crucial mediators of inflammatory, coagulation, and tissue repair responses that discriminate between disease outcomes. It has been demonstrated that these responses are largely impacted by the immune cell influx into the site of infection. We examined the contribution of resident and lung-infiltrating immune cell quantity and activation state to disease phenotypes by implementing digital cell quantification (DCQ), which allowed us to model the in vivo dynamics of immune cells across viral strains based on time and strain-dependent infected lung transcriptional profiles. Through linear regression modeling, we revealed that distinct monocyte, granulocyte, DC, and T cell populations predict disease morbidity. Additionally, we demonstrated that both the kinetics of CD8+ T cell infiltration and the extent of their activation state, serve as a key predictor of immunopathological consequences of influenza infection. Overall, we provide a multidimensional analysis that delineates how innate immune responses relate to aberrant tissue repair responses and drive subsequent host adaptive immune response that control the severity of influenza virus disease in the murine model.