An esophagus cell atlas reveals dynamic rewiring during active eosinophilic esophagitis and remission

Coordinated cell interactions within the esophagus maintain homeostasis, and disruption can lead to eosinophilic esophagitis (EoE), a chronic inflammatory disease with poorly understood pathogenesis. We profile 421,312 individual cells from the esophageal mucosa of 7 healthy and 15 EoE participants, revealing 60 cell subsets and functional alterations in cell states, compositions, and interactions that highlight previously unclear features of EoE. Active disease displays enrichment of ALOX15+ macrophages, PRDM16+ dendritic cells expressing the EoE risk gene ATP10A, and cycling mast cells, with concomitant reduction of TH17 cells. Ligand-receptor expression uncovers eosinophil recruitment programs, increased fibroblast interactions in disease, and IL-9+IL-4+IL-13+ TH2 and endothelial cells as potential mast cell interactors. Resolution of inflammation-associated signatures includes mast and CD4+ TRM cell contraction and cell type-specific downregulation of eosinophil chemoattractant, growth, and survival factors. These cellular alterations in EoE and remission advance our understanding of eosinophilic inflammation and opportunities for therapeutic intervention.

Multi-modal skin atlas identifies a multicellular immune-stromal community associated with altered cornification and specific T cell expansion in atopic dermatitis

In healthy skin, a cutaneous immune system maintains the balance between tolerance towards innocuous environmental antigens and immune responses against pathological agents. In atopic dermatitis (AD), barrier and immune dysfunction result in chronic tissue inflammation. Our understanding of the skin tissue ecosystem in AD remains incomplete with regard to the hallmarks of pathological barrier formation, and cellular state and clonal composition of disease-promoting cells. Here, we generated a multi-modal cell census of 310,691 cells spanning 86 cell subsets from whole skin tissue of 19 adult individuals, including non-lesional and lesional skin from 11 AD patients, and integrated it with 396,321 cells from four studies into a comprehensive human skin cell atlas in health and disease. Reconstruction of human keratinocyte differentiation from basal to cornified layers revealed a disrupted cornification trajectory in AD. This disrupted epithelial differentiation was associated with signals from a unique immune and stromal multicellular community comprised of MMP12 + dendritic cells (DCs), mature migratory DCs, cycling ILCs, NK cells, inflammatory CCL19 + IL4I1 + fibroblasts, and clonally expanded IL13 + IL22 + IL26 + T cells with overlapping type 2 and type 17 characteristics. Cell subsets within this immune and stromal multicellular community were connected by multiple inter-cellular positive feedback loops predicted to impact community assembly and maintenance. AD GWAS gene expression was enriched both in disrupted cornified keratinocytes and in cell subsets from the lesional immune and stromal multicellular community including IL13 + IL22 + IL26 + T cells and ILCs, suggesting that epithelial or immune dysfunction in the context of the observed cellular communication network can initiate and then converge towards AD. Our work highlights specific, disease-associated cell subsets and interactions as potential targets in progression and resolution of chronic inflammation.

Remodeling of colon plasma cell repertoire within ulcerative colitis patients

Plasma cells (PCs) constitute a significant fraction of colonic mucosal cells and contribute to inflammatory infiltrates in ulcerative colitis (UC). While gut PCs secrete bacteria-targeting IgA antibodies, their role in UC pathogenesis is unknown. We performed single-cell V(D)J- and RNA-seq on sorted B cells from the colon of healthy individuals and patients with UC. A large fraction of B cell clones is shared between different colon regions, but inflammation in UC broadly disrupts this landscape, causing transcriptomic changes characterized by an increase in the unfolded protein response (UPR) and antigen presentation genes, clonal expansion, and isotype skewing from IgA1 and IgA2 to IgG1. We also directly expressed and assessed the specificity of 152 mAbs from expanded PC clones. These mAbs show low polyreactivity and autoreactivity and instead target both shared bacterial antigens and specific bacterial strains. Altogether, our results characterize the microbiome-specific colon PC response and how its disruption might contribute to inflammation in UC.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

