Aging of the blood system impacts systemic health and can be traced to hematopoietic stem cells (HSCs). Despite multiple reports on human HSC aging, a unified map detailing their molecular age-related changes is lacking. We developed a consensus map of gene expression in HSCs by integrating seven single-cell datasets. This map reveals previously unappreciated heterogeneity within the HSC population. It also links inflammatory pathway activation (TNF/NFκB, AP-1) and quiescence within a single gene expression program. This program dominates an inflammatory HSC subpopulation that increases with age, highlighting a potential target for further experimental studies and anti-aging interventions.

T cells recognize antigens and induce specialized gene expression programs (GEPs), enabling functions like proliferation, cytotoxicity and cytokine production. Traditionally, different T cell classes are thought to exhibit mutually exclusive responses, including TH1, TH2 and TH17 programs. However, single-cell RNA sequencing has revealed a continuum of T cell states without clearly distinct subsets, necessitating new analytical frameworks. Here, we introduce T-CellAnnoTator (TCAT), a pipeline that improves T cell characterization by simultaneously quantifying predefined GEPs capturing activation states and cellular subsets. Analyzing 1,700,000 T cells from 700 individuals spanning 38 tissues and five disease contexts, we identify 46 reproducible GEPs reflecting core T cell functions including proliferation, cytotoxicity, exhaustion and effector states. We experimentally demonstrate new activation programs and apply TCAT to characterize activation GEPs that predict immune checkpoint inhibitor response across multiple tumor types. Our software package starCAT generalizes this framework, enabling reproducible annotation in other cell types and tissues.

Infection with Lassa virus (LASV) can cause Lassa fever, a haemorrhagic illness with an estimated fatality rate of 29.7%, but causes no or mild symptoms in many individuals. Here, to investigate whether human genetic variation underlies the heterogeneity of LASV infection, we carried out genome-wide association studies (GWAS) as well as seroprevalence surveys, human leukocyte antigen typing and high-throughput variant functional characterization assays. We analysed Lassa fever susceptibility and fatal outcomes in 533 cases of Lassa fever and 1,986 population controls recruited over a 7 year period in Nigeria and Sierra Leone. We detected genome-wide significant variant associations with Lassa fever fatal outcomes near GRM7 and LIF in the Nigerian cohort. We also show that a haplotype bearing signatures of positive selection and overlapping LARGE1, a required LASV entry factor, is associated with decreased risk of Lassa fever in the Nigerian cohort but not in the Sierra Leone cohort. Overall, we identified variants and genes that may impact the risk of severe Lassa fever, demonstrating how GWAS can provide insight into viral pathogenesis. GWAS in difficult-to-recruit populations identifies variants associated with Lassa fever outcome and susceptibility at loci proximal to LIF, GRM7 and LARGE1.

Precision medicine in immune-mediated inflammatory diseases (IMIDs) requires an understanding of how cellular networks change following therapy. We describe a therapeutic atlas for Crohn’s disease (CD) and ulcerative colitis (UC) following anti-tumour necrosis factor (TNF) therapy. We generated ~1 million single-cell transcriptomes, organised into 109 cell states, from 216 gut biopsies from 38 patients and three controls, revealing disease- and therapy-specific differences. A systems-biology analysis identified distinct spatially-resolved cellular microenvironments: granuloma signatures in CD and interferon (IFN)-response signatures localising to T-cell aggregates and epithelial damage in CD and UC. Longitudinal comparisons demonstrated that disease progression in non-responders associated with myeloid and stromal cell perturbations in CD and increased multi-cellular IFN signalling in UC. IFN signalling was also observed in rheumatoid arthritis (RA) synovium with a lymphoid pathotype. Our therapeutic atlas informs drug positioning across IMIDs, and suggests a rationale for the use of janus kinase (JAK) inhibition following anti-TNF resistance.

