Defining the T cell transcriptional landscape in pediatric liver transplant rejection at single cell resolution.

Acute cellular rejection (ACR) affects >80% of pediatric liver transplant recipients within 5 years, and late ACR is associated with graft failure. Traditional anti-rejection therapy for late ACR is ineffective and has remained unchanged for six decades. Although CD8+ T cells promote late ACR, little has been done to define their specificity and gene expression. Here, we used single-cell sequencing and immune repertoire profiling (10X Genomics) on 30 cryopreserved 16G liver biopsies from 14 patients (5 pre-transplant or with no ACR, 9 with ACR). We identified expanded intragraft CD8+ T cell clonotypes (CD8EXP) and their gene expression profiles in response to anti-rejection treatment. Notably, we found that expanded CD8+ clonotypes (CD8EXP) bore markers of effector and CD56hiCD161- 'NK-like' T cells, retaining their clonotype identity and phenotype in subsequent biopsies from the same patients despite histologic ACR resolution. CD8EXP clonotypes localized to portal infiltrates during active ACR, and persisted in the lobule after histologic ACR resolution. CellPhoneDB analysis revealed differential crosstalk between KC and CD8EXP during late ACR, with activation of the LTB-LTBR pathway and downregulation of TGFß signaling. Therefore, persistently-detected intragraft CD8EXP clones remain active despite ACR treatment and may contribute to long-term allograft fibrosis and failure of operational tolerance.

Disease diagnostics using machine learning of immune receptors.

Clinical diagnosis typically incorporates physical examination, patient history, and various laboratory tests and imaging studies, but makes limited use of the human system's own record of antigen exposures encoded by receptors on B cells and T cells. We analyzed immune receptor datasets from 593 individuals to develop MAchine Learning for Immunological Diagnosis (Mal-ID) , an interpretive framework to screen for multiple illnesses simultaneously or precisely test for one condition. This approach detects specific infections, autoimmune disorders, vaccine responses, and disease severity differences. Human-interpretable features of the model recapitulate known immune responses to SARS-CoV-2, Influenza, and HIV, highlight antigen-specific receptors, and reveal distinct characteristics of Systemic Lupus Erythematosus and Type-1 Diabetes autoreactivity. This analysis framework has broad potential for scientific and clinical interpretation of human immune responses.

VExD: a curated resource for human gene expression alterations following viral infection.

Much of the host antiviral response is mediated through changes to host gene expression levels. Likewise, viruses induce changes to host gene expression levels in order to promote the viral life cycle and evade the host immune system. However, there is no resource that specifically collects human gene expression levels pre- and post-virus infection. Further, public gene expression repositories do not contain enough specialized metadata to easily find relevant experiments. Here, we present the Virus Expression Database (VExD), a freely available website and database, that collects human gene expression datasets in response to viral infection. VExD contains ∼8,000 uniformly processed samples obtained from 289 studies examining 51 distinct human viruses. We show that the VExD processing pipeline captures known antiviral responses in the form of interferon-stimulated genes. We further show that the datasets collected in VExD can be used to quickly identify supporting data for experiments performed in human cells or model organisms. VExD is freely available at https://vexd.cchmc.org/.

Single-cell transcriptomic analysis of renal allograft rejection reveals insights into intragraft TCR clonality.

Bulk analysis of renal allograft biopsies (rBx) identified RNA transcripts associated with acute cellular rejection (ACR); however, these lacked cellular context critical to mechanistic understanding of how rejection occurs despite immunosuppression (IS). We performed combined single-cell RNA transcriptomic and TCR-α/ß sequencing on rBx from patients with ACR under differing IS drugs: tacrolimus, iscalimab, and belatacept. We found distinct CD8+ T cell phenotypes (e.g., effector, memory, exhausted) depending upon IS type, particularly within expanded CD8+ T cell clonotypes (CD8EXP). Gene expression of CD8EXP identified therapeutic targets that were influenced by IS type. TCR analysis revealed a highly restricted number of CD8EXP, independent of HLA mismatch or IS type. Subcloning of TCR-α/ß cDNAs from CD8EXP into Jurkat 76 cells (TCR-/-) conferred alloreactivity by mixed lymphocyte reaction. Analysis of sequential rBx samples revealed persistence of CD8EXP that decreased, but were not eliminated, after successful antirejection therapy. In contrast, CD8EXP were maintained in treatment-refractory rejection. Finally, most rBx-derived CD8EXP were also observed in matching urine samples, providing precedent for using urine-derived CD8EXP as a surrogate for those found in the rejecting allograft. Overall, our data define the clonal CD8+ T cell response to ACR, paving the next steps for improving detection, assessment, and treatment of rejection.

