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.

Salmonella Infection Drives Promiscuous B Cell Activation Followed by Extrafollicular Affinity Maturation.

The B cell response to Salmonella typhimurium (STm) occurs massively at extrafollicular sites, without notable germinal centers (GCs). Little is known in terms of its specificity. To expand the knowledge of antigen targets, we screened plasmablast (PB)-derived monoclonal antibodies (mAbs) for Salmonella specificity, using ELISA, flow cytometry, and antigen microarray. Only a small fraction (0.5%-2%) of the response appeared to be Salmonella-specific. Yet, infection of mice with limited B cell receptor (BCR) repertoires impaired the response, suggesting that BCR specificity was important. We showed, using laser microdissection, that somatic hypermutation (SHM) occurred efficiently at extrafollicular sites leading to affinity maturation that in turn led to detectable STm Ag-binding. These results suggest a revised vision of how clonal selection and affinity maturation operate in response to Salmonella. Clonal selection initially is promiscuous, activating cells with virtually undetectable affinity, yet SHM and selection occur during the extrafollicular response yielding higher affinity, detectable antibodies.

Breast cancer is marked by specific, Public T-cell receptor CDR3 regions shared by mice and humans.

The partial success of tumor immunotherapy induced by checkpoint blockade, which is not antigen-specific, suggests that the immune system of some patients contain antigen receptors able to specifically identify tumor cells. Here we focused on T-cell receptor (TCR) repertoires associated with spontaneous breast cancer. We studied the alpha and beta chain CDR3 domains of TCR repertoires of CD4 T cells using deep sequencing of cell populations in mice and applied the results to published TCR sequence data obtained from human patients. We screened peripheral blood T cells obtained monthly from individual mice spontaneously developing breast tumors by 5 months. We then looked at identical TCR sequences in published human studies; we used TCGA data from tumors and healthy tissues of 1,256 breast cancer resections and from 4 focused studies including sequences from tumors, lymph nodes, blood and healthy tissues, and from single cell dataset of 3 breast cancer subjects. We now report that mice spontaneously developing breast cancer manifest shared, Public CDR3 regions in both their alpha and beta and that a significant number of women with early breast cancer manifest identical CDR3 sequences. These findings suggest that the development of breast cancer is associated, across species, with biomarker, exclusive TCR repertoires.

Hierarchical Clustering Can Identify B Cell Clones with High Confidence in Ig Repertoire Sequencing Data.

Adaptive immunity is driven by the expansion, somatic hypermutation, and selection of B cell clones. Each clone is the progeny of a single B cell responding to Ag, with diversified Ig receptors. These receptors can now be profiled on a large scale by next-generation sequencing. Such data provide a window into the microevolutionary dynamics that drive successful immune responses and the dysregulation that occurs with aging or disease. Clonal relationships are not directly measured, but they must be computationally inferred from these sequencing data. Although several hierarchical clustering-based methods have been proposed, they vary in distance and linkage methods and have not yet been rigorously compared. In this study, we use a combination of human experimental and simulated data to characterize the performance of hierarchical clustering-based methods for partitioning sequences into clones. We find that single linkage clustering has high performance, with specificity, sensitivity, and positive predictive value all >99%, whereas other linkages result in a significant loss of sensitivity. Surprisingly, distance metrics that incorporate the biases of somatic hypermutation do not outperform simple Hamming distance. Although errors were more likely in sequences with short junctions, using the entire dataset to choose a single distance threshold for clustering is near optimal. Our results suggest that hierarchical clustering using single linkage with Hamming distance identifies clones with high confidence and provides a fully automated method for clonal grouping. The performance estimates we develop provide important context to interpret clonal analysis of repertoire sequencing data and allow for rigorous testing of other clonal grouping algorithms.

Dysregulation of B Cell Repertoire Formation in Myasthenia Gravis Patients Revealed through Deep Sequencing.