While genome-wide association studies (GWAS) and positive selection scans identify genomic loci driving human phenotypic diversity, functional validation is required to discover the variant(s) responsible. We dissected the IVD gene locus-which encodes the isovaleryl-CoA dehydrogenase enzyme-implicated by selection statistics, multiple GWAS, and clinical genetics as important to function and fitness. We combined luciferase assays, CRISPR/Cas9 genome-editing, massively parallel reporter assays (MPRA), and a deletion tiling MPRA strategy across regulatory loci. We identified three regulatory variants, including an indel, that may underpin GWAS signals for pulmonary fibrosis and testosterone, and that are linked on a positively selected haplotype in the Japanese population. These regulatory variants exhibit synergistic and opposing effects on IVD expression experimentally. Alleles at these variants lie on a haplotype tagged by the variant most strongly associated with IVD expression and metabolites, but with no functional evidence itself. This work demonstrates how comprehensive functional investigation and multiple technologies are needed to discover the true genetic drivers of phenotypic diversity.

Ebola virus (EBOV) causes Ebola virus disease (EVD), marked by severe hemorrhagic fever; however, the mechanisms underlying the disease remain unclear. To assess the molecular basis of EVD across time, we performed RNA sequencing on 17 tissues from a natural history study of 21 rhesus monkeys, developing new methods to characterize host-pathogen dynamics. We identified alterations in host gene expression with previously unknown tissue-specific changes, including downregulation of genes related to tissue connectivity. EBOV was widely disseminated throughout the body; using a new, broadly applicable deconvolution method, we found that viral load correlated with increased monocyte presence. Patterns of viral variation between tissues differentiated primary infections from compartmentalized infections, and several variants impacted viral fitness in a EBOV/Kikwit minigenome system, suggesting that functionally significant variants can emerge during early infection. This comprehensive portrait of host-pathogen dynamics in EVD illuminates new features of pathogenesis and establishes resources to study other emerging pathogens.

BACKGROUND:Keratinocyte carcinoma (KC) is the most common type of nonmelanoma skin cancer. Currently, prophylactic treatment options are limited. Nonablative fractional lasers (NAFL) have received the Food and Drug Administration approval for the treatment of actinic damage; however, their role in KC prophylaxis is not known.OBJECTIVE:The aim of this study is to determine whether NAFL treatment is associated with a decrease in subsequent facial KC development.MATERIALS AND METHODS:A retrospective cohort study of patients with a history of facial KC treated at the Massachusetts General Hospital Dermatology Laser and Cosmetic Center between 2005 and 2021 was conducted.RESULTS:Forty-three NAFL-treated patients with a history of facial KC and 52 matched control subjects were included in the study. The rate of subsequent facial KC development was 20.9% in NAFL-treated patients and 40.4% in control subjects (RR 0.52, p = .049). Control subjects developed new facial KC significantly sooner than NAFL-treated patients (p = .033). When controlling for age, gender, and skin type, control subjects were more likely to develop new facial KC than NAFL-treated patients (hazard ratio 2.65, p = .0169).CONCLUSION:NAFL treatment was associated with a decreased risk of subsequent facial KC development and may have a benefit for KC prophylaxis.

ABSTRACT While genome-wide association studies (GWAS) and selection scans identify genomic loci driving human phenotypic diversity, functional validation is required to discover the variant(s) responsible. We dissected the IVD locus, implicated by selection statistics, multiple GWAS, and clinical genetics as important to function and fitness. We combined luciferase assays, CRISPR/Cas9 genome-editing, massively parallel reporter assays (MPRA), and bashing of regulatory loci. We identified three regulatory variants, including an indel, that may underpin GWAS signals for pulmonary fibrosis and testosterone, and that are linked on a positively selected haplotype in the Japanese population. These regulatory variants exhibit synergistic and opposing effects on IVD expression experimentally. Alleles at these variants lie on a haplotype tagged by the variant most strongly associated with IVD expression and metabolites, but with no functional evidence itself. This work demonstrates how comprehensive functional investigation and multiple technologies are needed to discover the true genetic drivers of phenotypic diversity.