Effects of invivo CXCR4 blockade and proteasome inhibition on bone marrow plasma cells in HLA-sensitized kidney transplant candidates.

To date, plasma cell (PC)-targeted therapies have been limited by suboptimal PC depletion and antibody rebound. We hypothesized this is partly because of PC residence in protective bone marrow (BM) microenvironments. The purpose of this proof-of-concept study was to examine the effects of the CXCR4 antagonist, plerixafor, on PC BM residence; its safety profile (alone and in combination with a proteasome inhibitor, bortezomib); and the transcriptional effect on BMPCs in HLA-sensitized kidney transplant candidates. Participants were enrolled into 3 groups: group A (n = 4), plerixafor monotherapy; and groups B (n = 4) and C (n = 4), plerixafor and bortezomib combinations. CD34+ stem cell and PC levels increased in the blood after plerixafor treatment. PC recovery from BM aspirates varied depending on the dose of plerixafor and bortezomib. Single-cell RNA sequencing on BMPCs from 3 group C participants pretreatment and posttreatment revealed multiple populations of PCs, with a posttreatment enrichment of oxidative phosphorylation, proteasome assembly, cytoplasmic translation, and autophagy-related genes. Murine studies demonstrated dually inhibiting the proteasome and autophagy resulted in greater BMPC death than did monotherapies. In conclusion, this pilot study revealed anticipated effects of combined plerixafor and bortezomib on BMPCs, an acceptable safety profile, and suggests the potential for autophagy inhibitors in desensitization regimens.

B cell repertoire in children with skin barrier dysfunction supports altered IgE maturation associated with allergic food sensitization.

The skin is a major immune organ and skin barrier dysfunction is a major risk factor for the development of the inappropriate immune response seen in allergic disease. Skin barrier disruption alters the landscape of antigens experienced by the immune system and the downstream impacts on the antibody repertoire remain poorly characterized, particularly for the IgE isotype responsible for allergic specificity and in early life, when allergic disease is developing. In this study, we sequenced antibody gene repertoires from a large and well-characterized cohort of children with atopic dermatitis and found that food sensitization was associated with lower mutation frequencies in the IgE compartment. This trend was abrogated in children living with pets during the first year of life. These results elucidate potential molecular mechanisms underlying the protective effects of pet ownership and non-antiseptic environs reported for allergic disease, and the hygiene hypothesis more broadly. We also observed increased IgE diversity and increased isotype-switching to the IgE isotype, suggesting that B cell development, particularly isotype-switching, is heavily altered in the those with food allergen sensitizations relative to those without food allergen sensitizations. Unlike for food antigens, aeroallergen sensitization exhibited no effect on IgE mutation or diversity. Consistent patterns of antibody rearrangement were associated with food allergen sensitization in subjects with atopic dermatitis. Thus, we propose the Immune Repertoire in Atopic Disease (IRAD) score, to quantify this repertoire shift and to aid clinically in patient diagnosis and risk stratification.

Multi-omic characterization of pediatric ARDS via nasal brushings.