Myasthenia gravis (MG) is a prototypical B cell-mediated autoimmune disease affecting 20-50 people per 100,000. The majority of patients fall into two clinically distinguishable types based on whether they produce autoantibodies targeting the acetylcholine receptor (AChR-MG) or muscle specific kinase (MuSK-MG). The autoantibodies are pathogenic, but whether their generation is associated with broader defects in the B cell repertoire is unknown. To address this question, we performed deep sequencing of the BCR repertoire of AChR-MG, MuSK-MG, and healthy subjects to generate ∼518,000 unique VH and VL sequences from sorted naive and memory B cell populations. AChR-MG and MuSK-MG subjects displayed distinct gene segment usage biases in both VH and VL sequences within the naive and memory compartments. The memory compartment of AChR-MG was further characterized by reduced positive selection of somatic mutations in the VH CDR and altered VH CDR3 physicochemical properties. The VL repertoire of MuSK-MG was specifically characterized by reduced V-J segment distance in recombined sequences, suggesting diminished VL receptor editing during B cell development. Our results identify large-scale abnormalities in both the naive and memory B cell repertoires. Particular abnormalities were unique to either AChR-MG or MuSK-MG, indicating that the repertoires reflect the distinct properties of the subtypes. These repertoire abnormalities are consistent with previously observed defects in B cell tolerance checkpoints in MG, thereby offering additional insight regarding the impact of tolerance defects on peripheral autoimmune repertoires. These collective findings point toward a deformed B cell repertoire as a fundamental component of MG.

A Model of Somatic Hypermutation Targeting in Mice Based on High-Throughput Ig Sequencing Data.

Analyses of somatic hypermutation (SHM) patterns in B cell Ig sequences have important basic science and clinical applications, but they are often confounded by the intrinsic biases of SHM targeting on specific DNA motifs (i.e., hot and cold spots). Modeling these biases has been hindered by the difficulty in identifying mutated Ig sequences in vivo in the absence of selection pressures, which skew the observed mutation patterns. To generate a large number of unselected mutations, we immunized B1-8 H chain transgenic mice with nitrophenyl to stimulate nitrophenyl-specific λ+ germinal center B cells and sequenced the unexpressed κ L chains using next-generation methods. Most of these κ sequences had out-of-frame junctions and were presumably uninfluenced by selection. Despite being nonfunctionally rearranged, they were targeted by SHM and displayed a higher mutation frequency than functional sequences. We used 39,173 mutations to construct a quantitative SHM targeting model. The model showed targeting biases that were consistent with classic hot and cold spots, yet revealed additional highly mutable motifs. We observed comparable targeting for functional and nonfunctional sequences, suggesting similar biological processes operate at both loci. However, we observed species- and chain-specific targeting patterns, demonstrating the need for multiple SHM targeting models. Interestingly, the targeting of C/G bases and the frequency of transition mutations at C/G bases was higher in mice compared with humans, suggesting lower levels of DNA repair activity in mice. Our models of SHM targeting provide insights into the SHM process and support future analyses of mutation patterns.

High-resolution antibody dynamics of vaccine-induced immune responses.

The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.

Integrating B cell lineage information into statistical tests for detecting selection in Ig sequences.

Detecting selection in B cell Ig sequences is critical to understanding affinity maturation and can provide insights into Ag-driven selection in normal and pathologic immune responses. The most common sequence-based methods for detecting selection analyze the ratio of replacement and silent mutations using a binomial statistical analysis. However, these approaches have been criticized for low sensitivity. An alternative method is based on the analysis of lineage trees constructed from sets of clonally related Ig sequences. Several tree shape measures have been proposed as indicators of selection that can be statistically compared across cohorts. However, we show that tree shape analysis is confounded by underlying experimental factors that are difficult to control for in practice, including the sequencing depth and number of generations in each clone. Thus, although lineage tree shapes may reflect selection, their analysis alone is an unreliable measure of in vivo selection. To usefully capture the information provided by lineage trees, we propose a new method that applies the binomial statistical framework to mutations identified based on lineage tree structure. This hybrid method is able to detect selection with increased sensitivity in both simulated and experimental data sets. We anticipate that this approach will be especially useful in the analysis of large-scale Ig sequencing data sets generated by high-throughput sequencing technologies.

Systematic identification of edited microRNAs in the human brain.

Adenosine-to-inosine (A-to-I) editing modifies RNA transcripts from their genomic blueprint. A prerequisite for this process is a double-stranded RNA (dsRNA) structure. Such dsRNAs are formed as part of the microRNA (miRNA) maturation process, and it is therefore expected that miRNAs are affected by A-to-I editing. Editing of miRNAs has the potential to add another layer of complexity to gene regulation pathways, especially if editing occurs within the miRNA-mRNA recognition site. Thus, it is of interest to study the extent of this phenomenon. Current reports in the literature disagree on its extent; while some reports claim that it may be widespread, others deem the reported events as rare. Utilizing a next-generation sequencing (NGS) approach supplemented by an extensive bioinformatic analysis, we were able to systematically identify A-to-I editing events in mature miRNAs derived from human brain tissues. Our algorithm successfully identified many of the known editing sites in mature miRNAs and revealed 17 novel human sites, 12 of which are in the recognition sites of the miRNAs. We confirmed most of the editing events using in vitro ADAR overexpression assays. The editing efficiency of most sites identified is very low. Similar results are obtained for publicly available data sets of mouse brain-regions tissues. Thus, we find that A-to-I editing does alter several miRNAs, but it is not widespread.

pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires.

UNLABELLED: Driven by dramatic technological improvements, large-scale characterization of lymphocyte receptor repertoires via high-throughput sequencing is now feasible. Although promising, the high germline and somatic diversity, especially of B-cell immunoglobulin repertoires, presents challenges for analysis requiring the development of specialized computational pipelines. We developed the REpertoire Sequencing TOolkit (pRESTO) for processing reads from high-throughput lymphocyte receptor studies. pRESTO processes raw sequences to produce error-corrected, sorted and annotated sequence sets, along with a wealth of metrics at each step. The toolkit supports multiplexed primer pools, single- or paired-end reads and emerging technologies that use single-molecule identifiers. pRESTO has been tested on data generated from Roche and Illumina platforms. It has a built-in capacity to parallelize the work between available processors and is able to efficiently process millions of sequences generated by typical high-throughput projects. AVAILABILITY AND IMPLEMENTATION: pRESTO is freely available for academic use. The software package and detailed tutorials may be downloaded from http://clip.med.yale.edu/presto.

The effects of somatic hypermutation on neutralization and binding in the PGT121 family of broadly neutralizing HIV antibodies.

Broadly neutralizing HIV antibodies (bnAbs) are typically highly somatically mutated, raising doubts as to whether they can be elicited by vaccination. We used 454 sequencing and designed a novel phylogenetic method to model lineage evolution of the bnAbs PGT121-134 and found a positive correlation between the level of somatic hypermutation (SHM) and the development of neutralization breadth and potency. Strikingly, putative intermediates were characterized that show approximately half the mutation level of PGT121-134 but were still capable of neutralizing roughly 40-80% of PGT121-134 sensitive viruses in a 74-virus panel at median titers between 15- and 3-fold higher than PGT121-134. Such antibodies with lower levels of SHM may be more amenable to elicitation through vaccination while still providing noteworthy coverage. Binding characterization indicated a preference of inferred intermediates for native Env binding over monomeric gp120, suggesting that the PGT121-134 lineage may have been selected for binding to native Env at some point during maturation. Analysis of glycan-dependent neutralization for inferred intermediates identified additional adjacent glycans that comprise the epitope and suggests changes in glycan dependency or recognition over the course of affinity maturation for this lineage. Finally, patterns of neutralization of inferred bnAb intermediates suggest hypotheses as to how SHM may lead to potent and broad HIV neutralization and provide important clues for immunogen design.

Barcoding bias in high-throughput multiplex sequencing of miRNA.

Second-generation sequencing is gradually becoming the method of choice for miRNA detection and expression profiling. Given the relatively small number of miRNAs and improvements in DNA sequencing technology, studying miRNA expression profiles of multiple samples in a single flow cell lane becomes feasible. Multiplexing strategies require marking each miRNA library with a DNA barcode. Here we report that barcodes introduced through adapter ligation confer significant bias on miRNA expression profiles. This bias is much higher than the expected Poisson noise and masks significant expression differences between miRNA libraries. This bias can be eliminated by adding barcodes during PCR amplification of libraries. The accuracy of miRNA expression measurement in multiplexed experiments becomes a function of sample number.

Transcriptional profiling of CD4 T cells identifies distinct subgroups of HIV-1 elite controllers.

Human immunodeficiency virus type 1 (HIV-1) elite controllers maintain undetectable levels of viral replication in the absence of antiretroviral therapy (ART), but their underlying immunological and virological characteristics may vary. Here, we used a whole-genome transcriptional profiling approach to characterize gene expression signatures of CD4 T cells from an unselected cohort of elite controllers. The transcriptional profiles for the majority of elite controllers were similar to those of ART-treated patients but different from those of HIV-1-negative persons. Yet, a smaller proportion of elite controllers showed an alternative gene expression pattern that was indistinguishable from that of HIV-1-negative persons but different from that of highly active antiretroviral therapy (HAART)-treated individuals. Elite controllers with the latter gene expression signature had significantly higher CD4 T cell counts and lower levels of HIV-1-specific CD8(+) T cell responses but did not significantly differ from other elite controllers in terms of HLA class I alleles, HIV-1 viral loads determined by ultrasensitive single-copy PCR assays, or chemokine receptor polymorphisms. Thus, these data identify a specific subgroup of elite controllers whose immunological and gene expression characteristics approximate those of HIV-1-negative persons.

