Induction of bronchus-associated lymphoid tissue is an early life adaptation for promoting human B cell immunity.

Infants and young children are more susceptible to common respiratory pathogens than adults but can fare better against novel pathogens like severe acute respiratory syndrome coronavirus 2. The mechanisms by which infants and young children mount effective immune responses to respiratory pathogens are unknown. Through investigation of lungs and lung-associated lymph nodes from infant and pediatric organ donors aged 0-13 years, we show that bronchus-associated lymphoid tissue (BALT), containing B cell follicles, CD4+ T cells and functionally active germinal centers, develop during infancy. BALT structures are prevalent around lung airways during the first 3 years of life, and their numbers decline through childhood coincident with the accumulation of memory T cells. Single-cell profiling and repertoire analysis reveals that early life lung B cells undergo differentiation, somatic hypermutation and immunoglobulin class switching and exhibit a more activated profile than lymph node B cells. Moreover, B cells in the lung and lung-associated lymph nodes generate biased antibody responses to multiple respiratory pathogens compared to circulating antibodies, which are mostly specific for vaccine antigens in the early years of life. Together, our findings provide evidence for BALT as an early life adaptation for mobilizing localized immune protection to the diverse respiratory challenges during this formative life stage.

Targeting BMI-1 in B cells restores effective humoral immune responses and controls chronic viral infection.

Ineffective antibody-mediated responses are a key characteristic of chronic viral infection. However, our understanding of the intrinsic mechanisms that drive this dysregulation are unclear. Here, we identify that targeting the epigenetic modifier BMI-1 in mice improves humoral responses to chronic lymphocytic choriomeningitis virus. BMI-1 was upregulated by germinal center B cells in chronic viral infection, correlating with changes to the accessible chromatin landscape, compared to acute infection. B cell-intrinsic deletion of Bmi1 accelerated viral clearance, reduced splenomegaly and restored splenic architecture. Deletion of Bmi1 restored c-Myc expression in B cells, concomitant with improved quality of antibody and coupled with reduced antibody-secreting cell numbers. Specifically, BMI-1-deficiency induced antibody with increased neutralizing capacity and enhanced antibody-dependent effector function. Using a small molecule inhibitor to murine BMI-1, we could deplete antibody-secreting cells and prohibit detrimental immune complex formation in vivo. This study defines BMI-1 as a crucial immune modifier that controls antibody-mediated responses in chronic infection.

Liver Is a Generative Site for the B Cell Response to Ehrlichia muris.

The B cell response to Ehrlichia muris is dominated by plasmablasts (PBs), with few-if any-germinal centers (GCs), yet it generates protective immunoglobulin M (IgM) memory B cells (MBCs) that express the transcription factor T-bet and harbor V-region mutations. Because Ehrlichia prominently infects the liver, we investigated the nature of liver B cell response and that of the spleen. B cells within infected livers proliferated and underwent somatic hypermutation (SHM). Vh-region sequencing revealed trafficking of clones between the spleen and liver and often subsequent local clonal expansion and intraparenchymal localization of T-bet+ MBCs. T-bet+ MBCs expressed MBC subset markers CD80 and PD-L2. Many T-bet+ MBCs lacked CD11b or CD11c expression but had marginal zone (MZ) B cell phenotypes and colonized the splenic MZ, revealing T-bet+ MBC plasticity. Hence, liver and spleen are generative sites of B cell responses, and they include V-region mutation and result in liver MBC localization.

Applying systems-level spectral imaging and analysis to reveal the organelle interactome.