RATIONALE: While nasal brushing transcriptomics can identify disease subtypes in chronic pulmonary diseases, it is unknown whether this is true in pediatric acute respiratory distress syndrome (PARDS). OBJECTIVES: Determine whether nasal transcriptomics and methylomics can identify clinically meaningful PARDS subgroups that reflect important pathobiological processes. METHODS: Nasal brushings and serum were collected on days 1, 3, 7, and 14 from control and PARDS subjects from two centers. PARDS duration was the primary endpoint. MEASUREMENTS AND MAIN RESULTS: Twenty-four control and 39 PARDS subjects were enrolled. Two nasal methylation patterns were identified. Compared to Methyl Subgroup 1, Subgroup 2 had hypomethylation of inflammatory genes and was enriched for immunocompromised subjects. Four transcriptomic patterns were identified with temporal patterns indicating injury, repair, and regeneration. Over time, both inflammatory (Subgroup B) and cell injury (Subgroup D) patterns transitioned to repair (Subgroup A) and eventually homeostasis (Subgroup C). When control specimens were included, they were largely Subgroup C. In comparison with 17 serum biomarkers, the nasal transcriptome was more predictive of prolonged PARDS. Subjects with initial Transcriptomic Subgroup B or D assignment had median PARDS duration of 8 days compared to 2 in A or C (p = 0.02). For predicting PARDS duration ≥ 3 days, nasal transcriptomics was more sensitive and serum biomarkers more specific. CONCLUSIONS: PARDS nasal transcriptome may reflect distal lung injury, repair, and regeneration. A combined nasal PCR and serum biomarker assay could be useful for predictive and diagnostic enrichment. Trial registration Clinicaltrials.gov NCT03539783 May 29, 2018.

Lymphoid blast transformation in an MPN with BCR-JAK2 treated with ruxolitinib: putative mechanisms of resistance.

The basis for acquired resistance to JAK inhibition in patients with JAK2-driven hematologic malignancies is not well understood. We report a patient with a myeloproliferative neoplasm (MPN) with a BCR activator of RhoGEF and GTPase (BCR)-JAK2 fusion with initial hematologic response to ruxolitinib who rapidly developed B-lymphoid blast transformation. We analyzed pre-ruxolitinib and blast transformation samples using genome sequencing, DNA mate-pair sequencing (MPseq), RNA sequencing (RNA-seq), and chromosomal microarray to characterize possible mechanisms of resistance. No resistance mutations in the BCR-JAK2 fusion gene or transcript were identified, and fusion transcript expression levels remained stable. However, at the time of blast transformation, MPseq detected a new IKZF1 copy-number loss, which is predicted to result in loss of normal IKZF1 protein translation. RNA-seq revealed significant upregulation of genes negatively regulated by IKZF1, including IL7R and CRLF2. Disease progression was also characterized by adaptation to an activated B-cell receptor (BCR)-like signaling phenotype, with marked upregulation of genes such as CD79A, CD79B, IGLL1, VPREB1, BLNK, ZAP70, RAG1, and RAG2. In summary, IKZF1 deletion and a switch from cytokine dependence to activated BCR-like signaling phenotype represent putative mechanisms of ruxolitinib resistance in this case, recapitulating preclinical data on resistance to JAK inhibition in CRLF2-rearranged Philadelphia chromosome-like acute lymphoblastic leukemia.

Shared B cell memory to coronaviruses and other pathogens varies in human age groups and tissues.

Vaccination and infection promote the formation, tissue distribution, and clonal evolution of B cells, which encode humoral immune memory. We evaluated pediatric and adult blood and deceased adult organ donor tissues to determine convergent antigen-specific antibody genes of similar sequences shared between individuals. B cell memory varied for different pathogens. Polysaccharide antigen-specific clones were not exclusive to the spleen. Adults had higher clone frequencies and greater class switching in lymphoid tissues than blood, while pediatric blood had abundant class-switched convergent clones. Consistent with reported serology, prepandemic children had class-switched convergent clones to severe acute respiratory syndrome coronavirus 2 with weak cross-reactivity to other coronaviruses, while adult blood or tissues showed few such clones. These results highlight the prominence of early childhood B cell clonal expansions and cross-reactivity for future responses to novel pathogens.

Modeling human adaptive immune responses with tonsil organoids.

