Search

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 post-treatment 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.

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.One Sentence SummaryFood allergic sensitization is associated with altered B cell development in children with skin barrier dysfunction.

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.

Kids armed with anti-coronavirus B cells It remains unclear whether B cell repertoires against coronaviruses and other pathogens differ between adults and children and how important these distinctions are. Yang et al. analyzed blood samples from young children and adults, as well as tissues from deceased organ donors, characterizing the B cell receptor (BCR) repertoires specific to six common pathogens and two viruses that they had not seen before: Ebola virus and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Children had higher frequencies of B cells with convergent BCR heavy chains against previously encountered pathogens and higher frequencies of class-switched convergent B cell clones against SARS-CoV-2 and related coronaviruses. These findings suggest that encounters with coronaviruses in early life may produce cross-reactive memory B cell populations that contribute to divergent COVID-19 susceptibilities. Science, this issue p. 738, Blood taken from children before the COVID-19 pandemic contains memory B cells that bind SARS-CoV-2., 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 antigenspecific 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.

Clinical diagnoses rely on a wide variety of laboratory tests and imaging studies, interpreted alongside physical examination findings and the patient’s history and symptoms. Currently, the tools of diagnosis make limited use of the immune system’s internal record of specific disease exposures encoded by the antigen-specific receptors of memory B cells and T cells, and there has been little integration of the combined information from B cell and T cell receptor sequences. Here, we analyze extensive receptor sequence datasets with three different machine learning representations of immune receptor repertoires to develop an interpretive framework,MAchine Learning for Immunological Diagnosis (Mal-ID), that screens for multiple illnesses simultaneously. This approach is effective in identifying a variety of disease states, including acute and chronic infections and autoimmune disorders. It is able to do so even when there are other differences present in the immune repertoires, such as between pediatric or adult patient groups. Importantly, many features of the model of immune receptor sequences are human-interpretable. They independently recapitulate known biology of the responses to infection by SARS-CoV-2 and HIV, provide evidence of receptor antigen specificity, and reveal common features of autoreactive immune receptor repertoires, indicating that machine learning on immune repertoires can yield new immunological knowledge. This framework could be useful in identifying immune responses to new infectious diseases as they emerge.

Plasma cells (PCs) are the major source of pathogenic allo- and autoantibodies and have historically demonstrated resistance to therapeutic targeting. However, significant recent clinical progress has been made with the use of second-generation proteasome inhibitors (PIs). PIs provide efficient elimination of plasmablast-mediated humoral responses; however, long-lived bone marrow (BM) resident PCs (LLPCs) demonstrate therapeutic resistance, particularly to first-generation PIs. In addition, durability of antibody (Ab) reduction still requires improvement. More recent clinical trials have focused on conditions mediated by LLPCs and have included mechanistic studies of LLPCs from PI-treated patients. A recent clinical trial of carfilzomib (a second-generation irreversible PI) demonstrated improved efficacy in eliminating BM PCs and reducing anti-HLA Abs in chronically HLA-sensitized patients; however, Ab rebound was observed over several weeks to months following PI therapy. Importantly, recent murine studies have provided substantial insights into PC biology, thereby further enhancing our understanding of PC populations. It is now clear that BMPC populations, where LLPCs are thought to primarily reside, are heterogeneous and have distinct gene expression, metabolic, and survival signatures that enable identification and characterization of PC subsets. This review highlights recent advances in PC biology and clinical trials in transplant populations.

RationaleWhile nasal brushing transcriptomics can identify disease subtypes in chronic pulmonary diseases, it is unknown whether this is true in pediatric ARDS (PARDS).ObjectivesDetermine whether nasal transcriptomics and methylomics can identify clinically meaningful PARDS subgroups that reflect important pathobiological processes.MethodsNasal 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 ResultsTwenty-four control and 39 PARDS subjects were enrolled. Two nasal methylation patterns were identified. Methyl Subgroup 1 had hypomethylation of cell adhesion genes and had only PARDS subjects. Methyl Subgroup 2 contained PARDS and control subjects with hypomethylation of cell metabolism genes. Neither methylation pattern had longer PARDS duration. 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). Control specimens were largely Subgroup C. In comparison with 17 serum biomarkers, the nasal transcriptome was more predictive of prolonged PARDS. Initial Transcriptomic Subgroup B or D subjects 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. ConclusionsPARDS 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.RegistrationClinicaltrials.gov NCT#03539783 May 29, 2018

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.

