CD8+ lymphocytes are critical for early control of tuberculosis in macaques.

The functional role of CD8+ lymphocytes in tuberculosis remains poorly understood. We depleted innate and/or adaptive CD8+ lymphocytes in macaques and showed that loss of all CD8α+ cells (using anti-CD8α antibody) significantly impaired early control of Mycobacterium tuberculosis (Mtb) infection, leading to increased granulomas, lung inflammation, and bacterial burden. Analysis of barcoded Mtb from infected macaques demonstrated that depletion of all CD8+ lymphocytes allowed increased establishment of Mtb in lungs and dissemination within lungs and to lymph nodes, while depletion of only adaptive CD8+ T cells (with anti-CD8ß antibody) worsened bacterial control in lymph nodes. Flow cytometry and single-cell RNA sequencing revealed polyfunctional cytotoxic CD8+ lymphocytes in control granulomas, while CD8-depleted animals were unexpectedly enriched in CD4 and γδ T cells adopting incomplete cytotoxic signatures. Ligand-receptor analyses identified IL-15 signaling in granulomas as a driver of cytotoxic T cells. These data support that CD8+ lymphocytes are required for early protection against Mtb and suggest polyfunctional cytotoxic responses as a vaccine target.

Protective intravenous BCG vaccination induces enhanced immune signaling in the airways.

Intradermal (ID) Bacillus Calmette-Guérin (BCG) is the most widely administered vaccine in the world. However, ID-BCG fails to achieve the level of protection needed in adults to alter the course of the tuberculosis epidemic. Recent studies in non-human primates have demonstrated high levels of protection against Mycobacterium tuberculosis ( Mtb ) following intravenous (IV) administration of BCG. However, the protective immune features that emerge following IV BCG vaccination remain incompletely defined. Here we used single-cell RNA-sequencing (scRNAseq) to transcriptionally profile 157,114 unstimulated and purified protein derivative (PPD)-stimulated bronchoalveolar lavage (BAL) cells from 29 rhesus macaques immunized with BCG across routes of administration and doses to uncover cell composition-, gene expression-, and biological network-level signatures associated with IV BCG-mediated protection. Our analyses revealed that high-dose IV BCG drove an influx of polyfunctional T cells and macrophages into the airways. These macrophages exhibited a basal activation phenotype even in the absence of PPD-stimulation, defined in part by IFN and TNF-α signaling up to 6 months following BCG immunization. Furthermore, intercellular immune signaling pathways between key myeloid and T cell subsets were enhanced following PPD-stimulation in high-dose IV BCG-vaccinated macaques. High-dose IV BCG also engendered quantitatively and qualitatively stronger transcriptional responses to PPD-stimulation, with a robust Th1-Th17 transcriptional phenotype in T cells, and augmented transcriptional signatures of reactive oxygen species production, hypoxia, and IFN-γ response within alveolar macrophages. Collectively, this work supports that IV BCG immunization creates a unique cellular ecosystem in the airways, which primes and enables local myeloid cells to effectively clear Mtb upon challenge.

CD4+ T cells are homeostatic regulators during Mtb reinfection.

Immunological priming - either in the context of prior infection or vaccination - elicits protective responses against subsequent Mycobacterium tuberculosis (Mtb) infection. However, the changes that occur in the lung cellular milieu post-primary Mtb infection and their contributions to protection upon reinfection remain poorly understood. Here, using clinical and microbiological endpoints in a non-human primate reinfection model, we demonstrate that prior Mtb infection elicits a long-lasting protective response against subsequent Mtb exposure and that the depletion of CD4+ T cells prior to Mtb rechallenge significantly abrogates this protection. Leveraging microbiologic, PET-CT, flow cytometric, and single-cell RNA-seq data from primary infection, reinfection, and reinfection-CD4+ T cell depleted granulomas, we identify differential cellular and microbial features of control. The data collectively demonstrate that the presence of CD4+ T cells in the setting of reinfection results in a reduced inflammatory lung milieu characterized by reprogrammed CD8+ T cell activity, reduced neutrophilia, and blunted type-1 immune signaling among myeloid cells, mitigating Mtb disease severity. These results open avenues for developing vaccines and therapeutics that not only target CD4+ and CD8+ T cells, but also modulate innate immune cells to limit Mtb disease.