Summary Visualizing two-dimensional (2D) embeddings (e.g. UMAP or tSNE) is a key step in interrogating single-cell RNA sequencing (scRNA-Seq) data. Subsequently, users typically iterate between programmatic analyses (e.g. clustering and differential expression) and visual exploration (e.g. coloring cells by interesting features) to uncover biological signals in the data. Interactive tools exist to facilitate visual exploration of embeddings such as performing differential expression on user-selected cells. However, the practical utility of these tools is limited because they don’t support rapid movement of data and results to and from the programming environments where the bulk of data analysis takes place, interrupting the iterative process. Here, we present the Single-cell Interactive Viewer (Sciviewer), a tool that overcomes this limitation by allowing interactive visual interrogation of embeddings from within Python. Beyond differential expression analysis of user-selected cells, Sciviewer implements a novel method to identify genes varying locally along any user-specified direction on the embedding. Sciviewer enables rapid and flexible iteration between interactive and programmatic modes of scRNA-Seq exploration, illustrating a useful approach for analyzing high-dimensional data. Availability and implementation Code and examples are provided at https://github.com/colabobio/sciviewer

While investigating a signal of adaptive evolution in humans at the gene LARGE, we encountered an intriguing finding by Dr. Stefan Kunz that the gene plays a critical role in Lassa virus binding and entry. This led us to pursue field work to test our hypothesis that natural selection acting on LARGE-detected in the Yoruba population of Nigeria-conferred resistance to Lassa Fever in some West African populations. As we delved further, we conjectured that the "emerging" nature of recently discovered diseases like Lassa fever is related to a newfound capacity for detection, rather than a novel viral presence, and that humans have in fact been exposed to the viruses that cause such diseases for much longer than previously suspected. Dr. Stefan Kunz's critical efforts not only laid the groundwork for this discovery, but also inspired and catalyzed a series of events that birthed Sentinel, an ambitious and large-scale pandemic prevention effort in West Africa. Sentinel aims to detect and characterize deadly pathogens before they spread across the globe, through implementation of its three fundamental pillars: Detect, Connect, and Empower. More specifically, Sentinel is designed to detect known and novel infections rapidly, connect and share information in real time to identify emerging threats, and empower the public health community to improve pandemic preparedness and response anywhere in the world. We are proud to dedicate this work to Stefan Kunz, and eagerly invite new collaborators, experts, and others to join us in our efforts.

The high morbidity and mortality associated with clinical cases of Lassa fever, together with the lack of licensed vaccines and limited and partially effective interventions, make Lassa virus (LASV) an important health concern in its regions of endemicity in West Africa. Previous infection with LASV protects from disease after subsequent exposure, providing a framework for designing vaccines to elicit similar protective immunity. Multiple major lineages of LASV circulate in West Africa, and therefore, ideal vaccine candidates should elicit immunity to all lineages. We therefore sought to identify common T cell epitopes between Lassa fever survivors from Sierra Leone and Nigeria, where distinct lineages circulate. We identified three such epitopes derived from highly conserved regions within LASV proteins. In this process, we also identified nine other T cell epitopes. These data should help in the design of an effective pan-LASV vaccine.

Lassa virus infects hundreds of thousands of people each year across rural West Africa, resulting in a high number of cases of Lassa fever (LF), a febrile disease associated with high morbidity and significant mortality. The lack of approved treatments or interventions underscores the need for an effective vaccine. At least four viral lineages circulate in defined regions throughout West Africa with substantial interlineage nucleotide and amino acid diversity. An effective vaccine should be designed to elicit Lassa virus specific humoral and cell mediated immunity across all lineages. Most current vaccine candidates use only lineage IV antigens encoded by Lassa viruses circulating around Sierra Leone, Liberia, and Guinea but not Nigeria where lineages I-III are found. As previous infection is known to protect against disease from subsequent exposure, we sought to determine whether LF survivors from Nigeria and Sierra Leone harbor memory T cells that respond to lineage IV antigens. Our results indicate a high degree of cross-reactivity of CD8+ T cells from Nigerian LF survivors to lineage IV antigens. In addition, we identified regions within the Lassa virus glycoprotein complex and nucleoprotein that contributed to these responses while T cell epitopes were not widely conserved across our study group. These data are important for current efforts to design effective and efficient vaccine candidates that can elicit protective immunity across all Lassa virus lineages.