Vaccination and infection promote the formation, tissue distribution, and clonal evolution of B cells encoding humoral immune memory. We evaluated convergent antigen-specific antibody genes of similar sequences shared between individuals in pediatric and adult blood, and deceased organ donor tissues. B cell memory varied for different pathogens. Polysaccharide antigen-specific clones were not exclusive to the spleen. Adults’ convergent clones often express mutated IgM or IgD in blood and are class-switched in lymphoid tissues; in contrast, children have abundant class-switched convergent clones in blood. Consistent with serological reports, pre-pandemic children had class-switched convergent clones to SARS-CoV-2, enriched in cross-reactive clones for seasonal coronaviruses, while adults showed few such clones in blood or lymphoid tissues. These results extend age-related and anatomical mapping of human humoral pathogen-specific immunity.One Sentence SummaryChildren have elevated frequencies of pathogen-specific class-switched memory B cells, including SARS-CoV-2-binding clones.

Background HIV associated neurocognitive disorders cause significant morbidity and mortality despite the advent of highly active antiretroviral therapy. A deeper understanding of fundamental mechanisms underlying HIV infection and pathogenesis in the central nervous system is warranted. Microglia are resident myeloid cells of the brain that are readily infected by HIV and may constitute a CNS reservoir. We evaluated two microglial model cell lines (C20, HMC3) and two sources of primary cell-derived microglia (monocyte-derived microglia [MMG] and induced pluripotent stem cell-derived microglia [iPSC-MG]) as potential model systems for studying HIV-microglia interactions. Results All four microglial model cells expressed typical myeloid markers with the exception of low or absent CD45 and CD11b expression by C20 and HMC3, and all four expressed the microglia-specific markers P2RY12 and TMEM119. Marked differences were observed upon gene expression profiling, however, indicating that MMG and iPSC-MG cluster closely together with primary human microglial cells, while C20 and HMC3 were similar to each other but very different from primary microglia. Expression of HIV-relevant genes also revealed important differences, with iPSC-MG and MMG expressing relevant genes at levels more closely resembling primary microglia. iPSC-MG and MMG were readily infected with R5-tropic HIV, while C20 and HMC3 lack CD4 and require pseudotyping for infection. Despite many similarities, HIV replication dynamics and HIV-1 particle capture by Siglec-1 differed markedly between the MMG and iPSC-MG. Conclusions MMG and iPSC-MG appear to be viable microglial models that are susceptible to HIV infection and bear more similarities to authentic microglia than two transformed microglia cell lines. The observed differences in HIV replication and particle capture between MMG and iPSC-MG warrant further study.

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.

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., A subset of HIV-infected individuals can develop broadly neutralizing antibodies. Boyd and colleagues studied such HIV-infected individuals and found significant perturbations in their antibody repertoires, including increased frequencies of B cells expressing antibodies with features associated with autoreactivity.

The development of obesity-associated insulin resistance is associated with B-lymphocyte accumulation in visceral adipose tissue (VAT) and is prevented by B-cell ablation. To characterize potentially pathogenic B-cell repertoires in this disorder, we performed high-throughput immunoglobulin (Ig) sequencing from multiple tissues of mice fed high-fat diet (HFD) and regular diet (RD). HFD significantly changed the biochemical properties of Ig heavy-chain complementarity-determining region-3 (CDRH3) sequences, selecting for IgA antibodies with shorter and more hydrophobic CDRH3 in multiple tissues. A set of convergent antibodies of highly similar sequences found in the VAT of HFD mice but not RD mice showed significant somatic mutation, suggesting a response shared between mice to a common antigen or antigens. These findings indicate that a simple high-fat dietary intervention has a major impact on mouse B-cell repertoires, particularly in adipose tissues.

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.