HIV specific CD8+ TRM-like cells in tonsils express exhaustive signatures in the absence of natural HIV control.

Lymphoid tissues are an important HIV reservoir site that persists in the face of antiretroviral therapy and natural immunity. Targeting these reservoirs by harnessing the antiviral activity of local tissue-resident memory (TRM) CD8+ T-cells is of great interest, but limited data exist on TRM-like cells within lymph nodes of people living with HIV (PLWH). Here, we studied tonsil CD8+ T-cells obtained from PLWH and uninfected controls from South Africa. We show that these cells are preferentially located outside the germinal centers (GCs), the main reservoir site for HIV, and display a low cytolytic and a transcriptionally TRM-like profile distinct from blood CD8+ T-cells. In PLWH, CD8+ TRM-like cells are expanded and adopt a more cytolytic, activated, and exhausted phenotype not reversed by antiretroviral therapy (ART). This phenotype was enhanced in HIV-specific CD8+ T-cells from tonsils compared to matched blood suggesting a higher antigen burden in tonsils. Single-cell transcriptional and clonotype resolution showed that these HIV-specific CD8+ T-cells in the tonsils express heterogeneous signatures of T-cell activation, clonal expansion, and exhaustion ex-vivo. Interestingly, this signature was absent in a natural HIV controller, who expressed lower PD-1 and CXCR5 levels and reduced transcriptional evidence of T-cell activation, exhaustion, and cytolytic activity. These data provide important insights into lymphoid tissue-derived HIV-specific CD8+ TRM-like phenotypes in settings of HIV remission and highlight their potential for immunotherapy and targeting of the HIV reservoirs.

Identification and characterization of the T cell receptor (TCR) repertoire of the cynomolgus macaque (Macaca Fascicularis).

BACKGROUND: Cynomolgus macaque (Macaca fascicularis) is an attractive animal model for the study of human disease and is extensively used in biomedical research. Cynomolgus macaques share behavioral, physiological, and genomic traits with humans and recapitulate human disease manifestations not observed in other animal species. To improve the use of the cynomolgus macaque model to investigate immune responses, we defined and characterized the T cell receptor (TCR) repertoire. RESULT: We identified and analyzed the alpha (TRA), beta (TRB), gamma (TRG), and delta (TRD) TCR loci of the cynomolgus macaque. The expressed repertoire was determined using 22 unique lung samples from Mycobacterium tuberculosis infected cynomolgus macaques by single cell RNA sequencing. Expressed TCR alpha (TRAV) and beta (TRBV) variable region genes were enriched and identified using gene specific primers, which allowed their functional status to be determined. Analysis of the primers used for cynomolgus macaque TCR variable region gene enrichment showed they could also be used to amplify rhesus macaque (M. mulatta) variable region genes. CONCLUSION: The genomic organization of the cynomolgus macaque has great similarity with the rhesus macaque and they shared > 90% sequence similarity with the human TCR repertoire. The identification of the TCR repertoire facilitates analysis of T cell immunity in cynomolgus macaques.

Enriching and Characterizing T Cell Repertoires from 3' Barcoded Single-Cell Whole Transcriptome Amplification Products.

Antigen-specific T cells play an essential role in immunoregulation and many diseases such as cancer. Characterizing the T cell receptor (TCR) sequences that encode T cell specificity is critical for elucidating the antigenic determinants of immunological diseases and designing therapeutic remedies. However, methods of obtaining single-cell TCR sequencing data are labor and cost intensive, typically requiring both cell sorting and full-length single-cell RNA-sequencing (scRNA-seq). New high-throughput 3' cell-barcoding scRNA-seq methods can simplify and scale this process; however, they do not routinely capture TCR sequences during library preparation and sequencing. While 5' cell-barcoding scRNA-seq methods can be used to examine TCR repertoire at single-cell resolution, doing so requires specialized reagents which cannot be applied to samples previously processed using 3' cell-barcoding methods.Here, we outline a method for sequencing TCRα and TCRß transcripts from samples already processed using 3' cell-barcoding scRNA-seq platforms, ensuring TCR recovery at a single-cell resolution. In short, a fraction of the 3' barcoded whole transcriptome amplification (WTA) product typically used to generate a massively parallel 3' scRNA-seq library is enriched for TCR transcripts using biotinylated probes and further amplified using the same universal primer sequence from WTA. Primer extension using TCR V-region primers and targeted PCR amplification using a second universal primer result in a 3' barcoded single-cell CDR3-enriched library that can be sequenced with custom sequencing primers. Coupled with 3' scRNA-seq of the same WTA, this method enables simultaneous analysis of single-cell transcriptomes and TCR sequences which can help interpret inherent heterogeneity among antigen-specific T cells and salient disease biology. The method presented here can also be adapted readily to enrich and sequence other transcripts of interest from both 3' and 5' barcoded scRNA-seq WTA libraries.