Most of what we know about adaptive immunity has come from inbred mouse studies, using methods that are often difficult or impossible to confirm in humans. In addition, vaccine responses in mice are often poorly predictive of responses to those same vaccines in humans. Here we use human tonsils, readily available lymphoid organs, to develop a functional organotypic system that recapitulates key germinal center features in vitro, including the production of antigen-specific antibodies, somatic hypermutation and affinity maturation, plasmablast differentiation and class-switch recombination. We use this system to define the essential cellular components necessary to produce an influenza vaccine response. We also show that it can be used to evaluate humoral immune responses to two priming antigens, rabies vaccine and an adenovirus-based severe acute respiratory syndrome coronavirus 2 vaccine, and to assess the effects of different adjuvants. This system should prove useful for studying critical mechanisms underlying adaptive immunity in much greater depth than previously possible and to rapidly test vaccine candidates and adjuvants in an entirely human system.

Origins and clonal convergence of gastrointestinal IgE+ B cells in human peanut allergy.

B cells in human food allergy have been studied predominantly in the blood. Little is known about IgE+ B cells or plasma cells in tissues exposed to dietary antigens. We characterized IgE+ clones in blood, stomach, duodenum, and esophagus of 19 peanut-allergic patients, using high-throughput DNA sequencing. IgE+ cells in allergic patients are enriched in stomach and duodenum, and have a plasma cell phenotype. Clonally related IgE+ and non-IgE-expressing cell frequencies in tissues suggest local isotype switching, including transitions between IgA and IgE isotypes. Highly similar antibody sequences specific for peanut allergen Ara h 2 are shared between patients, indicating that common immunoglobulin genetic rearrangements may contribute to pathogenesis. These data define the gastrointestinal tract as a reservoir of IgE+ B lineage cells in food allergy.

Aberrant B cell repertoire selection associated with HIV neutralizing antibody breadth.

A goal of HIV vaccine development is to elicit antibodies with neutralizing breadth. Broadly neutralizing antibodies (bNAbs) to HIV often have unusual sequences with long heavy-chain complementarity-determining region loops, high somatic mutation rates and polyreactivity. A subset of HIV-infected individuals develops such antibodies, but it is unclear whether this reflects systematic differences in their antibody repertoires or is a consequence of rare stochastic events involving individual clones. We sequenced antibody heavy-chain repertoires in a large cohort of HIV-infected individuals with bNAb responses or no neutralization breadth and uninfected controls, identifying consistent features of bNAb repertoires, encompassing thousands of B cell clones per individual, with correlated T cell phenotypes. These repertoire features were not observed during chronic cytomegalovirus infection in an independent cohort. Our data indicate that the development of numerous B cell lineages with antibody features associated with autoreactivity may be a key aspect in the development of HIV neutralizing antibody breadth.

B cell clonal expansion within immune infiltrates in human cardiac allograft vasculopathy.

Cardiac allograft vasculopathy (CAV) is associated with intragraft B cell infiltrates. Here, we studied the clonal composition of B cell infiltrates using 4 graft specimens with CAV. Using deep sequencing, we analyzed the immunoglobulin heavy chain variable region repertoire in both graft and blood. Results showed robust B cell clonal expansion in the graft but not in the blood for all cases. Several expanded B cell clones, characterized by their uniquely rearranged complementarity-determining region 3, were detected in different locations in the graft. Sequences from intragraft B cells also showed elevated levels of mutated rearrangements in the graft compared to blood B cells. The number of somatic mutations per rearrangement was also higher in the graft than in the blood, suggesting that B cells continued maturing in situ. Overall, our studies demonstrated B cell clonal expansion in human cardiac allografts with CAV. This local B cell response may contribute to the pathophysiology of CAV through a mechanism that needs to be identified.

Shaping of infant B cell receptor repertoires by environmental factors and infectious disease.

Antigenic exposures at epithelial sites in infancy and early childhood are thought to influence the maturation of humoral immunity and modulate the risk of developing immunoglobulin E (IgE)-mediated allergic disease. How different kinds of environmental exposures influence B cell isotype switching to IgE, IgG, or IgA, and the somatic mutation maturation of these antibody pools, is not fully understood. We sequenced antibody repertoires in longitudinal blood samples in a birth cohort from infancy through the first 3 years of life and found that, whereas IgG and IgA show linear increases in mutational maturation with age, IgM and IgD mutations are more closely tied to pathogen exposure. IgE mutation frequencies are primarily increased in children with impaired skin barrier conditions such as eczema, suggesting that IgE affinity maturation could provide a mechanistic link between epithelial barrier failure and allergy development.