Details about IgE-producing B cells in the gut in the context of food allergy are scarce, despite the frequent exposure of the gut and its associated lymphoid tissues to dietary antigens. A new study finds that IgE-producing B cells are enriched in gut tissues and are probably generated from local antibody isotype switching.

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 IgG 4 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 IgG 4 .

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.

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.

SummaryB cells produce a diverse antibody repertoire by undergoing gene rearrangements. Pathogen exposure induces the clonal expansion of B cells expressing antibodies that can bind the infectious agent. To assess human B cell responses to trivalent seasonal influenza and monovalent pandemic H1N1 vaccination, we sequenced gene rearrangements encoding the immunoglobulin heavy chain, a major determinant of epitope recognition. The magnitude of B cell clonal expansions correlates with an individual’s secreted antibody response to the vaccine, and the expanded clones are enriched with those expressing influenza-specific monoclonal antibodies. Additionally, B cell responses to pandemic influenza H1N1 vaccination and infection in different people show a prominent family of convergent antibody heavy chain gene rearrangements specific to influenza antigens. These results indicate that microbes can induce specific signatures of immunoglobulin gene rearrangements and that pathogen exposure can potentially be assessed from B cell repertoires.

Elderly humans show decreased humoral immunity to pathogens and vaccines, yet the effects of aging on B cells are not fully known. Chronic viral infection by cytomegalovirus (CMV) is implicated as a driver of clonal T cell proliferations in some aging humans, but whether CMV or Epstein-Barr virus (EBV) infection contributes to alterations in the B cell repertoire with age is unclear. We have used high-throughput DNA sequencing of immunoglobulin heavy chain (IGH) gene rearrangements to study the B cell receptor repertoires over two successive years in 27 individuals ranging in age from 20 to 89 years. Some features of the B cell repertoire remain stable with age, but elderly subjects show increased numbers of B cells with long CDR3 regions, a trend toward accumulation of more highly mutated IgM and IgG immunoglobulin genes, and persistent clonal B cell populations in the blood. Seropositivity for CMV or EBV infection alters B cell repertoires, regardless of the individual's age: EBV infection correlates with the presence of persistent clonal B cell expansions, while CMV infection correlates with the proportion of highly mutated antibody genes. These findings isolate effects of aging from those of chronic viral infection on B cell repertoires, and provide a baseline for understanding human B cell responses to vaccination or infectious stimuli.

Current human immunodeficiency virus-1 (HIV-1) vaccines elicit strain-specific neutralizing antibodies. However, cross-reactive neutralizing antibodies arise in approximately 20% of HIV-1-infected individuals, and details of their generation could provide a blueprint for effective vaccination. Here we report the isolation, evolution and structure of a broadly neutralizing antibody from an African donor followed from the time of infection. The mature antibody, CH103, neutralized approximately 55% of HIV-1 isolates, and its co-crystal structure with the HIV-1 envelope protein gp120 revealed a new loop-based mechanism of CD4-binding-site recognition. Virus and antibody gene sequencing revealed concomitant virus evolution and antibody maturation. Notably, the unmutated common ancestor of the CH103 lineage avidly bound the transmitted/founder HIV-1 envelope glycoprotein, and evolution of antibody neutralization breadth was preceded by extensive viral diversification in and near the CH103 epitope. These data determine the viral and antibody evolution leading to induction of a lineage of HIV-1 broadly neutralizing antibodies, and provide insights into strategies to elicit similar antibodies by 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.

