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 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

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

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.

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.

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.

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

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.

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.

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

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.

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

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.

We present an algorithm for simplifying terrain data that preserves topology. We use a decimation algorithm that simplifies the given data set using hierarchical clustering. Topology constraints, along with local error metrics, are used to ensure topology-preserving simplification and to compute precise error bounds in the simplified data. The earth's mover distance is used as a global metric to compute the degradation in topology as the simplification proceeds. Experiments with both analytic and real terrain data are presented. Results indicate that one can obtain significant simplification with low errors without losing topology information.

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

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

We construct several score functions for use in locating unusually conserved regions in a genome-wide 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.

We describe three systems that use natural or event-based sounds as means of data delivery. In these systems we have mapped data to natural sounds using metaphors. In the first system we evaluate the use of sounds of air, horn, and train to convey ordered numeric values between 1 to 6. An example of the metaphor used here is the association of speed values to the sound of a moving train at different speeds. In the second system, we use sounds of ocean waves to convey whether the exposure in a protein structural alignment is buried, partially exposed or fully exposed. The metaphor used here is the association of sound with how exposed the user is with respect to the ocean. In the third system, we map animal sounds such as the sound of a roaring lion or a chirping bird to certain stocks based on user preferences. The behavior of the stocks are then sounded by the use of whistles and car crash to signify the movement in process of the stocks. An up whistling sound can be clearly associated with an uptrend. We present and discuss the results of user evaluation studies for all the three systems.

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.