Multimodal profiling of lung granulomas in macaques reveals cellular correlates of tuberculosis control.

Mycobacterium tuberculosis lung infection results in a complex multicellular structure: the granuloma. In some granulomas, immune activity promotes bacterial clearance, but in others, bacteria persist and grow. We identified correlates of bacterial control in cynomolgus macaque lung granulomas by co-registering longitudinal positron emission tomography and computed tomography imaging, single-cell RNA sequencing, and measures of bacterial clearance. Bacterial persistence occurred in granulomas enriched for mast, endothelial, fibroblast, and plasma cells, signaling amongst themselves via type 2 immunity and wound-healing pathways. Granulomas that drove bacterial control were characterized by cellular ecosystems enriched for type 1-type 17, stem-like, and cytotoxic T cells engaged in pro-inflammatory signaling networks involving diverse cell populations. Granulomas that arose later in infection displayed functional characteristics of restrictive granulomas and were more capable of killing Mtb. Our results define the complex multicellular ecosystems underlying (lack of) granuloma resolution and highlight host immune targets that can be leveraged to develop new vaccine and therapeutic strategies for TB.

Cellular and transcriptional diversity over the course of human lactation.

Human breast milk (hBM) is a dynamic fluid that contains millions of cells, but their identities and phenotypic properties are poorly understood. We generated and analyzed single-cell RNA-sequencing (scRNA-seq) data to characterize the transcriptomes of cells from hBM across lactational time from 3 to 632 d postpartum in 15 donors. We found that the majority of cells in hBM are lactocytes, a specialized epithelial subset, and that cell-type frequencies shift over the course of lactation, yielding greater epithelial diversity at later points. Analysis of lactocytes reveals a continuum of cell states characterized by transcriptional changes in hormone-, growth factor-, and milk production-related pathways. Generalized additive models suggest that one subcluster, LC1 epithelial cells, increases as a function of time postpartum, daycare attendance, and the use of hormonal birth control. We identify several subclusters of macrophages in hBM that are enriched for tolerogenic functions, possibly playing a role in protecting the mammary gland during lactation. Our description of the cellular components of breast milk, their association with maternal­infant dyad metadata, and our quantification of alterations at the gene and pathway levels provide a detailed longitudinal picture of hBM cells across lactational time. This work paves the way for future investigations of how a potential division of cellular labor and differential hormone regulation might be leveraged therapeutically to support healthy lactation and potentially aid in milk production.

HIV infection drives interferon signaling within intestinal SARS-CoV-2 target cells.

SARS-CoV-2 infects epithelial cells of the human gastrointestinal (GI) tract and causes related symptoms. HIV infection impairs gut homeostasis and is associated with an increased risk of COVID-19 fatality. To investigate the potential link between these observations, we analyzed single-cell transcriptional profiles and SARS-CoV-2 entry receptor expression across lymphoid and mucosal human tissue from chronically HIV-infected individuals and uninfected controls. Absorptive gut enterocytes displayed the highest coexpression of SARS-CoV-2 receptors ACE2, TMPRSS2, and TMPRSS4, of which ACE2 expression was associated with canonical interferon response and antiviral genes. Chronic treated HIV infection was associated with a clear antiviral response in gut enterocytes and, unexpectedly, with a substantial reduction of ACE2 and TMPRSS2 target cells. Gut tissue from SARS-CoV-2-infected individuals, however, showed abundant SARS-CoV-2 nucleocapsid protein in both the large and small intestine, including an HIV-coinfected individual. Thus, upregulation of antiviral response genes and downregulation of ACE2 and TMPRSS2 in the GI tract of HIV-infected individuals does not prevent SARS-CoV-2 infection in this compartment. The impact of these HIV-associated intestinal mucosal changes on SARS-CoV-2 infection dynamics, disease severity, and vaccine responses remains unclear and requires further investigation.