Common variable immune deficiency (CVID) is the most common symptomatic primary immune deficiency, affecting ~1 in 25,000 persons. These patients suffer from impaired antibody responses, autoimmunity, and susceptibility to lymphoid cancers. To explore the cellular basis for these clinical phenotypes, we conducted high-throughput DNA sequencing of immunoglobulin heavy chain gene rearrangements from 93 CVID patients and 105 control subjects and sorted naïve and memory B cells from 13 of the CVID patients and 10 of the control subjects. The CVID patients showed abnormal VDJ rearrangement and abnormal formation of complementarity-determining region 3 (CDR3). We observed a decreased selection against antibodies with long CDR3s in memory repertoires and decreased variable gene replacement, offering possible mechanisms for increased patient autoreactivity. Our data indicate that patient immunodeficiency might derive from both decreased diversity of the naïve B cell pool and decreased somatic hypermutation in memory repertoires. The CVID patients also exhibited an abnormal clonal expansion of unmutated B cells relative to the controls. Although impaired B cell germinal center activation is commonly viewed as causative in CVID, these data indicate that CVID B cells diverge from controls as early as the pro-B stage, cell and suggest possible explanations for the increased incidence of autoimmunity, immunodeficiency, and lymphoma CVID patients.

The human and mouse antibody repertoires are formed by identical processes, but like all small animals, mice only have sufficient lymphocytes to express a small part of the potential antibody repertoire. In this study, we determined how the heavy chain repertoires of two mouse strains are generated. Analysis of IgM- and IgG-associated VDJ rearrangements generated by high-throughput sequencing confirmed the presence of 99 functional immunoglobulin heavy chain variable (IGHV) genes in the C57BL/6 genome, and inferred the presence of 164 IGHV genes in the BALB/c genome. Remarkably, only five IGHV sequences were common to both strains. Compared with humans, little N nucleotide addition was seen in the junctions of mouse VDJ genes. Germline human IgG-associated IGHV genes are rare, but many murine IgG-associated IGHV genes were unmutated. Together these results suggest that the expressed mouse repertoire is more germline-focused than the human repertoire. The apparently divergent germline repertoires of the mouse strains are discussed with reference to reports that inbred mouse strains carry blocks of genes derived from each of the three subspecies of the house mouse. We hypothesize that the germline genes of BALB/c and C57BL/6 mice may originally have evolved to generate distinct germline-focused antibody repertoires in the different mouse subspecies.

We define a “threaded blockset,” which is a novel generalization of the classic notion of a multiple alignment. A new computer program called TBA (for “threaded blockset aligner”) builds a threaded blockset under the assumption that all matching segments occur in the same order and orientation in the given sequences; inversions and duplications are not addressed. TBA is designed to be appropriate for aligning many, but by no means all, megabase-sized regions of multiple mammalian genomes. The output of TBA can be projected onto any genome chosen as a reference, thus guaranteeing that different projections present consistent predictions of which genomic positions are orthologous. This capability is illustrated using a new visualization tool to view TBA-generated alignments of vertebrate Hox clusters from both the mammalian and fish perspectives. Experimental evaluation of alignment quality, using a program that simulates evolutionary change in genomic sequences, indicates that TBA is more accurate than earlier programs. To perform the dynamic-programming alignment step, TBA runs a stand-alone program called MULTIZ, which can be used to align highly rearranged or incompletely sequenced genomes. We describe our use of MULTIZ to produce the whole-genome multiple alignments at the Santa Cruz Genome Browser.