The organization of the eukaryotic cell into discrete membrane-bound organelles allows for the separation of incompatible biochemical processes, but the activities of these organelles must be coordinated. For example, lipid metabolism is distributed between the endoplasmic reticulum for lipid synthesis, lipid droplets for storage and transport, mitochondria and peroxisomes for ß-oxidation, and lysosomes for lipid hydrolysis and recycling. It is increasingly recognized that organelle contacts have a vital role in diverse cellular functions. However, the spatial and temporal organization of organelles within the cell remains poorly characterized, as fluorescence imaging approaches are limited in the number of different labels that can be distinguished in a single image. Here we present a systems-level analysis of the organelle interactome using a multispectral image acquisition method that overcomes the challenge of spectral overlap in the fluorescent protein palette. We used confocal and lattice light sheet instrumentation and an imaging informatics pipeline of five steps to achieve mapping of organelle numbers, volumes, speeds, positions and dynamic inter-organelle contacts in live cells from a monkey fibroblast cell line. We describe the frequency and locality of two-, three-, four- and five-way interactions among six different membrane-bound organelles (endoplasmic reticulum, Golgi, lysosome, peroxisome, mitochondria and lipid droplet) and show how these relationships change over time. We demonstrate that each organelle has a characteristic distribution and dispersion pattern in three-dimensional space and that there is a reproducible pattern of contacts among the six organelles, that is affected by microtubule and cell nutrient status. These live-cell confocal and lattice light sheet spectral imaging approaches are applicable to any cell system expressing multiple fluorescent probes, whether in normal conditions or when cells are exposed to disturbances such as drugs, pathogens or stress. This methodology thus offers a powerful descriptive tool and can be used to develop hypotheses about cellular organization and dynamics.

Caspase-8 promotes c-Rel-dependent inflammatory cytokine expression and resistance against Toxoplasma gondii.

Caspase-8 is a key integrator of cell survival and cell death decisions during infection and inflammation. Following engagement of tumor necrosis factor superfamily receptors or certain Toll-like receptors (TLRs), caspase-8 initiates cell-extrinsic apoptosis while inhibiting RIPK3-dependent programmed necrosis. In addition, caspase-8 has an important, albeit less well understood, role in cell-intrinsic inflammatory gene expression. Macrophages lacking caspase-8 or the adaptor FADD have defective inflammatory cytokine expression and inflammasome priming in response to bacterial infection or TLR stimulation. How caspase-8 regulates cytokine gene expression, and whether caspase-8-mediated gene regulation has a physiological role during infection, remain poorly defined. Here we demonstrate that both caspase-8 enzymatic activity and scaffolding functions contribute to inflammatory cytokine gene expression. Caspase-8 enzymatic activity was necessary for maximal expression of Il1b and Il12b, but caspase-8 deficient cells exhibited a further decrease in expression of these genes. Furthermore, the ability of TLR stimuli to induce optimal IκB kinase phosphorylation and nuclear translocation of the nuclear factor kappa light chain enhancer of activated B cells family member c-Rel required caspase activity. Interestingly, overexpression of c-Rel was sufficient to restore expression of IL-12 and IL-1ß in caspase-8-deficient cells. Moreover, Ripk3-/-Casp8-/- mice were unable to control infection by the intracellular parasite Toxoplasma gondii, which corresponded to defects in monocyte recruitment to the peritoneal cavity, and exogenous IL-12 restored monocyte recruitment and protection of caspase-8-deficient mice during acute toxoplasmosis. These findings provide insight into how caspase-8 controls inflammatory gene expression and identify a critical role for caspase-8 in host defense against eukaryotic pathogens.

Human Lymph Nodes Maintain TCF-1hi Memory T Cells with High Functional Potential and Clonal Diversity throughout Life.

