Models of Somatic Hypermutation Targeting and Substitution Based on Synonymous Mutations from High-Throughput Immunoglobulin Sequencing Data

Analyses of somatic hypermutation (SHM) patterns in B cell immunoglobulin (Ig) sequences contribute to our basic understanding of adaptive immunity, and have broad applications not only for understanding the immune response to pathogens, but also to determining the role of SHM in autoimmunity and B cell cancers. Although stochastic, SHM displays intrinsic biases that can confound statistical analysis, especially when combined with the particular codon usage and base composition in Ig sequences. Analysis of B cell clonal expansion, diversification, and selection processes thus critically depends on an accurate background model for SHM micro-sequence targeting (i.e., hot/cold-spots) and nucleotide substitution. Existing models are based on small numbers of sequences/mutations, in part because they depend on data from non-coding regions or non-functional sequences to remove the confounding influences of selection. Here, we combine high-throughput Ig sequencing with new computational analysis methods to produce improved models of SHM targeting and substitution that are based only on synonymous mutations, and are thus independent of selection. The resulting “S5F” models are based on 806,860 Synonymous mutations in 5-mer motifs from 1,145,182 Functional sequences and account for dependencies on the adjacent four nucleotides (two bases upstream and downstream of the mutation). The estimated profiles can explain almost half of the variance in observed mutation patterns, and clearly show that both mutation targeting and substitution are significantly influenced by neighboring bases. While mutability and substitution profiles were highly conserved across individuals, the variability across motifs was found to be much larger than previously estimated. The model and method source code are made available at http://clip.med.yale.edu/SH

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

SummaryThe 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

Reproducibility and reuse of adaptive immune receptor repertoire data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1-3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community\u27s founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets ([email protected])

Altered metabolism in cancer

Cancer cells have different metabolic requirements from their normal counterparts. Understanding the consequences of this differential metabolism requires a detailed understanding of glucose metabolism and its relation to energy production in cancer cells. A recent study in BMC Systems Biology by Vasquez et al. developed a mathematical model to assess some features of this altered metabolism. Here, we take a broader look at the regulation of energy metabolism in cancer cells, considering their anabolic as well as catabolic needs

Hepcidin sequesters iron to sustain nucleotide metabolism and mitochondrial function in colorectal cancer epithelial cells

Colorectal cancer (CRC) requires massive iron stores, but the complete mechanisms by which CRC modulates local iron handling are poorly understood. Here, we demonstrate that hepcidin is activated ectopically in CRC. Mice deficient in hepcidin specifically in the colon tumour epithelium, compared with wild-type littermates, exhibit significantly diminished tumour number, burden and size in a sporadic model of CRC, whereas accumulation of intracellular iron by deletion of the iron exporter ferroportin exacerbates these tumour parameters. Metabolomic analysis of three-dimensional patient-derived CRC tumour enteroids indicates a prioritization of iron in CRC for the production of nucleotides, which is recapitulated in our hepcidin/ferroportin mouse CRC models. Mechanistically, our data suggest that iron chelation decreases mitochondrial function, thereby altering nucleotide synthesis, whereas exogenous supplementation of nucleosides or aspartate partially rescues tumour growth in patient-derived enteroids and CRC cell lines in the presence of an iron chelator. Collectively, these data suggest that ectopic hepcidin in the tumour epithelium establishes an axis to sequester iron in order to maintain the nucleotide pool and sustain proliferation in colorectal tumours

Identification of Subject-Specific Immunoglobulin Alleles From Expressed Repertoire Sequencing Data

The adaptive immune receptor repertoire (AIRR) contains information on an individuals' immune past, present and potential in the form of the evolving sequences that encode the B cell receptor (BCR) repertoire. AIRR sequencing (AIRR-seq) studies rely on databases of known BCR germline variable (V), diversity (D), and joining (J) genes to detect somatic mutations in AIRR-seq data via comparison to the best-aligning database alleles. However, it has been shown that these databases are far from complete, leading to systematic misidentification of mutated positions in subsets of sample sequences. We previously presented TIgGER, a computational method to identify subject-specific V gene genotypes, including the presence of novel V gene alleles, directly from AIRR-seq data. However, the original algorithm was unable to detect alleles that differed by more than 5 single nucleotide polymorphisms (SNPs) from a database allele. Here we present and apply an improved version of the TIgGER algorithm which can detect alleles that differ by any number of SNPs from the nearest database allele, and can construct subject-specific genotypes with minimal prior information. TIgGER predictions are validated both computationally (using a leave-one-out strategy) and experimentally (using genomic sequencing), resulting in the addition of three new immunoglobulin heavy chain V (IGHV) gene alleles to the IMGT repertoire. Finally, we develop a Bayesian strategy to provide a confidence estimate associated with genotype calls. All together, these methods allow for much higher accuracy in germline allele assignment, an essential step in AIRR-seq studies

Phosphoglycerate dehydrogenase diverts glycolytic flux and contributes to oncogenesis

Most tumors exhibit increased glucose metabolism to lactate, however, the extent to which glucose-derived metabolic fluxes are used for alternative processes is poorly understood [1, 2]. Using a metabolomics approach with isotope labeling, we found that in some cancer cells a relatively large amount of glycolytic carbon is diverted into serine and glycine metabolism through phosphoglycerate dehydrogenase (PHGDH). An analysis of human cancers showed that PHGDH is recurrently amplified in a genomic region of focal copy number gain most commonly found in melanoma. Decreasing PHGDH expression impaired proliferation in amplified cell lines. Increased expression was also associated with breast cancer subtypes, and ectopic expression of PHGDH in mammary epithelial cells disrupted acinar morphogenesis and induced other phenotypic alterations that may predispose cells to transformation. Our findings show that the diversion of glycolytic flux into a specific alternate pathway can be selected during tumor development and may contribute to the pathogenesis of human cancer.National Institutes of Health (U.S.)National Cancer Institute (U.S.)Smith Family FoundationDamon Runyon Cancer Research FoundationBurroughs Wellcome Fun

Compromised fidelity of B-cell tolerance checkpoints in AChR and MuSK myasthenia gravis

Myasthenia gravis () is an autoimmune condition in which neurotransmission is impaired by binding of autoantibodies to acetylcholine receptors (hR) or, in a minority of patients, to muscle specific kinase (Mu). There are differences in the dominant IgG subclass, pathogenic mechanisms, and treatment responses between the two subtypes (hR or Mu). The antibodies are thought to be T-cell dependent, but the mechanisms underlying their production are not well understood. One aspect not previously described is whether defects in central and peripheral tolerance checkpoints, which allow autoreactive B cells to accumulate in the naive repertoire, are found in both or either form of . An established set of assays that measure the frequency of both polyreactive and autoreactive B cell receptors () in naive populations was applied to specimens collected from patients with either hR or Mu and healthy controls. Radioimmuno- and cell-based assays were used to measure binding to hR and Mu. The frequency of polyreactive and autoreactive s (n = 262) was higher in both hR and Mu patients than in healthy controls. None of the -derived s bound hR or Mu. The results indicate that both these subtypes harbor defects in central and peripheral B cell tolerance checkpoints. Defective B cell tolerance may represent a fundamental contributor to autoimmunity in and is of particular importance when considering the durability of myasthenia gravis treatment strategies, particularly biologics that eliminate B cells