The University of California Santa Cruz (UCSC) Table Browser (http://genome.ucsc.edu/cgi-bin/hgText) provides text-based access to a large collection of genome assemblies and annotation data stored in the Genome Browser Database. A flexible alternative to the graphical-based Genome Browser, this tool offers an enhanced level of query support that includes restrictions based on field values, free-form SQL queries and combined queries on multiple tables. Output can be filtered to restrict the fields and lines returned, and may be organized into one of several formats, including a simple tab- delimited file that can be loaded into a spreadsheet or database as well as advanced formats that may be uploaded into the Genome Browser as custom annotation tracks. The Table Browser User's Guide located on the UCSC website provides instructions and detailed examples for constructing queries and configuring output.

The University of California Santa Cruz (UCSC) Genome Browser Database is an up to date source for genome sequence data integrated with a large collection of related annotations. The database is optimized to support fast interactive performance with the web-based UCSC Genome Browser, a tool built on top of the database for rapid visualization and querying of the data at many levels. The annotations for a given genome are displayed in the browser as a series of tracks aligned with the genomic sequence. Sequence data and annotations may also be viewed in a text-based tabular format or downloaded as tab-delimited flat files. The Genome Browser Database, browsing tools and downloadable data files can all be found on the UCSC Genome Bioinformatics website (http://genome.ucsc.edu), which also contains links to documentation and related technical information.

SummaryMonoclonal antibodies derived from blood plasma cells of acute HIV-1-infected individuals are predominantly targeted to the HIV Env gp41 and cross-reactive with commensal bacteria. To understand this phenomenon, we examined anti-HIV responses in ileum B cells using recombinant antibody technology and probed their relationship to commensal bacteria. The dominant ileum B cell response was to Env gp41. Remarkably, a majority (82%) of the ileum anti-gp41 antibodies cross-reacted with commensal bacteria, and of those, 43% showed non-HIV-1 antigen polyreactivity. Pyrosequencing revealed shared HIV-1 antibody clonal lineages between ileum and blood. Mutated immunoglobulin G antibodies cross-reactive with both Env gp41 and microbiota could also be isolated from the ileum of HIV-1 uninfected individuals. Thus, the gp41 commensal bacterial antigen cross-reactive antibodies originate in the intestine, and the gp41 Env response in HIV-1 infection can be derived from a preinfection memory B cell pool triggered by commensal bacteria that cross-react with Env.

In order to identify novel somatic mutations associated with classic BCR/ABL1-negative myeloproliferative neoplasms, we performed high-coverage genome sequencing of DNA from peripheral blood granulocytes and cultured skin fibroblasts from a patient with MPL W515K-positive primary myelofibrosis. The primary myelofibrosis genome had a low somatic mutation rate, consistent with that observed in similar hematopoietic tumor genomes. Interfacing of whole-genome DNA sequence data with RNA expression data identified three somatic mutations of potential functional significance: i) a nonsense mutation in CARD6, implicated in modulation of NF-kappaB activation; ii) a 19-base pair deletion involving a potential regulatory region in the 5'-untranslated region of BRD2, implicated in transcriptional regulation and cell cycle control; and iii) a non-synonymous point mutation in KIAA0355, an uncharacterized protein. Additional mutations in three genes (CAP2, SOX30, and MFRP) were also evident, albeit with no support for expression at the RNA level. Re-sequencing of these six genes in 178 patients with polycythemia vera, essential thrombocythemia, and myelofibrosis did not identify recurrent somatic mutations in these genes. Finally, we describe methods for reducing false-positive variant calls in the analysis of hematologic malignancies with a low somatic mutation rate. This trial is registered with ClinicalTrials.gov (NCT01108159).

SummaryDengue is the most prevalent mosquito-borne viral disease in humans, and the lack of early prognostics, vaccines, and therapeutics contributes to immense disease burden. To identify patterns that could be used for sequence-based monitoring of the antibody response to dengue, we examined antibody heavy-chain gene rearrangements in longitudinal peripheral blood samples from 60 dengue patients. Comparing signatures between acute dengue, postrecovery, and healthy samples, we found increased expansion of B cell clones in acute dengue patients, with higher overall clonality in secondary infection. Additionally, we observed consistent antibody sequence features in acute dengue in the highly variable major antigen-binding determinant, complementarity-determining region 3 (CDR3), with specific CDR3 sequences highly enriched in acute samples compared to postrecovery, healthy, or non-dengue samples. Dengue thus provides a striking example of a human viral infection where convergent immune signatures can be identified in multiple individuals. Such signatures could facilitate surveillance of immunological memory in communities.

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 IgG 1 -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.

