Single-Cell and Bulk RNA-Sequencing Reveal Differences in Monocyte Susceptibility to Influenza A Virus Infection Between Africans and Europeans

Shaping immunity: The influence of natural selection on population immune diversity

Humans exhibit considerable variability in their immune responses to the same immune challenges. Such variation is widespread and affects individual and population-level susceptibility to infectious diseases and immune disorders. Although the factors influencing immune response diversity are partially understood, what mechanisms lead to the wide range of immune traits in healthy individuals remain largely unexplained. Here, we discuss the role that natural selection has played in driving phenotypic differences in immune responses across populations and present-day susceptibility to immune-related disorders. Further, we touch on future directions in the field of immunogenomics, highlighting the value of expanding this work to human populations globally, the utility of modeling the immune response as a dynamic process, and the importance of considering the potential polygenic nature of natural selection. Identifying loci acted upon by evolution may further pinpoint variants critically involved in disease etiology, and designing studies to capture these effects will enrich our understanding of the genetic contributions to immunity and immune dysregulation.

Variation in the basal immune state and implications for disease

Analysis of pre-existing immunity and its effects on acute infection often focus on memory responses associated with a prior infectious exposure. However, memory responses occur in the context of the overall immune state and leukocytes must interact with their microenvironment and other immune cells. Thus, it is important to also consider non-antigen-specific factors which shape the composite basal state and functional capacity of the immune system, termed here as I0 ('I naught'). In this review, we discuss the determinants of I0. Utilizing influenza virus as a model, we then consider the effect of I0 on susceptibility to infection and disease severity. Lastly, we outline a mathematical framework and demonstrate how researchers can build and tailor models to specific needs. Understanding how diverse factors uniquely and collectively impact immune competence will provide valuable insights into mechanisms of immune variation, aid in screening for high-risk populations, and promote the development of broadly applicable prophylactic and therapeutic treatments.

Dissecting human population variation in single-cell responses to SARS-CoV-2

Humans display substantial interindividual clinical variability after SARS-CoV-2 infection1-3, the genetic and immunological basis of which has begun to be deciphered4. However, the extent and drivers of population differences in immune responses to SARS-CoV-2 remain unclear. Here we report single-cell RNA-sequencing data for peripheral blood mononuclear cells-from 222 healthy donors of diverse ancestries-that were stimulated with SARS-CoV-2 or influenza A virus. We show that SARS-CoV-2 induces weaker, but more heterogeneous, interferon-stimulated gene activity compared with influenza A virus, and a unique pro-inflammatory signature in myeloid cells. Transcriptional responses to viruses display marked population differences, primarily driven by changes in cell abundance including increased lymphoid differentiation associated with latent cytomegalovirus infection. Expression quantitative trait loci and mediation analyses reveal a broad effect of cell composition on population disparities in immune responses, with genetic variants exerting a strong effect on specific loci. Furthermore, we show that natural selection has increased population differences in immune responses, particularly for variants associated with SARS-CoV-2 response in East Asians, and document the cellular and molecular mechanisms by which Neanderthal introgression has altered immune functions, such as the response of myeloid cells to viruses. Finally, colocalization and transcriptome-wide association analyses reveal an overlap between the genetic basis of immune responses to SARS-CoV-2 and COVID-19 severity, providing insights into the factors contributing to current disparities in COVID-19 risk.

Baseline innate and T cell populations are correlates of protection against symptomatic influenza virus infection independent of serology

Evidence suggests that innate and adaptive cellular responses mediate resistance to the influenza virus and confer protection after vaccination. However, few studies have resolved the contribution of cellular responses within the context of preexisting antibody titers. Here, we measured the peripheral immune profiles of 206 vaccinated or unvaccinated adults to determine how baseline variations in the cellular and humoral immune compartments contribute independently or synergistically to the risk of developing symptomatic influenza. Protection correlated with diverse and polyfunctional CD4+ and CD8+ T, circulating T follicular helper, T helper type 17, myeloid dendritic and CD16+ natural killer (NK) cell subsets. Conversely, increased susceptibility was predominantly attributed to nonspecific inflammatory populations, including γδ T cells and activated CD16- NK cells, as well as TNFα+ single-cytokine-producing CD8+ T cells. Multivariate and predictive modeling indicated that cellular subsets (1) work synergistically with humoral immunity to confer protection, (2) improve model performance over demographic and serologic factors alone and (3) comprise the most important predictive covariates. Together, these results demonstrate that preinfection peripheral cell composition improves the prediction of symptomatic influenza susceptibility over vaccination, demographics or serology alone.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Transposable elements are associated with the variable response to influenza infection

Influenza A virus (IAV) infections are frequent every year and result in a range of disease severity. Here, we wanted to explore the potential contribution of transposable elements (TEs) to the variable human immune response. Transcriptome profiling in monocyte-derived macrophages from 39 individuals following IAV infection revealed significant inter-individual variation in viral load post-infection. Using transposase-accessible chromatin using sequencing (ATAC-seq), we identified a set of TE families with either enhanced or reduced accessibility upon infection. Of the enhanced families, 15 showed high variability between individuals and had distinct epigenetic profiles. Motif analysis showed an association with known immune regulators (e.g., BATFs, FOSs/JUNs, IRFs, STATs, NFkBs, NFYs, and RELs) in stably enriched families and with other factors in variable families, including KRAB-ZNFs. We showed that TEs and host factors regulating TEs were predictive of viral load post-infection. Our findings shed light on the role TEs and KRAB-ZNFs may play in inter-individual variation in immunity.