Human airway mast cells proliferate and acquire distinct inflammation-driven phenotypes during type 2 inflammation.

Mast cells (MCs) play a pathobiologic role in type 2 (T2) allergic inflammatory diseases of the airway, including asthma and chronic rhinosinusitis with nasal polyposis (CRSwNP). Distinct MC subsets infiltrate the airway mucosa in T2 disease, including subepithelial MCs expressing the proteases tryptase and chymase (MCTC) and epithelial MCs expressing tryptase without chymase (MCT). However, mechanisms underlying MC expansion and the transcriptional programs underlying their heterogeneity are poorly understood. Here, we use flow cytometry and single-cell RNA-sequencing (scRNA-seq) to conduct a comprehensive analysis of human MC hyperplasia in CRSwNP, a T2 cytokine-mediated inflammatory disease. We link discrete cell surface phenotypes to the distinct transcriptomes of CRSwNP MCT and MCTC, which represent polarized ends of a transcriptional gradient of nasal polyp MCs. We find a subepithelial population of CD38highCD117high MCs that is markedly expanded during T2 inflammation. These CD38highCD117high MCs exhibit an intermediate phenotype relative to the expanded MCT and MCTC subsets. CD38highCD117high MCs are distinct from circulating MC progenitors and are enriched for proliferation, which is markedly increased in CRSwNP patients with aspirin-exacerbated respiratory disease, a severe disease subset characterized by increased MC burden and elevated MC activation. We observe that MCs expressing a polyp MCT-like effector program are also found within the lung during fibrotic diseases and asthma, and further identify marked differences between MCTC in nasal polyps and skin. These results indicate that MCs display distinct inflammation-associated effector programs and suggest that in situ MC proliferation is a major component of MC hyperplasia in human T2 inflammation.

Improved haplotype inference by exploiting long-range linking and allelic imbalance in RNA-seq datasets.

Haplotype reconstruction of distant genetic variants remains an unsolved problem due to the short-read length of common sequencing data. Here, we introduce HapTree-X, a probabilistic framework that utilizes latent long-range information to reconstruct unspecified haplotypes in diploid and polyploid organisms. It introduces the observation that differential allele-specific expression can link genetic variants from the same physical chromosome, thus even enabling using reads that cover only individual variants. We demonstrate HapTree-X's feasibility on in-house sequenced Genome in a Bottle RNA-seq and various whole exome, genome, and 10X Genomics datasets. HapTree-X produces more complete phases (up to 25%), even in clinically important genes, and phases more variants than other methods while maintaining similar or higher accuracy and being up to 10×  faster than other tools. The advantage of HapTree-X's ability to use multiple lines of evidence, as well as to phase polyploid genomes in a single integrative framework, substantially grows as the amount of diverse data increases.

SARS-CoV-2 Receptor ACE2 Is an Interferon-Stimulated Gene in Human Airway Epithelial Cells and Is Detected in Specific Cell Subsets across Tissues.

There is pressing urgency to understand the pathogenesis of the severe acute respiratory syndrome coronavirus clade 2 (SARS-CoV-2), which causes the disease COVID-19. SARS-CoV-2 spike (S) protein binds angiotensin-converting enzyme 2 (ACE2), and in concert with host proteases, principally transmembrane serine protease 2 (TMPRSS2), promotes cellular entry. The cell subsets targeted by SARS-CoV-2 in host tissues and the factors that regulate ACE2 expression remain unknown. Here, we leverage human, non-human primate, and mouse single-cell RNA-sequencing (scRNA-seq) datasets across health and disease to uncover putative targets of SARS-CoV-2 among tissue-resident cell subsets. We identify ACE2 and TMPRSS2 co-expressing cells within lung type II pneumocytes, ileal absorptive enterocytes, and nasal goblet secretory cells. Strikingly, we discovered that ACE2 is a human interferon-stimulated gene (ISG) in vitro using airway epithelial cells and extend our findings to in vivo viral infections. Our data suggest that SARS-CoV-2 could exploit species-specific interferon-driven upregulation of ACE2, a tissue-protective mediator during lung injury, to enhance infection.