Antibody lineages with evidence of somatic hypermutation persisting for >4 years in a South African subject with broad neutralizing activity M Moody, AM Trama, M Bonsignori, C Tsao, MS Drinker, TC Gurley, JD Amos, JA Eudailey, LC Armand, R Parks, KE Lloyd, S Wang, K Seo, J Lee, KJ Jackson, R Hoh, T Pham, KM Roskin, SD Boyd, AZ Fire, ES Gray, L Morris, H Liao, GD Tomaras, TB Kepler, G Kelsoe, BF Haynes

The mission of the Encyclopedia of DNA Elements (ENCODE) Project is to enable the scientific and medical communities to interpret the human genome sequence and apply it to understand human biology and improve health. The ENCODE Consortium is integrating multiple technologies and approaches in a collective effort to discover and define the functional elements encoded in the human genome, including genes, transcripts, and transcriptional regulatory regions, together with their attendant chromatin states and DNA methylation patterns. In the process, standards to ensure high-quality data have been implemented, and novel algorithms have been developed to facilitate analysis. Data and derived results are made available through a freely accessible database. Here we provide an overview of the project and the resources it is generating and illustrate the application of ENCODE data to interpret the human genome., National Human Genome Research Institute (U.S.), National Institutes of Health (U.S.)

Background Continuing research into the global multiple sequence alignment problem has resulted in more sophisticated and principled alignment methods. Unfortunately these new algorithms often require large amounts of time and memory to run, making it nearly impossible to run these algorithms on large datasets. As a solution, we present two general methods, Crumble and Prune, for breaking a phylogenetic alignment problem into smaller, more tractable sub-problems. We call Crumble and Prune meta-alignment methods because they use existing alignment algorithms and can be used with many current alignment programs. Crumble breaks long alignment problems into shorter sub-problems. Prune divides the phylogenetic tree into a collection of smaller trees to reduce the number of sequences in each alignment problem. These methods are orthogonal: they can be applied together to provide better scaling in terms of sequence length and in sequence depth. Both methods partition the problem such that many of the sub-problems can be solved independently. The results are then combined to form a solution to the full alignment problem. Results Crumble and Prune each provide a significant performance improvement with little loss of accuracy. In some cases, a gain in accuracy was observed. Crumble and Prune were tested on real and simulated data. Furthermore, we have implemented a system called Job-tree that allows hierarchical sub-problems to be solved in parallel on a compute cluster, significantly shortening the run-time. Conclusions These methods enabled us to solve gigabase alignment problems. These methods could enable a new generation of biologically realistic alignment algorithms to be applied to real world, large scale alignment problems.

We construct several score functions for use in locating unusually conserved regions in a genomewide search of aligned DNA from two species. We test these functions on regions of the human genome aligned to the mouse genome. These score functions are derived from properties of neutrally evolving sites on the mouse and human genome and can be adjusted to the local background rate of conservation. The aim of these functions is to try to identify regions of the human genome that are conserved by evolutionary selection because they have an important function, rather than by chance. We use them to get a very rough estimate of the amount of DNA in the human genome that is under selection.

Levels of recombination vary among species, among chromosomes within species, and among regions within chromosomes in mammals. This heterogeneity may affect levels of diversity, efficiency of selection, and genome composition, as well as have practical consequences for the genetic mapping of traits. We compared the genetic maps to the genome sequence assemblies of rat, mouse, and human to estimate local recombination rates across these genomes. Humans have greater overall levels of recombination, as well as greater variance. In rat and mouse, the size of the chromosome and proximity to telomere have less effect on local recombination rate than in human. At the chromosome level, rat and mouse X chromosomes have the lowest recombination rates, whereas human chromosome X does not show the same pattern. In all species, local recombination rate is significantly correlated with several sequence variables, including GC%, CpG density, repetitive elements, and the neutral mutation rate, with some pronounced differences between species. Recombination rate in one species is not strongly correlated with the rate in another, when comparing homologous syntenic blocks of the genome. This comparative approach provides additional insight into the causes and consequences of genomic heterogeneity in recombination.