Identification of differentially expressed genes in uveal melanoma using suppressive subtractive hybridization.

PURPOSE: Uveal melanoma (UM) is the most common primary cancer of the eye, resulting not only in vision loss, but also in metastatic death. This study attempts to identify changes in the patterns of gene expression that lead to malignant transformation and proliferation of normal uveal melanocytes (UVM) using the Suppressive Subtractive Hybridization (SSH) technique. METHODS: The SSH technique was used to isolate genes that are differentially expressed in the TP31 cell line derived from a primary UM compared to UVM. The expression level of selected genes was further validated by microarray, semi-quantitative RT-PCR and western blot analyses. RESULTS: Analysis of the subtracted libraries revealed that 37 and 36 genes were, respectively, up- and downregulated in TP31 cells compared to UVM. Differential expression of the majority of these genes was confirmed by comparing UM cells with UVM by microarray. The expression pattern of selected genes was analyzed by semi-quantitative RT-PCR and western blot, and was found to be consistent with the SSH findings. CONCLUSIONS: We demonstrated that the SSH technique is efficient to detect differentially expressed genes in UM. The genes identified in this study represent valuable candidates for further functional analysis in UM and should be informative in studying the biology of this tumor.

In cellulo DNA analysis (LMPCR footprinting).

The in cellulo analysis of DNA protein interactions and chromatin structure is very important to better understand the mechanisms involved in the regulation of gene expression. The nuclease-hypersensitive sites and sequences bound by transcription factors often correspond to genetic regulatory elements. Using the Ligation-mediated polymerase chain reaction (LMPCR) technology, it is possible to precisely analyze these DNA sequences to demonstrate the existence of DNA-protein interactions or unusual DNA structures directly in living cells. Indeed, the ideal chromatin substrate is, of course, found inside intact cells. LMPCR, a genomic-sequencing, technique that map DNA single-strand breaks at the sequence level of resolution, is the method of choice for in cellulo footprinting and DNA structure studies because it can be used to investigate any complex genomes, including human. The detailed conventional and automated LMPCR protocols are presented in this chapter.

Single T Cell Sequencing Demonstrates the Functional Role of αß TCR Pairing in Cell Lineage and Antigen Specificity.

Although structural studies of individual T cell receptors (TCRs) have revealed important roles for both the α and ß chain in directing MHC and antigen recognition, repertoire-level immunogenomic analyses have historically examined the ß chain alone. To determine the amount of useful information about TCR repertoire function encoded within αß pairings, we analyzed paired TCR sequences from nearly 100,000 unique CD4+ and CD8+ T cells captured using two different high-throughput, single-cell sequencing approaches. Our results demonstrate little overlap in the healthy CD4+ and CD8+ repertoires, with shared TCR sequences possessing significantly shorter CDR3 sequences corresponding to higher generation probabilities. We further utilized tools from information theory and machine learning to show that while α and ß chains are only weakly associated with lineage, αß pairings appear to synergistically drive TCR-MHC interactions. Vαß gene pairings were found to be the TCR feature most informative of T cell lineage, supporting the existence of germline-encoded paired αß TCR-MHC interaction motifs. Finally, annotating our TCR pairs using a database of sequences with known antigen specificities, we demonstrate that approximately a third of the T cells possess α and ß chains that each recognize different known antigens, suggesting that αß pairing is critical for the accurate inference of repertoire functionality. Together, these findings provide biological insight into the functional implications of αß pairing and highlight the utility of single-cell sequencing in immunogenomics.

Mosaic deletion patterns of the human antibody heavy chain gene locus shown by Bayesian haplotyping.