Integrated single-cell analysis of multicellular immune dynamics during hyperacute HIV-1 infection.

Cellular immunity is critical for controlling intracellular pathogens, but individual cellular dynamics and cell-cell cooperativity in evolving human immune responses remain poorly understood. Single-cell RNA-sequencing (scRNA-seq) represents a powerful tool for dissecting complex multicellular behaviors in health and disease1,2 and nominating testable therapeutic targets3. Its application to longitudinal samples could afford an opportunity to uncover cellular factors associated with the evolution of disease progression without potentially confounding inter-individual variability4. Here, we present an experimental and computational methodology that uses scRNA-seq to characterize dynamic cellular programs and their molecular drivers, and apply it to HIV infection. By performing scRNA-seq on peripheral blood mononuclear cells from four untreated individuals before and longitudinally during acute infection5, we were powered within each to discover gene response modules that vary by time and cell subset. Beyond previously unappreciated individual- and cell-type-specific interferon-stimulated gene upregulation, we describe temporally aligned gene expression responses obscured in bulk analyses, including those involved in proinflammatory T cell differentiation, prolonged monocyte major histocompatibility complex II upregulation and persistent natural killer (NK) cell cytolytic killing. We further identify response features arising in the first weeks of infection, for example proliferating natural killer cells, which potentially may associate with future viral control. Overall, our approach provides a unified framework for characterizing multiple dynamic cellular responses and their coordination.

Allergic inflammatory memory in human respiratory epithelial progenitor cells.

Barrier tissue dysfunction is a fundamental feature of chronic human inflammatory diseases1. Specialized subsets of epithelial cells-including secretory and ciliated cells-differentiate from basal stem cells to collectively protect the upper airway2-4. Allergic inflammation can develop from persistent activation5 of type 2 immunity6 in the upper airway, resulting in chronic rhinosinusitis, which ranges in severity from rhinitis to severe nasal polyps7. Basal cell hyperplasia is a hallmark of severe disease7-9, but it is not known how these progenitor cells2,10,11 contribute to clinical presentation and barrier tissue dysfunction in humans. Here we profile primary human surgical chronic rhinosinusitis samples (18,036 cells, n = 12) that span the disease spectrum using Seq-Well for massively parallel single-cell RNA sequencing12, report transcriptomes for human respiratory epithelial, immune and stromal cell types and subsets from a type 2 inflammatory disease, and map key mediators. By comparison with nasal scrapings (18,704 cells, n = 9), we define signatures of core, healthy, inflamed and polyp secretory cells. We reveal marked differences between the epithelial compartments of the non-polyp and polyp cellular ecosystems, identifying and validating a global reduction in cellular diversity of polyps characterized by basal cell hyperplasia, concomitant decreases in glandular cells, and phenotypic shifts in secretory cell antimicrobial expression. We detect an aberrant basal progenitor differentiation trajectory in polyps, and propose cell-intrinsic13, epigenetic14,15 and extrinsic factors11,16,17 that lock polyp basal cells into this uncommitted state. Finally, we functionally demonstrate that ex vivo cultured basal cells retain intrinsic memory of IL-4/IL-13 exposure, and test the potential for clinical blockade of the IL-4 receptor α-subunit to modify basal and secretory cell states in vivo. Overall, we find that reduced epithelial diversity stemming from functional shifts in basal cells is a key characteristic of type 2 immune-mediated barrier tissue dysfunction. Our results demonstrate that epithelial stem cells may contribute to the persistence of human disease by serving as repositories for allergic memories.

De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds.

The Zika outbreak, spread by the Aedes aegypti mosquito, highlights the need to create high-quality assemblies of large genomes in a rapid and cost-effective way. Here we combine Hi-C data with existing draft assemblies to generate chromosome-length scaffolds. We validate this method by assembling a human genome, de novo, from short reads alone (67× coverage). We then combine our method with draft sequences to create genome assemblies of the mosquito disease vectors Aeaegypti and Culex quinquefasciatus, each consisting of three scaffolds corresponding to the three chromosomes in each species. These assemblies indicate that almost all genomic rearrangements among these species occur within, rather than between, chromosome arms. The genome assembly procedure we describe is fast, inexpensive, and accurate, and can be applied to many species.