Impact of somatic and germline mutations on the outcome of systemic mastocytosis.

Systemic mastocytosis (SM) is a highly heterogeneous disease with indolent and aggressive forms, with the mechanisms leading to malignant transformation still remaining to be elucidated. Here, we investigated the presence and frequency of genetic variants in 34 SM patients with multilineal KIT D816V mutations. Initial screening was performed by targeted sequencing of 410 genes in DNA extracted from purified bone marrow cells and hair from 12 patients with nonadvanced SM and 8 patients with advanced SM, followed by whole-genome sequencing (WGS) in 4 cases. Somatic mutations were further investigated in another 14 patients with advanced SM. Despite the fact that no common mutation other than KIT D816V was found in WGS analyses, targeted next-generation sequencing identified 67 nonsynonymous genetic variants involving 39 genes. Half of the mutations were somatic (mostly multilineal), whereas the other half were germline variants. The presence of ≥1 multilineal somatic mutation involving genes other than KIT D816V, ≥3 germline variants, and ≥1 multilineal mutation in the SRSF2, ASXL1, RUNX1, and/or EZH2 genes (S/A/R/E genes), in addition to skin lesions, splenomegaly, thrombocytopenia, low hemoglobin levels, and increased alkaline phosphatase and ß2-microglobulin serum levels, were associated with a poorer patient outcome. However, the presence of ≥1 multilineal mutation, particularly involving S/A/R/E genes, was the only independent predictor for progression-free survival and overall survival in our cohort.

A novel TRIP11-FLT3 fusion in a patient with a myeloid/lymphoid neoplasm with eosinophilia.

FLT3 fusions are associated with myeloid and lymphoid neoplasms with eosinophilia. We describe a patient presenting with clinicopathologic features of both chronic eosinophilic leukemia, not otherwise specified (CEL, NOS) and systemic mastocytosis (SM). The bone marrow demonstrated a myeloproliferative neoplasm with eosinophilia and aggregates of atypical mast cells. Cytogenetic analysis revealed a t(13;14)(q12;q32), which was subsequently molecularly characterized as a novel TRIP11-FLT3 rearrangement. A KIT D816V mutation was also identified. The patient rapidly transformed to T-lymphoblastic leukemia/lymphoma and expired shortly after diagnosis. This is the fifth FLT3 fusion gene described in the literature; the presence of both myeloid and lymphoid neoplasms implicates involvement of an early hematopoietic progenitor by rearranged FLT3. We suggest that leukemias and lymphomas with FLT3 fusion genes exhibit similar clinicopathologic features to, and should be included in, the WHO category of "Myeloid and lymphoid neoplasms with eosinophilia and abnormalities of PDGFRA, PDGFRB, or FGFR1, or with PCM1-JAK2."

Defining antigen-specific plasmablast and memory B cell subsets in human blood after viral infection or vaccination.

Antigen-specific B cells bifurcate into antibody-secreting cells (ASCs) and memory B cells (MBCs) after infection or vaccination. ASCs (plasmablasts) have been extensively studied in humans, but less is known about B cells that become activated but do not differentiate into plasmablasts. Here we have defined the phenotype and transcriptional program of a subset of antigen-specific B cells, which we have called 'activated B cells' (ABCs), that were distinct from ASCs and were committed to the MBC lineage. We detected ABCs in humans after infection with Ebola virus or influenza virus and also after vaccination. By simultaneously analyzing antigen-specific ASCs and ABCs in human blood after vaccination against influenza virus, we investigated the clonal overlap and extent of somatic hypermutation (SHM) in the ASC (effector) and ABC (memory) lineages. Longitudinal tracking of vaccination-induced hemagglutinin (HA)-specific clones revealed no overall increase in SHM over time, which suggested that repeated annual immunization might have limitations in enhancing the quality of influenza-virus-specific antibody.