Translating studies on T cell function and modulation from mouse models to humans requires extrapolating in vivo results on mouse T cell responses in lymphoid organs (spleen and lymph nodes [LN]) to human peripheral blood T cells. However, our understanding of T cell responses in human lymphoid sites and their relation to peripheral blood remains sparse. In this study, we used a unique human tissue resource to study human T cells in different anatomical compartments within individual donors and identify a subset of memory CD8+ T cells in LN, which maintain a distinct differentiation and functional profile compared with memory CD8+ T cells in blood, spleen, bone marrow, and lungs. Whole-transcriptome and high-dimensional cytometry by time-of-flight profiling reveals that LN memory CD8+ T cells express signatures of quiescence and self-renewal compared with corresponding populations in blood, spleen, bone marrow, and lung. LN memory T cells exhibit a distinct transcriptional signature, including expression of stem cell-associated transcription factors TCF-1 and LEF-1, T follicular helper cell markers CXCR5 and CXCR4, and reduced expression of effector molecules. LN memory T cells display high homology to a subset of mouse CD8+ T cells identified in chronic infection models that respond to checkpoint blockade immunotherapy. Functionally, human LN memory T cells exhibit increased proliferation to TCR-mediated stimulation and maintain higher TCR clonal diversity compared with memory T cells from blood and other sites. These findings establish human LN as reservoirs for memory T cells with high capacities for expansion and diverse recognition and important targets for immunotherapies.

ImmuneDB: a system for the analysis and exploration of high-throughput adaptive immune receptor sequencing data.

As high-throughput sequencing of B cells becomes more common, the need for tools to analyze the large quantity of data also increases. This article introduces ImmuneDB, a system for analyzing vast amounts of heavy chain variable region sequences and exploring the resulting data. It can take as input raw FASTA/FASTQ data, identify genes, determine clones, construct lineages, as well as provide information such as selection pressure and mutation analysis. It uses an industry leading database, MySQL, to provide fast analysis and avoid the complexities of using error prone flat-files. AVAILABILITY AND IMPLEMENTATION: ImmuneDB is freely available at http://immunedb.comA demo of the ImmuneDB web interface is available at: http://immunedb.com/demo CONTACT: Uh25@drexel.eduSupplementary information: Supplementary data are available at Bioinformatics online.

T Cell Receptor-induced Nuclear Factor κB (NF-κB) Signaling and Transcriptional Activation Are Regulated by STIM1- and Orai1-mediated Calcium Entry.

T cell activation following antigen binding to the T cell receptor (TCR) involves the mobilization of intracellular Ca(2+) to activate the key transcription factors nuclear factor of activated T lymphocytes (NFAT) and NF-κB. The mechanism of NFAT activation by Ca(2+) has been determined. However, the role of Ca(2+) in controlling NF-κB signaling is poorly understood, and the source of Ca(2+) required for NF-κB activation is unknown. We demonstrate that TCR- but not TNF-induced NF-κB signaling upstream of IκB kinase activation absolutely requires the influx of extracellular Ca(2+) via STIM1-dependent Ca(2+) release-activated Ca(2+)/Orai channels. We further show that Ca(2+) influx controls phosphorylation of the NF-κB protein p65 on Ser-536 and that this posttranslational modification controls its nuclear localization and transcriptional activation. Notably, our data reveal that this role for Ca(2+) is entirely separate from its upstream control of IκBα degradation, thereby identifying a novel Ca(2+)-dependent distal step in TCR-induced NF-κB activation. Finally, we demonstrate that this control of distal signaling occurs via Ca(2+)-dependent PKCα-mediated phosphorylation of p65. Thus, we establish the source of Ca(2+) required for TCR-induced NF-κB activation and define a new distal Ca(2+)-dependent checkpoint in TCR-induced NF-κB signaling that has broad implications for the control of immune cell development and T cell functional specificity.

T cell-independent and toll-like receptor-dependent antigen-driven activation of autoreactive B cells.