Summary Ebola virus (EBOV) causes epidemics with high case fatality rates, yet remains understudied due to the challenge of experimentation in high-containment and outbreak settings. To better understand EBOV infection in vivo , we used single-cell transcriptomics and CyTOF-based single-cell protein quantification to characterize peripheral immune cell activity during EBOV infection in rhesus monkeys. We obtained 100,000 transcriptomes and 15,000,000 protein profiles, providing insight into pathogenesis. We find that immature, proliferative monocyte-lineage cells with reduced antigen presentation capacity replace conventional circulating monocyte subsets within days of infection, while lymphocytes upregulate apoptosis genes and decline in abundance. By quantifying viral RNA abundance in individual cells, we identify molecular determinants of tropism and examine temporal dynamics in viral and host gene expression. Within infected cells, we observe that EBOV down-regulates STAT1 mRNA and interferon signaling, and up-regulates putative pro-viral genes (e.g., DYNLL1 and HSPA5 ), nominating cellular pathways the virus manipulates for its replication. Overall, this study sheds light on EBOV tropism, replication dynamics, and elicited immune response, and provides a framework for characterizing interactions between hosts and emerging viruses in a maximum containment setting.

Ebola virus (EBOV) causes epidemics with high mortality yet remains understudied due to the challenge of experimentation in high-containment and outbreak settings. Here, we used single-cell transcriptomics and CyTOF-based single-cell protein quantification to characterize peripheral immune cells during EBOV infection in rhesus monkeys. We obtained 100,000 transcriptomes and 15,000,000 protein profiles, finding that immature, proliferative monocyte-lineage cells with reduced antigen-presentation capacity replace conventional monocyte subsets, while lymphocytes upregulate apoptosis genes and decline in abundance. By quantifying intracellular viral RNA, we identify molecular determinants of tropism among circulating immune cells and examine temporal dynamics in viral and host gene expression. Within infected cells, EBOV downregulates STAT1 mRNA and interferon signaling, and it upregulates putative pro-viral genes (e.g., DYNLL1 and HSPA5), nominating pathways the virus manipulates for its replication. This study sheds light on EBOV tropism, replication dynamics, and elicited immune response and provides a framework for characterizing host-virus interactions under maximum containment.

Many factors contribute to the selection of epitope-specific T cells. Often, candidate epitopes are generated through in silico prediction of various biological processes, mainly peptide generation and MHC binding, that are then tested in relevant biological assays. This is especially true of epitopes derived from BSL-4 pathogens or where relevant patient samples are difficult to obtain. Two previous studies identified CD8+ T cell epitopes from the Lassa virus glycoprotein through in silico prediction, experimental MHC binding assays, and epitope generation in HLA-A2 transgenic mice. Using samples from ten HLA-A*02:01 Lassa fever survivors, we tested whether these previously described epitope-specific CD8+ T cells were present and their relation to the broader Lassa virus-specific T cell response. Using overnight stimulation assays, we detected robust LASV-specific responses to the glycoprotein and nucleoprotein, but only one of the three epitopes (GP60–68)shown to be present and protective in mice made a substantial contribution to the overall LASV-specific response. Using a more sensitive proliferation assay, we detected the remaining two epitopes in some individuals at a very low frequency. Overall, this study shows the limitations of epitope discovery through in silico prediction, MHC binding and transgenic mouse models and highlights the complex nature of T cell selection during natural infection of humans.