The CD4+ T cell response to a commensal-derived epitope transitions from a tolerant to an inflammatory state in Crohn's disease

Reciprocal interactions between host T helper cells and gut microbiota enforce local immunological tolerance and modulate extra-intestinal immunity. However, our understanding of antigen-specific tolerance to the microbiome is limited. Here, we developed a systematic approach to predict HLA class-II-specific epitopes using the humanized bacteria-originated T cell antigen (hBOTA) algorithm. We identified a diverse set of microbiome epitopes spanning all major taxa that are compatible with presentation by multiple HLA-II alleles. In particular, we uncovered an immunodominant epitope from the TonB-dependent receptor SusC that was universally recognized and ubiquitous among Bacteroidales. In healthy human subjects, SusC-reactive T cell responses were characterized by IL-10-dominant cytokine profiles, whereas in patients with active Crohn's disease, responses were associated with elevated IL-17A. Our results highlight the potential of targeted antigen discovery within the microbiome to reveal principles of tolerance and functional transitions during inflammation.

Stepwise-edited, human melanoma models reveal mutations' effect on tumor and microenvironment

Establishing causal relationships between genetic alterations of human cancers and specific phenotypes of malignancy remains a challenge. We sequentially introduced mutations into healthy human melanocytes in up to five genes spanning six commonly disrupted melanoma pathways, forming nine genetically distinct cellular models of melanoma. We connected mutant melanocyte genotypes to malignant cell expression programs in vitro and in vivo, replicative immortality, malignancy, rapid tumor growth, pigmentation, metastasis, and histopathology. Mutations in malignant cells also affected tumor microenvironment composition and cell states. Our melanoma models shared genotype-associated expression programs with patient melanomas, and a deep learning model showed that these models partially recapitulated genotype-associated histopathological features as well. Thus, a progressive series of genome-edited human cancer models can causally connect genotypes carrying multiple mutations to phenotype.

Stem-like intestinal Th17 cells give rise to pathogenic effector T cells during autoimmunity

While intestinal Th17 cells are critical for maintaining tissue homeostasis, recent studies have implicated their roles in the development of extra-intestinal autoimmune diseases including multiple sclerosis. However, the mechanisms by which tissue Th17 cells mediate these dichotomous functions remain unknown. Here, we characterized the heterogeneity, plasticity, and migratory phenotypes of tissue Th17 cells in vivo by combined fate mapping with profiling of the transcriptomes and TCR clonotypes of over 84,000 Th17 cells at homeostasis and during CNS autoimmune inflammation. Inter- and intra-organ single-cell analyses revealed a homeostatic, stem-like TCF1+ IL-17+ SLAMF6+ population that traffics to the intestine where it is maintained by the microbiota, providing a ready reservoir for the IL-23-driven generation of encephalitogenic GM-CSF+ IFN-γ+ CXCR6+ T cells. Our study defines a direct in vivo relationship between IL-17+ non-pathogenic and GM-CSF+ and IFN-γ+ pathogenic Th17 populations and provides a mechanism by which homeostatic intestinal Th17 cells direct extra-intestinal autoimmune disease.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure^1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63⁺ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

Single-cell clonotype tracing of tissue Th17 cells during CNS autoimmunity reveals the conversion of intestinal homeostatic cells to pathogenic cells

At homeostasis, most Th17 cells are found in the lamina propria of the intestine, where they contribute to tissue homeostasis by inhibiting microbiota from tissue invasion and promote barrier functions. Recent studies in humans and mice have implicated intestinal Th17 cells in extra-intestinal autoimmune diseases, including multiple sclerosis, but the mechanism in which intestinal Th17 cells mediate these dichotomous functions remains unknown. Here, we combined single-cell RNA- and TCR-seq (scRNA/TCR-seq) with fate mapping of tissue Th17 cells to profile over 84,000 tissue Th17 cells and characterize their heterogeneity, plasticity, and migration at homeostasis and during CNS autoimmunity. We discovered a homeostatic Th17 cell population, that is induced by the intestinal microbiota, is present in both lymphoid organs and the intestine, and expresses IL-17. Joint scRNA- and TCR-seq showed that during EAE this homeostatic population gives rise to a pathogenic Th17 cell population, that migrates specifically through the draining lymph nodes and the spleen to the CNS, and highly expresses GM-CSF and IFNγ, but expresses little or low levels of IL-17. This conversion of homeostatic Th17 cells into encephalitogenic Th17 cells depended on IL-23R signaling. Our work characterizes Th17 heterogeneity, plasticity, and migration during homeostasis and CNS autoimmune disease, identifies a mechanism by which intestinal Th17 cells influence the generation of an autoimmune disease-inducing population, and provides further evidence for a direct relationship between the IL-17-producing and GM-CSF- and IFNγ- producing pathogenic T cells.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries

Genome targeting methods enable cost-effective capture of specific subsets of the genome for sequencing. We present here an automated, highly scalable method for carrying out the Solution Hybrid Selection capture approach that provides a dramatic increase in scale and throughput of sequence-ready libraries produced. Significant process improvements and a series of in-process quality control checkpoints are also added. These process improvements can also be used in a manual version of the protocol.

Pharmacological subtypes of the gamma-aminobutyric acidA receptors defined by a gamma-aminobutyric acid analogue 4,5,6,7-tetrahydroisoxazolo[5,4-c] pyridin-3-ol and allosteric coupling: characterization using subunit-specific antibodies

Various alpha and beta 3 subunit-specific antibodies were used to characterize some of the heterogeneous ligand-binding properties of gamma-aminobutyric acidA receptors. Polyclonal antibodies that were raised against the cytoplasmic amino acid sequence (380-392) of the rat beta 3 subunit recognized a single polypeptide of molecular mass of 58 kDa in Western blots with Ro7-1986 affinity-purified GABAA receptors from the rat brain, and a doublet of molecular mass of 54 kDa and 56 kDa in receptors from the bovine cortex, hippocampus, and cerebellum. Deglycosylation of purified receptors from the bovine cortex with N-glycanase resulted in a single band immunostained at molecular mass of 52 kDa. These anti-beta 3 subunit antibodies immunoprecipitated approximately 50% of [3H]flunitrazepam binding sites from soluble extracts of bovine cortex, whereas beta cyto antibodies, which probably recognize all beta subunit isoforms, precipitated almost 100% of benzodiazepine binding sites. These results indicate heterogeneity of GABAA receptor subunit composition with respect to the nature of beta subunits. The GABA analogue 4,5,6,7-tetrahydroisoxazolo[5,4-c]pyridin-3-ol (THIP), like GABA, shows heterogeneous binding affinities in brain homogenates. The higher affinity sites were previously suggested as corresponding to a 58-kDa polypeptide in rat that is photoaffinity-labeled with [3H]muscimol, a band that comigrates with the one stained by anti-beta 3 antibodies. However, THIP affinity was not significantly different between receptors containing beta 3 subunits and those lacking beta 3, as demonstrated by similar affinities in receptors that ere immunoprecipitated by anti-beta 3 antibodies and those that were not. Also, THIP displaced [3H]muscimol binding with similar multiple affinities across brain regions where different beta subunit variants are expressed with varying abundances. These observations suggest that the property of high affinity THIP binding cannot be explained solely by beta 3 subunits. The coupling efficiency between GABA and benzodiazepine binding sites appears to be determined by the nature of alpha subunits rather than of beta subunits. GABA enhanced [3H]flunitrazepam binding with different efficacies and potencies in receptors immunoprecipitated by anti-alpha 1, -alpha 2, and -alpha 3 subunit antibodies. In contrast, beta 3 subunit-enriched and disenriched receptors did not differ in this property. [3H]Flunitrazepam binding in GABAA receptors containing alpha 2 and alpha 3 subunits was enhanced to a significantly greater extent than were those with alpha 1. In addition, receptors containing alpha 1 and alpha 3 subunits had higher potencies of enhancement than did those with alpha 2 subunits.(ABSTRACT TRUNCATED AT 400 WORDS)