On the lupus-prone MRL-lpr/lpr (MRL-lpr) background, AM14 rheumatoid factor (RF) B cells are activated, differentiate into plasmablasts, and undergo somatic hypermutation outside of follicles. Using multiple strategies to impair T cells, we found that such AM14 B cell activation did not require T cells but could be modulated by them. In vitro, the signaling adaptor MyD88 is required for IgG anti-chromatin to stimulate AM14 B cell proliferation when T cells are absent. However, the roles of Toll-like receptors (TLRs) in AM14 B cell activation in vivo have not been investigated. We found that activation, expansion, and differentiation of AM14 B cells depended on MyD88; however, mice lacking either TLR7 or TLR9 displayed partial defects, indicating complex roles for these receptors. T cell-independent activation of certain autoreactive B cells, which gain stimuli via endogenous TLR ligands instead of T cells, may be the initial step in the generation of canonical autoantibodies.

Ig gene analysis reveals altered selective pressures on Ig-producing cells in parotid glands of primary Sjögren's syndrome patients.

In this study, we sought to understand the selective pressures shaping the Ig-producing cell repertoire in the parotid glands of primary Sjögren's syndrome (pSS) patients before and after rituximab treatment (RTX). In particular, we evaluated the role of potential N-glycosylation motifs acquired by somatic hypermutation (ac-Nglycs) within Ig H chain V region (IGHV) genes as alternative selective pressures for B cells in pSS. Five pSS patients received RTX. Sequential parotid salivary gland biopsies were taken before RTX, at 12 wk and at 36-52 wk after treatment. Parotid biopsies from four non-pSS patients served as controls. Sequence analysis was carried out on the IgA and IgG RNA transcripts expressing IGHV3 genes in all parotid biopsies. Both IgG and IgA sequences from pSS patients exhibited no evidence for positive Ag-driven selection pressure in their CDRs in contrast to non-pSS controls. The prevalence of IgG sequences with ac-Nglycs was significantly higher in pSS patients than in non-pSS controls. Selection pressures shaping the IgG and IgA repertoire within pSS patients' parotid glands are distinct from those in non-pSS controls, with very little evidence for positive (auto)antigen selection. The higher prevalence of ac-Nglycs on pSS-IgG compared with non-pSS IgG indicates that ac-Nglycs could be an alternative form of selection pressure. We speculate that B cell hyperproliferation within parotid glands of pSS patients may result from Ag-independent interactions such as that between glycosylated B cell receptors and lectins within the microenvironment rather than (auto)antigen-specific stimulation. Our study brings a new perspective into research on pSS.

VDJML: a file format with tools for capturing the results of inferring immune receptor rearrangements.

BACKGROUND: The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses. RESULTS: To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format. CONCLUSIONS: The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/ . We welcome participation from the community in developing the file format standard, as well as code contributions.

Utilization of HIV-1 envelope V3 to identify X4- and R5-specific Tat and LTR sequence signatures.

BACKGROUND: HIV-1 entry is a receptor-mediated process directed by the interaction of the viral envelope with the host cell CD4 molecule and one of two co-receptors, CCR5 or CXCR4. The amino acid sequence of the third variable (V3) loop of the HIV-1 envelope is highly predictive of co-receptor utilization preference during entry, and machine learning predictive algorithms have been developed to characterize sequences as CCR5-utilizing (R5) or CXCR4-utilizing (X4). It was hypothesized that while the V3 loop is predominantly responsible for determining co-receptor binding, additional components of the HIV-1 genome may contribute to overall viral tropism and display sequence signatures associated with co-receptor utilization. RESULTS: The accessory protein Tat and the HlV-1 long terminal repeat (LTR) were analyzed with respect to genetic diversity and compared by Jensen-Shannon divergence which resulted in a correlation with both mean genetic diversity as well as the absolute difference in genetic diversity between R5- and X4-genome specific trends. As expected, the V3 domain of the gp120 protein was enriched with statistically divergent positions. Statistically divergent positions were also identified in Tat amino acid sequences within the transactivation and TAR-binding domains, and in nucleotide positions throughout the LTR. We further analyzed LTR sequences for putative transcription factor binding sites using the JASPAR transcription factor binding profile database and found several putative differences in transcription factor binding sites between R5 and X4 HIV-1 genomes, specifically identifying the C/EBP sites I and II, and Sp site III to differ with respect to sequence configuration for R5 and X4 LTRs. CONCLUSION: These observations support the hypothesis that co-receptor utilization coincides with specific genetic signatures in HIV-1 Tat and the LTR, likely due to differing transcriptional regulatory mechanisms and selective pressures applied within specific cellular targets during the course of productive HIV-1 infection.