Human B-cell isotype switching origins of IgE.

BACKGROUND: B cells expressing IgE contribute to immunity against parasites and venoms and are the source of antigen specificity in allergic patients, yet the developmental pathways producing these B cells in human subjects remain a subject of debate. Much of our knowledge of IgE lineage development derives from model studies in mice rather than from human subjects. OBJECTIVE: We evaluate models for isotype switching to IgE in human subjects using immunoglobulin heavy chain (IGH) mutational lineage data. METHODS: We analyzed IGH repertoires in 9 allergic and 24 healthy adults using high-throughput DNA sequencing of 15,843,270 IGH rearrangements to identify clonal lineages of B cells containing members expressing IgE. Somatic mutations in IGH inherited from common ancestors within the clonal lineage are used to infer the relationships between B cells. RESULTS: Data from 613,641 multi-isotype B-cell clonal lineages, of which 592 include an IgE member, are consistent with indirect switching to IgE from IgG- or IgA-expressing lineage members in human subjects. We also find that these inferred isotype switching frequencies are similar in healthy and allergic subjects. CONCLUSIONS: We found evidence that secondary isotype switching of mutated IgG1-expressing B cells is the primary source of IgE in human subjects, with lesser contributions from precursors expressing other switched isotypes and rarely IgM or IgD, suggesting that IgE is derived from previously antigen-experienced B cells rather than naive B cells that typically express low-affinity unmutated antibodies. These data provide a basis from which to evaluate allergen-specific human antibody repertoires in healthy and diseased subjects.

Single B-cell deconvolution of peanut-specific antibody responses in allergic patients.

BACKGROUND: The frequencies, cellular phenotypes, epitope specificity, and clonal diversity of allergen-specific B cells in patients with food allergy are not fully understood but are of major pathogenic and therapeutic significance. OBJECTIVE: We sought to characterize peanut allergen-specific B-cell populations and the sequences and binding activities of their antibodies before and during immunotherapy. METHODS: B cells binding fluorescently labeled Ara h 1 or Ara h 2 were phenotyped and isolated by means of flow cytometric sorting from 18 patients at baseline and 13 patients during therapy. Fifty-seven mAbs derived from allergen-binding single B cells were evaluated by using ELISA, Western blotting, and peptide epitope mapping. Deep sequencing of the B-cell repertoires identified additional members of the allergen-specific B-cell clones. RESULTS: Median allergen-binding B-cell frequencies were 0.0097% (Ara h 1) or 0.029% (Ara h 2) of B cells in baseline blood from allergic patients and approximately 3-fold higher during immunotherapy. Five of 57 allergen-specific cells belonged to clones containing IgE-expressing members. Almost all allergen-specific antibodies were mutated, and binding to both conformational and linear allergen epitopes was detected. Increasing somatic mutation of IgG4 members of a clone was seen in immunotherapy, whereas IgE mutation levels in the clone did not increase. CONCLUSION: Most peanut allergen-binding B cells isolated by means of antigen-specific flow sorting express mutated and isotype-switched antibodies. Immunotherapy increases their frequency in the blood, and even narrowly defined allergen epitopes are recognized by numerous distinct B-cell clones in a patient. The results also suggest that oral immunotherapy can stimulate somatic mutation of allergen-specific IgG4.

Laboratory and Data Analysis Methods for Characterization of Human B Cell Repertoires by High-Throughput DNA Sequencing.

High-throughput DNA sequencing techniques have greatly accelerated the pace of research into the repertoires of antibody and T cell receptor gene rearrangements that confer antigen specificity to adaptive immune responses. Studies of aging-related changes in human B cell repertoires have benefited from the ability to detect and quantify thousands to millions of B cell clones in human samples, and study the mutational lineages and isotype switching relationships within each clonal lineage. Correlation of repertoire analysis with antibody gene data from antigen-specific B cells is poised to give much greater insight into clinically relevant B cell responses and memory storage. Here, we describe strategies for preparing and analyzing human antibody gene libraries for studying B cell repertoires.