The rates at which human genomic DNA changes by neutral substitution and insertion of certain families of transposable elements covary in large, megabase-sized segments. We used the rat, mouse, and human genomic DNA sequences to examine these processes in more detail in comparisons over both shorter (rat–mouse) and longer (rodent–primate) times, and demonstrated the generality of the covariation. Different families of transposable elements show distinctive insertion preferences and patterns of variation with substitution rates. SINEs are more abundant in GC-rich DNA, but the regional GC preference for insertion (monitored in young SINEs) differs between rodents and humans. In contrast, insertions in the rodent genomes are predominantly LINEs, which prefer to insert into AT-rich DNA in all three mammals. The insertion frequency of repeats other than SINEs correlates strongly positively with the frequency of substitutions in all species. However, correlations with SINEs show the opposite effects. The correlations are explained only in part by the GC content, indicating that other factors also contribute to the inherent tendency of DNA segments to change over evolutionary time.

The laboratory rat (Rattus norvegicus) is an indispensable tool in experimental medicine and drug development, having made inestimable contributions to human health. We report here the genome sequence of the Brown Norway (BN) rat strain. The sequence represents a high-quality 'draft' covering over 90% of the genome. The BN rat sequence is the third complete mammalian genome to be deciphered, and three-way comparisons with the human and mouse genomes resolve details of mammalian evolution. This first comprehensive analysis includes genes and proteins and their relation to human disease, repeated sequences, comparative genome-wide studies of mammalian orthologous chromosomal regions and rearrangement breakpoints, reconstruction of ancestral karyotypes and the events leading to existing species, rates of variation, and lineage-specific and lineage-independent evolutionary events such as expansion of gene families, orthology relations and protein evolution.