Lasting antibody responses are mediated by a combination of newly formed and established bone marrow plasma cells drawn from clonally distinct precursors.

Current models hold that serum Ab titers are maintained chiefly by long-lived bone marrow (BM) plasma cells (PCs). In this study, we characterize the role of subpopulations of BM PCs in long-term humoral responses to T cell-dependent Ag. Surprisingly, our results indicate that 40-50% of BM PCs are recently formed cells, defined, in part, by rapid steady-state turnover kinetics and secretion of low-affinity IgM Abs. Further, for months after immunization with a hapten-protein conjugate, newly formed Ag-induced, IgM-secreting BM PCs were detected in parallel with longer-lived IgG-secreting cells, suggesting ongoing and parallel input to the BM PC pool from two distinct pools of activated B cells. Consistent with this interpretation, IgM and IgG Abs secreted by cells within distinct PC subsets exhibited distinct L chain usage. We conclude that long-term Ab responses are maintained by a dynamic BM PC pool composed of both recently formed and long-lived PCs drawn from clonally disparate precursors.

Viral CD8 T cell epitope nucleotide composition shows evidence of short- and long-term evolutionary strategies.

Viral epitopes have a distinct codon usage that reflects their dual role in infection and immunity. On the one hand, epitopes are part of proteins important to viral function; on the other hand, they are targets of the immune response. Studies of selection are most commonly based on changes of amino acid and seen through the accumulation of non-synonymous mutations. An independent measure of selection is the codon usage and underlying changeability of the nucleotide sequences. We here use multiple tools and a large-scale analysis of viral genomes to demonstrate that viral epitopes have a distinct codon usage and that this codon usage reflects distinct short- and long-term types of selection during viral evolution. We show that CD8(+) T cell epitopes are encoded by codons more distant from stop codons and more changeable than codons outside epitopes. This biased codon usage reflects the viral population toggling back and forth from a wild-type sequence to an escape mode, which enable them to avoid immune detection when needed, and go back to the functionally favorable form when the threat is removed (i.e., in a new host).

Reconstruction of regulatory networks through temporal enrichment profiling and its application to H1N1 influenza viral infection.

BACKGROUND: H1N1 influenza viruses were responsible for the 1918 pandemic that caused millions of deaths worldwide and the 2009 pandemic that caused approximately twenty thousand deaths. The cellular response to such virus infections involves extensive genetic reprogramming resulting in an antiviral state that is critical to infection control. Identifying the underlying transcriptional network driving these changes, and how this program is altered by virally-encoded immune antagonists, is a fundamental challenge in systems immunology. RESULTS: Genome-wide gene expression patterns were measured in human monocyte-derived dendritic cells (DCs) infected in vitro with seasonal H1N1 influenza A/New Caledonia/20/1999. To provide a mechanistic explanation for the timing of gene expression changes over the first 12 hours post-infection, we developed a statistically rigorous enrichment approach integrating genome-wide expression kinetics and time-dependent promoter analysis. Our approach, TIme-Dependent Activity Linker (TIDAL), generates a regulatory network that connects transcription factors associated with each temporal phase of the response into a coherent linked cascade. TIDAL infers 12 transcription factors and 32 regulatory connections that drive the antiviral response to influenza. To demonstrate the generality of this approach, TIDAL was also used to generate a network for the DC response to measles infection. The software implementation of TIDAL is freely available at http://tsb.mssm.edu/primeportal/?q=tidal_prog. CONCLUSIONS: We apply TIDAL to reconstruct the transcriptional programs activated in monocyte-derived human dendritic cells in response to influenza and measles infections. The application of this time-centric network reconstruction method in each case produces a single transcriptional cascade that recapitulates the known biology of the response with high precision and recall, in addition to identifying potentially novel antiviral factors. The ability to reconstruct antiviral networks with TIDAL enables comparative analysis of antiviral responses, such as the differences between pandemic and seasonal influenza infections.