Analysis of antibody repertoires by high-throughput sequencing is of major importance in understanding adaptive immune responses. Our knowledge of variations in the genomic loci encoding immunoglobulin genes is incomplete, resulting in conflicting VDJ gene assignments and biased genotype and haplotype inference. Haplotypes can be inferred using IGHJ6 heterozygosity, observed in one third of the people. Here, we propose a robust novel method for determining VDJ haplotypes by adapting a Bayesian framework. Our method extends haplotype inference to IGHD- and IGHV-based analysis, enabling inference of deletions and copy number variations in the entire population. To test this method, we generated a multi-individual data set of naive B-cell repertoires, and found allele usage bias, as well as a mosaic, tiled pattern of deleted IGHD and IGHV genes. The inferred haplotypes may have clinical implications for genetic disease predispositions. Our findings expand the knowledge that can be extracted from antibody repertoire sequencing data.

Control of integrin genes expression in the eye.

The ending race for sequencing the human genome has left the scientists faced with new challenges. Indeed, now that almost every human gene has been sequenced and precisely positioned on the human chromosomes, one of the next, most burning task consist in understanding how transcription of these genes is ensured in any given cell. The integrins encoding genes are no exception. Integrins bridge the cell to the many components from the extracellular matrix (ECM), such as laminins (LM) and collagens, and thereby transduce intracellular signals that will alter many of the cell's properties such as adhesion, migration, proliferation and survival. As a much clearer picture of the many proteins that belong to this family has emerged over the last few years, tremendous efforts have been dedicated to the identification of the regulatory sequences that modulate their expression. This review provides an overview of the current state of knowledge about the organization of the regulatory elements and the transcription factors (TF) they bind that are used by the cell in order to ensure transcription of each of the integrins gene. A particular attention has been given to those reported to be expressed in the eye. It also explores how components from the ECM might participate in the control of integrins gene expression and establishes links to wound healing of the corneal epithelium, a process that transiently alter the composition of the basement membrane on which the epithelial cells lie.

Transcriptional regulation of the cyclin-dependent kinase inhibitor 1A (p21) gene by NFI in proliferating human cells.

The cyclin-dependent kinase inhibitor 1A (CDKN1A), also known as p21 (WAF1/CIP1) modulates cell cycle, apoptosis, senescence and differentiation via specific protein-protein interactions with the cyclins, cyclin-dependent kinase (Cdk), and many others. Expression of the p21 gene is mainly regulated at the transcriptional level. By conducting both ligation-mediated PCR (LMPCR) and chromatin immunoprecipitation (ChIP) in vivo, we identified a functional target site for the transcription factor, nuclear factor I (NFI), in the basal promoter from the p21 gene. Transfection of recombinant constructs bearing mutations in the p21 NFI site demonstrated that NFI acts as a repressor of p21 gene expression in various types of cultured cells. Inhibition of NFI in human skin fibroblasts through RNAi considerably increased p21 promoter activity suggesting that NFI is a key repressor of p21 transcription. Over-expression of each of the four NFI isoforms in HCT116 cells established that each of them contribute to various extend to the repression of the p21 gene. Most of all, over-expression of NFI-B in doxorubicin, growth-arrested HCT116 increased the proportion of cells in the S-phase of the cell cycle whereas NFI-A and NFI-X reduced it, thereby establishing a role for NFI in the cell cycle dependent expression of p21.

Antibody Repertoire Analysis of Hepatitis C Virus Infections Identifies Immune Signatures Associated With Spontaneous Clearance.

Hepatitis C virus (HCV) is a major public health concern, with over 70 million people infected worldwide, who are at risk for developing life-threatening liver disease. No vaccine is available, and immunity against the virus is not well-understood. Following the acute stage, HCV usually causes chronic infections. However, ~30% of infected individuals spontaneously clear the virus. Therefore, using HCV as a model for comparing immune responses between spontaneous clearer (SC) and chronically infected (CI) individuals may empower the identification of mechanisms governing viral infection outcomes. Here, we provide the first in-depth analysis of adaptive immune receptor repertoires in individuals with current or past HCV infection. We demonstrate that SC individuals, in contrast to CI patients, develop clusters of antibodies with distinct properties. These antibodies' characteristics were used in a machine learning framework to accurately predict infection outcome. Using combinatorial antibody phage display library technology, we identified HCV-specific antibody sequences. By integrating these data with the repertoire analysis, we constructed two antibodies characterized by high neutralization breadth, which are associated with clearance. This study provides insight into the nature of effective immune response against HCV and demonstrates an innovative approach for constructing antibodies correlating with successful infection clearance. It may have clinical implications for prognosis of the future status of infection, and the design of effective immunotherapies and a vaccine for HCV.