Identifying gene expression programs underlying both cell-type identity and cellular activities (e.g. life-cycle processes, responses to environmental cues) is crucial for understanding the organization of cells and tissues. Although single-cell RNA-Seq (scRNA-Seq) can quantify transcripts in individual cells, each cell's expression profile may be a mixture of both types of programs, making them difficult to disentangle. Here, we benchmark and enhance the use of matrix factorization to solve this problem. We show with simulations that a method we call consensus non-negative matrix factorization (cNMF) accurately infers identity and activity programs, including their relative contributions in each cell. To illustrate the insights this approach enables, we apply it to published brain organoid and visual cortex scRNA-Seq datasets; cNMF refines cell types and identifies both expected (e.g. cell cycle and hypoxia) and novel activity programs, including programs that may underlie a neurosecretory phenotype and synaptogenesis.

The recent Ebola epidemic exemplified the importance of understanding and controlling emerging infections. Despite the importance of T cells in clearing virus during acute infection, little is known about Ebola-specific CD8(+) T cell responses. We investigated immune responses of individuals infected with Ebola virus (EBOV) during the 2013-2016 West Africa epidemic in Sierra Leone, where the majority of the >28,000 EBOV disease (EVD) cases occurred. We examined T cell memory responses to seven of the eight Ebola proteins (GP, sGP, NP, VP24, VP30, VP35, and VP40) and associated HLA expression in survivors. Of the 30 subjects included in our analysis, CD8(+) T cells from 26 survivors responded to at least one EBOV antigen. A minority, 10 of 26 responders (38%), made CD8(+) T cell responses to the viral GP or sGP. In contrast, 25 of the 26 responders (96%) made response to viral NP, 77% to VP24 (20 of 26), 69% to VP40 (18 of 26), 42% (11 of 26) to VP35, with no response to VP30. Individuals making CD8(+) T cells to EBOV VP24, VP35, and VP40 also made CD8(+) T cells to NP, but rarely to GP. We identified 34 CD8(+ )T cell epitopes for Ebola. Our data indicate the immunodominance of the EBOV NP-specific T cell response and suggest that its inclusion in a vaccine along with the EBOV GP would best mimic survivor responses and help boost cell-mediated immunity during vaccination.

Biological interpretation of genome-wide association study data frequently involves assessing whether SNPs linked to a biological process, for example, binding of a transcription factor, show unsigned enrichment for disease signal. However, signed annotations quantifying whether each SNP allele promotes or hinders the biological process can enable stronger statements about disease mechanism. We introduce a method, signed linkage disequilibrium profile regression, for detecting genome-wide directional effects of signed functional annotations on disease risk. We validate the method via simulations and application to molecular quantitative trait loci in blood, recovering known transcriptional regulators. We apply the method to expression quantitative trait loci in 48 Genotype-Tissue Expression tissues, identifying 651 transcription factor-tissue associations including 30 with robust evidence of tissue specificity. We apply the method to 46 diseases and complex traits (average n = 290 K), identifying 77 annotation-trait associations representing 12 independent transcription factor-trait associations, and characterize the underlying transcriptional programs using gene-set enrichment analyses. Our results implicate new causal disease genes and new disease mechanisms.

Lassa fever, a hemorrhagic fever caused by Lassa virus (LASV), is endemic in West Africa. It is difficult to distinguish febrile illnesses that are common in West Africa from Lassa fever based solely on a patient's clinical presentation. The field performance of recombinant antigen-based Lassa fever immunoassays was compared to that of quantitative polymerase chain assays (qPCRs) using samples from subjects meeting the case definition of Lassa fever presenting to Kenema Government Hospital in Sierra Leone. The recombinant Lassa virus (ReLASV) enzyme-linked immunosorbant assay (ELISA) for detection of viral antigen in blood performed with 95% sensitivity and 97% specificity using a diagnostic standard that combined results of the immunoassays and qPCR. The ReLASV rapid diagnostic test (RDT), a lateral flow immunoassay based on paired monoclonal antibodies to the Josiah strain of LASV (lineage IV), performed with 90% sensitivity and 100% specificity. ReLASV immunoassays performed better than the most robust qPCR currently available, which had 82% sensitivity and 95% specificity. The performance characteristics of recombinant antigen-based Lassa virus immunoassays indicate that they can aid in the diagnosis of LASV Infection and inform the clinical management of Lassa fever patients.