Integrated Drivers of Basal and Acute Immunity in Diverse Human Populations

Prior studies have identified genetic, infectious, and biological associations with immune competence and disease severity; however, there have been few integrative analyses of these factors and study populations are often limited in demographic diversity. Utilizing samples from 1,705 individuals in 5 countries, we examined putative determinants of immunity, including: single nucleotide polymorphisms, ancestry informative markers, herpesvirus status, age, and sex. In healthy subjects, we found significant differences in cytokine levels, leukocyte phenotypes, and gene expression. Transcriptional responses also varied by cohort, and the most significant determinant was ancestry. In influenza infected subjects, we found two disease severity immunophenotypes, largely driven by age. Additionally, cytokine regression models show each determinant differentially contributes to acute immune variation, with unique and interactive, location-specific herpesvirus effects. These results provide novel insight into the scope of immune heterogeneity across diverse populations, the integrative effects of factors which drive it, and the consequences for illness outcomes.

High-throughput droplet-based analysis of influenza A virus genetic reassortment by single-virus RNA sequencing

The segmented RNA genome of influenza A viruses (IAVs) enables viral evolution through genetic reassortment after multiple IAVs coinfect the same cell, leading to viruses harboring combinations of eight genomic segments from distinct parental viruses. Existing data indicate that reassortant genotypes are not equiprobable; however, the low throughput of available virology techniques does not allow quantitative analysis. Here, we have developed a high-throughput single-cell droplet microfluidic system allowing encapsulation of IAV-infected cells, each cell being infected by a single progeny virion resulting from a coinfection process. Customized barcoded primers for targeted viral RNA sequencing enabled the analysis of 18,422 viral genotypes resulting from coinfection with two circulating human H1N1pdm09 and H3N2 IAVs. Results were highly reproducible, confirmed that genetic reassortment is far from random, and allowed accurate quantification of reassortants including rare events. In total, 159 out of the 254 possible reassortant genotypes were observed but with widely varied prevalence (from 0.038 to 8.45%). In cells where eight segments were detected, all 112 possible pairwise combinations of segments were observed. The inclusion of data from single cells where less than eight segments were detected allowed analysis of pairwise cosegregation between segments with very high confidence. Direct coupling analysis accurately predicted the fraction of pairwise segments and full genotypes. Overall, our results indicate that a large proportion of reassortant genotypes can emerge upon coinfection and be detected over a wide range of frequencies, highlighting the power of our tool for systematic and exhaustive monitoring of the reassortment potential of IAVs.

The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.

Single-cell genome-wide association reveals that a nonsynonymous variant in ERAP1 confers increased susceptibility to influenza virus

During pandemics, individuals exhibit differences in risk and clinical outcomes. Here, we developed single-cell high-throughput human in vitro susceptibility testing (scHi-HOST), a method for rapidly identifying genetic variants that confer resistance and susceptibility. We applied this method to influenza A virus (IAV), the cause of four pandemics since the start of the 20th century. scHi-HOST leverages single-cell RNA sequencing (scRNA-seq) to simultaneously assign genetic identity to cells in mixed infections of cell lines of European, African, and Asian origin, reveal associated genetic variants for viral burden, and identify expression quantitative trait loci. Integration of scHi-HOST with human challenge and experimental validation demonstrated that a missense variant in endoplasmic reticulum aminopeptidase 1 (ERAP1; rs27895) increased IAV burden in cells and human volunteers. rs27895 exhibits population differentiation, likely contributing to greater permissivity of cells from African populations to IAV. scHi-HOST is a broadly applicable method and resource for decoding infectious-disease genetics.

In Vivo Profiling of Individual Multiciliated Cells during Acute Influenza A Virus Infection

Influenza virus infections are thought to be initiated in a small number of cells; however, the heterogeneity across the cellular responses of the epithelial cells during establishment of disease is incompletely understood. Here, we used an H1N1 influenza virus encoding a fluorescent reporter gene, a cell lineage-labeling transgenic mouse line, and single-cell RNA sequencing to explore the range of responses in a susceptible epithelial cell population during an acute influenza A virus (IAV) infection. Focusing on multiciliated cells, we identified a subpopulation that basally expresses interferon-stimulated genes (ISGs), which we hypothesize may be important for the early response to infection. We subsequently found that a population of infected ciliated cells produce most of the ciliated cell-derived inflammatory cytokines, and nearly all bystander ciliated cells induce a broadly antiviral state. From these data together, we propose that variable preexisting gene expression patterns in the initial cells targeted by the virus may ultimately affect the establishment of viral disease. IMPORTANCE Influenza A virus poses a significant threat to public health, and each year, millions of people in the United States alone are exposed to the virus. We do not currently, however, fully understand why some individuals clear the infection asymptomatically and others become severely ill. Understanding how these divergent phenotypes arise could eventually be leveraged to design therapeutics that prevent severe disease. As a first step toward understanding these different infection states, we used a technology that allowed us to determine how thousands of individual murine lung epithelial cells behaved before and during IAV infection. We found that small subsets of epithelial cells exhibited an antiviral state prior to infection, and similarly, some cells made high levels of inflammatory cytokines during infection. We propose that different ratios of these individual cellular responses may contribute to the broader antiviral state of the lung and may ultimately affect disease severity.