Conserved variation: identifying patterns of stability and variability in BCR and TCR V genes with different diversity and richness metrics.

The immune system can detect most invading pathogens. The potential for detection of pathogens is dependent on the somatic diversity of the immune repertoires. While it is known that this somatic diversity is carefully generated, it is unclear how the diversity is distributed in the different genes encoding receptors of immune cells. Utilizing different metrics for richness and diversity at the level of small sequence fragments, we present here an analysis of the entire known human germline repertoire as represented by the sequences from the ImMunoGeneTics database of immune receptors. We have developed a fragment sequence quantification analysis to track variation of repertoires with different degrees of precision. Somatic diversity has previously been functionally characterized mostly by division of the V gene sequences into the more conserved and invariant framework (FR) of the receptor and more varied complementarity determining regions (CDR), that interact with the antigen. We find that CDR and FR can be explicitly identified with our sequence fragment diversity quantification technique. In terms of diversity, CDR and FR are especially distinct in B cell V genes. T cell V genes show less of the CDR/FR periodicity but are more diverse overall. Our analysis further shows that there are other areas of diversity outside the CDR and FR that are found widely dispersed in T cell receptor V genes and more tightly focused in FR1 and FR3 in the B cell receptor V genes. The diversity we observe is not dependent on allelic differences nor is this diversity segregated by individual V gene families. We would thus expect that each individual exhibit a diversity equivalent to that of the entire potential repertoire.

Establishing a relationship between bacteria in the human gut and complex regional pain syndrome.

Complex Regional Pain Syndrome (CRPS) is a serious and painful condition involving the peripheral and central nervous systems. Full comprehension of the disorder's pathophysiology remains incomplete, but research implicates the immune system as a contributor to chronic pain. Because of the impact gastrointestinal bacteria have in the development and behavior of the immune system, this study compares the GI microbial communities of 16 participants with CRPS (5 of whom have intestinal discomforts) and 16 healthy controls using 454 sequencing technology. CRPS subjects were found to have significantly less diversity than their healthy counterparts. Statistical analysis of the phylogenetic classifications revealed significantly increased levels of Proteobacteria and decreased levels of Firmicutes in CRPS subjects. Clustering analysis showed significant separation between healthy controls and CRPS subjects. These results support the hypothesis that the GI microbial communities of CRPS participants differ from those of their healthy counterparts. These variations may hold the key to understanding how CRPS develops and provide information that could yield a potential treatment.

Detecting selection in immunoglobulin sequences.

The ability to detect selection by analyzing mutation patterns in experimentally derived immunoglobulin (Ig) sequences is a critical part of many studies. Such techniques are useful not only for understanding the response to pathogens, but also to determine the role of antigen-driven selection in autoimmunity, B cell cancers and the diversification of pre-immune repertoires in certain species. Despite its importance, quantifying selection in experimentally derived sequences is fraught with difficulties. The necessary parameters for statistical tests (such as the expected frequency of replacement mutations in the absence of selection) are non-trivial to calculate, and results are not easily interpretable when analyzing more than a handful of sequences. We have developed a web server that implements our previously proposed Focused binomial test for detecting selection. Several features are integrated into the web site in order to facilitate analysis, including V(D)J germline segment identification with IMGT alignment, batch submission of sequences and integration of additional test statistics proposed by other groups. We also implement a Z-score-based statistic that increases the power of detecting selection while maintaining specificity, and further allows for the combined analysis of sequences from different germlines. The tool is freely available at http://clip.med.yale.edu/selection.