Six measures of evolutionary change in the human genome were studied, three derived from the aligned human and mouse genomes in conjunction with the Mouse Genome Sequencing Consortium, consisting of (1) nucleotide substitution per fourfold degenerate site in coding regions, (2) nucleotide substitution per site in relics of transposable elements active only before the human–mouse speciation, and (3) the nonaligning fraction of human DNA that is nonrepetitive or in ancestral repeats; and three derived from human genome data alone, consisting of (4) SNP density, (5) frequency of insertion of transposable elements, and (6) rate of recombination. Features 1 and 2 are measures of nucleotide substitutions at two classes of “neutral” sites, whereas 4 is a measure of recent mutations. Feature 3 is a measure dominated by deletions in mouse, whereas 5 represents insertions in human. It was found that all six vary significantly in megabase-sized regions genome-wide, and many vary together. This indicates that some regions of a genome change slowly by all processes that alter DNA, and others change faster. Regional variation in all processes is correlated with, but not completely accounted for, by GC content in human and the difference between GC content in human and mouse. [Supplemental material is available online atwww.genome.org andhttp://www.soe.ucsc.edu/research/compbio/covariation/.]

The sequence of the mouse genome is a key informational tool for understanding the contents of the human genome and a key experimental tool for biomedical research. Here, we report the results of an international collaboration to produce a high-quality draft sequence of the mouse genome. We also present an initial comparative analysis of the mouse and human genomes, describing some of the insights that can be gleaned from the two sequences. We discuss topics including the analysis of the evolutionary forces shaping the size, structure and sequence of the genomes; the conservation of large-scale synteny across most of the genomes; the much lower extent of sequence orthology covering less than half of the genomes; the proportions of the genomes under selection; the number of protein-coding genes; the expansion of gene families related to reproduction and immunity; the evolution of proteins; and the identification of intraspecies polymorphism.

As vertebrate genome sequences near completion and research refocuses to their analysis, the issue of effective genome annotation display becomes critical. A mature web tool for rapid and reliable display of any requested portion of the genome at any scale, together with several dozen aligned annotation tracks, is provided athttp://genome.ucsc.edu. This browser displays assembly contigs and gaps, mRNA and expressed sequence tag alignments, multiple gene predictions, cross-species homologies, single nucleotide polymorphisms, sequence-tagged sites, radiation hybrid data, transposon repeats, and more as a stack of coregistered tracks. Text and sequence-based searches provide quick and precise access to any region of specific interest. Secondary links from individual features lead to sequence details and supplementary off-site databases. One-half of the annotation tracks are computed at the University of California, Santa Cruz from publicly available sequence data; collaborators worldwide provide the rest. Users can stably add their own custom tracks to the browser for educational or research purposes. The conceptual and technical framework of the browser, its underlying MYSQL database, and overall use are described. The web site currently serves over 50,000 pages per day to over 3000 different users.

Abstract 2863 Background: JAK2 V617F is a recurrent, activating mutation in patients (pts) with BCR-ABL1-negative MPNs. Mutations in codon 515 of MPL occur in 1–5% and 5–10% of ET and PMF pts, respectively, and similar to JAK2 V617F, lead to constitutive JAK-STAT signaling. The acquisition of multiple mutational events affecting the JAK-STAT axis or components of the epigenetic machinery is common in MPNs and likely contributes to phenotypic diversity, including progression to acute myeloid leukemia. Mutated genes thus far implicated in MPN initiation and/or progression include TET2, CBL, SH2B3, ASXL1, DNMT3A, IDH1/2, IKZF1, EZH2, SRSF2, and TP53. In order to identify novel somatic mutations associated with classic BCR-ABL1-negative MPNs, we performed whole genome sequencing of DNA extracted from peripheral blood granulocytes and cultured skin fibroblasts from a patient with PMF and a known MPL W515K mutation. Methods: Whole genome sequencing (WGS) was undertaken in a 55 year-old man with untreated PMF four years after initial diagnosis. His DIPSS Plus risk group was low (score 0). His karyotype was normal, and molecular testing revealed wild-type JAK2 in addition to the MPL W515K mutation. WGS of purified granulocytes and paired cultured skin fibroblasts was performed using both Illumina HiSeq and Complete Genomics (CGI) platforms. The resulting data were analyzed using multiple independent aligners and variant callers. Amplicon-based targeted resequencing with the Illumina MiSeq platform was used to evaluate additional patient samples for recurrent mutations. Stanford institutional review board approval and informed pt consent was obtained for these analyses. Results: The PMF genome was sequenced to 88X (Illumina) and 128X (CGI) average fold coverage, and the cultured skin fibroblast genome was sequenced to 47X (Illumina) and 126X (CGI). The PMF genome had a low somatic mutation rate, consistent with that observed for other sequenced hematopoietic tumor genomes, with a low number of true somatic mutation calls. To definitively identify true mutations among various sequencing artefacts and germline variants, we use cultured skin fibroblasts which can be prepared with no contamination by neoplastic cells. In addition to re-identification of the MPL W515K mutation, this approach identified six additional somatic mutations that alter gene coding regions, splice sites, or known regulatory regions: a nonsense mutation in CARD6, implicated in modulation of NF-kappaB activation; a splice-site mutation in CAP2; three nonsynonymous point mutations in KIAA0355, SOX30, MFRP; and a 19-base pair (bp) deletion involving a regulatory region in the 5'-untranslated region (5'-UTR) of BRD2, a bromodomain-containing protein implicated in transcriptional regulation (Table). CARD6, BRD2, and KIAA0355, an uncharacterized protein, are expressed by the granulocytes derived from this patient, supporting a potential role in the development of PMF in this pt. Using massively parallel sequencing, we are currently examining the transcribed region of BRD2 and the coding region of the other five genes in additional pt samples. Analysis of the first 87 samples (MF=47, PV=20, ET=20) out of a cohort of 180 MPN pts has thus far not identified recurrent somatic mutations in these six genes. Conclusion: High-coverage genome sequencing of neoplastic and germline cells from a patient with MPL-mutated PMF identified six additional somatic mutations of potential functional significance. Work is ongoing to determine if somatic mutations in these genes are found in other patients with BCR-ABL1-negative MPNs, and their pathogenetic relevance to MPN biology. Disclosures: Snyder: Illumina: Consultancy; GenapSys: Membership on an entity's Board of Directors or advisory committees; Personalis: Consultancy, Founder Other.