Predictability of B cell clonal persistence and immunosurveillance in breast cancer.

B cells and T cells are important components of the adaptive immune system and mediate anticancer immunity. The T cell landscape in cancer is well characterized, but the contribution of B cells to anticancer immunosurveillance is less well explored. Here we show an integrative analysis of the B cell and T cell receptor repertoire from individuals with metastatic breast cancer and individuals with early breast cancer during neoadjuvant therapy. Using immune receptor, RNA and whole-exome sequencing, we show that both B cell and T cell responses seem to coevolve with the metastatic cancer genomes and mirror tumor mutational and neoantigen architecture. B cell clones associated with metastatic immunosurveillance and temporal persistence were more expanded and distinct from site-specific clones. B cell clonal immunosurveillance and temporal persistence are predictable from the clonal structure, with higher-centrality B cell antigen receptors more likely to be detected across multiple metastases or across time. This predictability was generalizable across other immune-mediated disorders. This work lays a foundation for prioritizing antibody sequences for therapeutic targeting in cancer.

Dynamic regulation of hypoxia-inducible factor-1α activity is essential for normal B cell development.

B lymphocyte development and selection are central to adaptive immunity and self-tolerance. These processes require B cell receptor (BCR) signaling and occur in bone marrow, an environment with variable hypoxia, but whether hypoxia-inducible factor (HIF) is involved is unknown. We show that HIF activity is high in human and murine bone marrow pro-B and pre-B cells and decreases at the immature B cell stage. This stage-specific HIF suppression is required for normal B cell development because genetic activation of HIF-1α in murine B cells led to reduced repertoire diversity, decreased BCR editing and developmental arrest of immature B cells, resulting in reduced peripheral B cell numbers. HIF-1α activation lowered surface BCR, CD19 and B cell-activating factor receptor and increased expression of proapoptotic BIM. BIM deletion rescued the developmental block. Administration of a HIF activator in clinical use markedly reduced bone marrow and transitional B cells, which has therapeutic implications. Together, our work demonstrates that dynamic regulation of HIF-1α is essential for normal B cell development.

Clonally expanded B cells in multiple sclerosis bind EBV EBNA1 and GlialCAM.

Multiple sclerosis (MS) is a heterogenous autoimmune disease in which autoreactive lymphocytes attack the myelin sheath of the central nervous system. B lymphocytes in the cerebrospinal fluid (CSF) of patients with MS contribute to inflammation and secrete oligoclonal immunoglobulins1,2. Epstein-Barr virus (EBV) infection has been epidemiologically linked to MS, but its pathological role remains unclear3. Here we demonstrate high-affinity molecular mimicry between the EBV transcription factor EBV nuclear antigen 1 (EBNA1) and the central nervous system protein glial cell adhesion molecule (GlialCAM) and provide structural and in vivo functional evidence for its relevance. A cross-reactive CSF-derived antibody was initially identified by single-cell sequencing of the paired-chain B cell repertoire of MS blood and CSF, followed by protein microarray-based testing of recombinantly expressed CSF-derived antibodies against MS-associated viruses. Sequence analysis, affinity measurements and the crystal structure of the EBNA1-peptide epitope in complex with the autoreactive Fab fragment enabled tracking of the development of the naive EBNA1-restricted antibody to a mature EBNA1-GlialCAM cross-reactive antibody. Molecular mimicry is facilitated by a post-translational modification of GlialCAM. EBNA1 immunization exacerbates disease in a mouse model of MS, and anti-EBNA1 and anti-GlialCAM antibodies are prevalent in patients with MS. Our results provide a mechanistic link for the association between MS and EBV and could guide the development of new MS therapies.

HIV/HBV coinfection remodels the immune landscape and natural killer cell ADCC functional responses.

BACKGROUND AND AIMS: HBV and HIV coinfection is a common occurrence globally, with significant morbidity and mortality. Both viruses lead to immune dysregulation including changes in natural killer (NK) cells, a key component of antiviral defense and a promising target for HBV cure strategies. Here we used high-throughput single-cell analysis to explore the immune cell landscape in people with HBV mono-infection and HIV/HBV coinfection, on antiviral therapy, with emphasis on identifying the distinctive characteristics of NK cell subsets that can be therapeutically harnessed. APPROACH AND RESULTS: Our data show striking differences in the transcriptional programs of NK cells. HIV/HBV coinfection was characterized by an over-representation of adaptive, KLRC2 -expressing NK cells, including a higher abundance of a chemokine-enriched ( CCL3/CCL4 ) adaptive cluster. The NK cell remodeling in HIV/HBV coinfection was reflected in enriched activation pathways, including CD3ζ phosphorylation and ZAP-70 translocation that can mediate stronger antibody-dependent cellular cytotoxicity responses and a bias toward chemokine/cytokine signaling. By contrast, HBV mono-infection imposed a stronger cytotoxic profile on NK cells and a more prominent signature of "exhaustion" with higher circulating levels of HBsAg. Phenotypic alterations in the NK cell pool in coinfection were consistent with increased "adaptiveness" and better capacity for antibody-dependent cellular cytotoxicity compared to HBV mono-infection. Overall, an adaptive NK cell signature correlated inversely with circulating levels of HBsAg and HBV-RNA in our cohort. CONCLUSIONS: This study provides new insights into the differential signature and functional profile of NK cells in HBV and HIV/HBV coinfection, highlighting pathways that can be manipulated to tailor NK cell-focused approaches to advance HBV cure strategies in different patient groups.

Predicting the risk of pancreatic cancer in adults with new-onset diabetes: development and internal-external validation of a clinical risk prediction model.

BACKGROUND: The National Institute for Health and Care Excellence (NICE) recommends that people aged 60+ years with newly diagnosed diabetes and weight loss undergo abdominal imaging to assess for pancreatic cancer. More nuanced stratification could lead to enrichment of these referral pathways. METHODS: Population-based cohort study of adults aged 30-85 years at type 2 diabetes diagnosis (2010-2021) using the QResearch primary care database in England linked to secondary care data, the national cancer registry and mortality registers. Clinical prediction models were developed to estimate risks of pancreatic cancer diagnosis within 2 years and evaluated using internal-external cross-validation. RESULTS: Seven hundred and sixty-seven of 253,766 individuals were diagnosed with pancreatic cancer within 2 years. Models included age, sex, BMI, prior venous thromboembolism, digoxin prescription, HbA1c, ALT, creatinine, haemoglobin, platelet count; and the presence of abdominal pain, weight loss, jaundice, heartburn, indigestion or nausea (previous 6 months). The Cox model had the highest discrimination (Harrell's C-index 0.802 (95% CI: 0.797-0.817)), the highest clinical utility, and was well calibrated. The model's highest 1% of predicted risks captured 12.51% of pancreatic cancer cases. NICE guidance had 3.95% sensitivity. DISCUSSION: A new prediction model could have clinical utility in identifying individuals with recent onset diabetes suitable for fast-track abdominal imaging.

Temporality of body mass index, blood tests, comorbidities and medication use as early markers for pancreatic ductal adenocarcinoma (PDAC): a nested case-control study.

OBJECTIVE: Prior studies identified clinical factors associated with increased risk of pancreatic ductal adenocarcinoma (PDAC). However, little is known regarding their time-varying nature, which could inform earlier diagnosis. This study assessed temporality of body mass index (BMI), blood-based markers, comorbidities and medication use with PDAC risk . DESIGN: We performed a population-based nested case-control study of 28 137 PDAC cases and 261 219 matched-controls in England. We described the associations of biomarkers with risk of PDAC using fractional polynomials and 5-year time trends using joinpoint regression. Associations with comorbidities and medication use were evaluated using conditional logistic regression. RESULTS: Risk of PDAC increased with raised HbA1c, liver markers, white blood cell and platelets, while following a U-shaped relationship for BMI and haemoglobin. Five-year trends showed biphasic BMI decrease and HbA1c increase prior to PDAC; early-gradual changes 2-3 years prior, followed by late-rapid changes 1-2 years prior. Liver markers and blood counts (white blood cell, platelets) showed monophasic rapid-increase approximately 1 year prior. Recent diagnosis of pancreatic cyst, pancreatitis, type 2 diabetes and initiation of certain glucose-lowering and acid-regulating therapies were associated with highest risk of PDAC. CONCLUSION: Risk of PDAC increased with raised HbA1c, liver markers, white blood cell and platelets, while followed a U-shaped relationship for BMI and haemoglobin. BMI and HbA1c derange biphasically approximately 3 years prior while liver markers and blood counts (white blood cell, platelets) derange monophasically approximately 1 year prior to PDAC. Profiling these in combination with their temporality could inform earlier PDAC diagnosis.

Analysis of T cell receptor clonotypes in tumor microenvironment identifies shared cancer-type-specific signatures.

Despite the conventional view that a truly random V(D)J recombination process should generate a highly diverse immune repertoire, emerging reports suggest that there is a certain bias toward the generation of shared/public immune receptor chains. These studies were performed in viral diseases where public T cell receptors (TCR) appear to confer better protective responses. Selective pressures generating common TCR clonotypes are currently not well understood, but it is believed that they confer a growth advantage. As very little is known about public TCR clonotypes in cancer, here we set out to determine the extent of shared TCR clonotypes in the intra-tumor microenvironments of virus- and non-virus-driven head and neck cancers using TCR sequencing. We report that tumor-infiltrating T cell clonotypes were indeed shared across individuals with the same cancer type, where the majority of shared sequences were specific to the cancer type (i.e., viral versus non-viral). These shared clonotypes were not particularly enriched in EBV-associated nasopharynx cancer but, in both cancers, exhibited distinct characteristics, namely shorter CDR3 lengths, restricted V- and J-gene usages, and also demonstrated convergent V(D)J recombination. Many of these shared TCRs were expressed in patients with a shared HLA background. Pattern recognition of CDR3 amino acid sequences revealed strong convergence to specific pattern motifs, and these motifs were uniquely found to each cancer type. This suggests that they may be enriched for specificity to common antigens found in the tumor microenvironment of different cancers. The identification of shared TCRs in infiltrating tumor T cells not only adds to our understanding of the tumor-adaptive immune recognition but could also serve as disease-specific biomarkers and guide the development of future immunotherapies.

Using de novo assembly to identify structural variation of eight complex immune system gene regions.

Driven by the necessity to survive environmental pathogens, the human immune system has evolved exceptional diversity and plasticity, to which several factors contribute including inheritable structural polymorphism of the underlying genes. Characterizing this variation is challenging due to the complexity of these loci, which contain extensive regions of paralogy, segmental duplication and high copy-number repeats, but recent progress in long-read sequencing and optical mapping techniques suggests this problem may now be tractable. Here we assess this by using long-read sequencing platforms from PacBio and Oxford Nanopore, supplemented with short-read sequencing and Bionano optical mapping, to sequence DNA extracted from CD14+ monocytes and peripheral blood mononuclear cells from a single European individual identified as HV31. We use this data to build a de novo assembly of eight genomic regions encoding four key components of the immune system, namely the human leukocyte antigen, immunoglobulins, T cell receptors, and killer-cell immunoglobulin-like receptors. Validation of our assembly using k-mer based and alignment approaches suggests that it has high accuracy, with estimated base-level error rates below 1 in 10 kb, although we identify a small number of remaining structural errors. We use the assembly to identify heterozygous and homozygous structural variation in comparison to GRCh38. Despite analyzing only a single individual, we find multiple large structural variants affecting core genes at all three immunoglobulin regions and at two of the three T cell receptor regions. Several of these variants are not accurately callable using current algorithms, implying that further methodological improvements are needed. Our results demonstrate that assessing haplotype variation in these regions is possible given sufficiently accurate long-read and associated data. Continued reductions in the cost of these technologies will enable application of these methods to larger samples and provide a broader catalogue of germline structural variation at these loci, an important step toward making these regions accessible to large-scale genetic association studies.

Research Recommendations Following the Discovery of Pain Sensitizing IgG Autoantibodies in Fibromyalgia Syndrome.

BACKGROUND: Fibromyalgia syndrome (FMS) is the most common chronic widespread pain condition in rheumatology. Until recently, no clear pathophysiological mechanism for fibromyalgia had been established, resulting in management challenges. Recent research has indicated that serum immunoglobulin Gs (IgGs) may play a role in FMS. We undertook a research prioritisation exercise to identify the most pertinent research approaches that may lead to clinically implementable outputs. METHODS: Research priority setting was conducted in five phases: situation analysis; design; expert group consultation; interim recommendations; consultation and revision. A dialogue model was used, and an international multi-stakeholder expert group was invited. Clinical, patient, industry, funder, and scientific expertise was represented throughout. Recommendation-consensus was determined via a voluntary closed eSurvey. Reporting guideline for priority setting of health research were employed to support implementation and maximise impact. RESULTS: Arising from the expert group consultation (n = 29 participants), 39 interim recommendations were defined. A response rate of 81.5% was achieved in the consensus survey. Six recommendations were identified as high priority- and 15 as medium level priority. The recommendations range from aspects of fibromyalgia features that should be considered in future autoantibody research, to specific immunological investigations, suggestions for trial design in FMS, and therapeutic interventions that should be assessed in trials. CONCLUSIONS: By applying the principles of strategic priority setting we directed research towards that which is implementable, thereby expediating the benefit to the FMS patient population. These recommendations are intended for patients, international professionals and grant-giving bodies concerned with research into causes and management of patients with fibromyalgia syndrome.

Distinctive binding properties of human monoclonal LGI1 autoantibodies determine pathogenic mechanisms.

Autoantibodies against leucine-rich glioma inactivated 1 (LGI1) are found in patients with limbic encephalitis and focal seizures. Here, we generate patient-derived monoclonal antibodies (mAbs) against LGI1. We explore their sequences and binding characteristics, plus their pathogenic potential using transfected HEK293T cells, rodent neuronal preparations, and behavioural and electrophysiological assessments in vivo after mAb injections into the rodent hippocampus. In live cell-based assays, LGI1 epitope recognition was examined with patient sera (n = 31), CSFs (n = 11), longitudinal serum samples (n = 15), and using mAbs (n = 14) generated from peripheral B cells of two patients. All sera and 9/11 CSFs bound both the leucine-rich repeat (LRR) and the epitempin repeat (EPTP) domains of LGI1, with stable ratios of LRR:EPTP antibody levels over time. By contrast, the mAbs derived from both patients recognized either the LRR or EPTP domain. mAbs against both domain specificities showed varied binding strengths, and marked genetic heterogeneity, with high mutation frequencies. LRR-specific mAbs recognized LGI1 docked to its interaction partners, ADAM22 and ADAM23, bound to rodent brain sections, and induced internalization of the LGI1-ADAM22/23 complex in both HEK293T cells and live hippocampal neurons. By contrast, few EPTP-specific mAbs bound to rodent brain sections or ADAM22/23-docked LGI1, but all inhibited the docking of LGI1 to ADAM22/23. After intrahippocampal injection, and by contrast to the LRR-directed mAbs, the EPTP-directed mAbs showed far less avid binding to brain tissue and were consistently detected in the serum. Post-injection, both domain-specific mAbs abrogated long-term potentiation induction, and LRR-directed antibodies with higher binding strengths induced memory impairment. Taken together, two largely dichotomous populations of LGI1 mAbs with distinct domain binding characteristics exist in the affinity matured peripheral autoantigen-specific memory pools of individuals, both of which have pathogenic potential. In human autoantibody-mediated diseases, the detailed characterization of patient mAbs provides a valuable method to dissect the molecular mechanisms within polyclonal populations.

Pre/pro-B cells generate macrophage populations during homeostasis and inflammation.

Most tissue-resident macrophages (Mφs) are believed to be derived prenatally and are assumed to maintain themselves throughout life by self-proliferation. However, in adult mice we identified a progenitor within bone marrow, early pro-B cell/fraction B, that differentiates into tissue Mφs. These Mφ precursors have non-rearranged B-cell receptor genes and coexpress myeloid (GR1, CD11b, and CD16/32) and lymphoid (B220 and CD19) lineage markers. During steady state, these precursors exit bone marrow, losing Gr1, and enter the systemic circulation, seeding the gastrointestinal system as well as pleural and peritoneal cavities but not the brain. While in these tissues, they acquire a transcriptome identical to embryonically derived tissue-resident Mφs. Similarly, these Mφ precursors also enter sites of inflammation, gaining CD115, F4/80, and CD16/32, and become indistinguishable from blood monocyte-derived Mφs. Thus, we have identified a population of cells within the bone marrow early pro-B cell compartment that possess functional plasticity to differentiate into either tissue-resident or inflammatory Mφs, depending on microenvironmental signals. We propose that these precursors represent an additional source of Mφ populations in adult mice during steady state and inflammation.

Unraveling the cellular origin and clinical prognostic markers of infant B-cell acute lymphoblastic leukemia using genome-wide analysis.

B-cell acute lymphoblastic leukemia is the commonest childhood cancer. In infants, B-cell acute lymphoblastic leukemia remains fatal, especially in patients with t(4;11), present in ~80% of cases. The pathogenesis of t(4;11)/KMT2A-AFF1+ (MLL-AF4+) infant B-cell acute lymphoblastic leukemia remains difficult to model, and the pathogenic contribution in cancer of the reciprocal fusions resulting from derivative translocated-chromosomes remains obscure. Here, "multi-layered" genome-wide analyses and validation were performed on a total of 124 de novo cases of infant B-cell acute lymphoblastic leukemia uniformly diagnosed and treated according to the Interfant 99/06 protocol. These patients showed the most silent mutational landscape reported so far for any sequenced pediatric cancer. Recurrent mutations were exclusively found in K-RAS and N-RAS, were subclonal and were frequently lost at relapse, despite a larger number of non-recurrent/non-silent mutations. Unlike non-MLL-rearranged B-cell acute lymphoblastic leukemias, B-cell receptor repertoire analysis revealed minor, non-expanded B-cell clones in t(4;11)+ infant B-cell acute lymphoblastic leukemia, and RNA-sequencing showed transcriptomic similarities between t(4;11)+ infant B-cell acute lymphoblastic leukemias and the most immature human fetal liver hematopoietic stem and progenitor cells, confirming a "pre-VDJ" fetal cellular origin for both t(4;11) and RAS mut The reciprocal fusion AF4-MLL was expressed in only 45% (19/43) of the t(4;11)+ patients, and HOXA cluster genes are exclusively expressed in AF4-MLL-expressing patients. Importantly, AF4-MLL/HOXA-expressing patients had a significantly better 4-year event-free survival (62.4% vs 11.7%, P=0.001), and overall survival (73.7 vs 25.2%, P=0.016). AF4-MLL expression retained its prognostic significance when analyzed in a Cox model adjusting for risk stratification according to the Interfant-06 protocol based on age at diagnosis, white blood cell count and response to prednisone. This study has clinical implications for disease outcome and diagnostic risk-stratification of t(4;11)+ infant B-cell acute lymphoblastic leukemia.

Antibody repertoire analysis in polygenic autoimmune diseases.

High-throughput sequencing of the DNA/RNA encoding antibody heavy- and light-chains is rapidly transforming the field of adaptive immunity. It can address key questions, including: (i) how the B-cell repertoire differs in health and disease; and (ii) if it does differ, the point(s) in B-cell development at which this occurs. The advent of technologies, such as whole-genome sequencing, offers the chance to link abnormalities in the B-cell antibody repertoire to specific genomic variants and polymorphisms. Here, we discuss the current research using B-cell antibody repertoire sequencing in three polygenic autoimmune diseases where there is good evidence for a pathological role for B-cells, namely systemic lupus erythematosus, multiple sclerosis and rheumatoid arthritis. These autoimmune diseases exhibit significantly skewed B-cell receptor repertoires compared with healthy controls. Interestingly, some common repertoire defects are shared between diseases, such as elevated IGHV4-34 gene usage. B-cell clones have effectively been characterized and tracked between different tissues and blood in autoimmune disease. It has been hypothesized that these differences may signify differences in B-cell tolerance; however, the mechanisms and implications of these defects are not clear.

Network properties derived from deep sequencing of human B-cell receptor repertoires delineate B-cell populations.

The adaptive immune response selectively expands B- and T-cell clones following antigen recognition by B- and T-cell receptors (BCR and TCR), respectively. Next-generation sequencing is a powerful tool for dissecting the BCR and TCR populations at high resolution, but robust computational analyses are required to interpret such sequencing. Here, we develop a novel computational approach for BCR repertoire analysis using established next-generation sequencing methods coupled with network construction and population analysis. BCR sequences organize into networks based on sequence diversity, with differences in network connectivity clearly distinguishing between diverse repertoires of healthy individuals and clonally expanded repertoires from individuals with chronic lymphocytic leukemia (CLL) and other clonal blood disorders. Network population measures defined by the Gini Index and cluster sizes quantify the BCR clonality status and are robust to sampling and sequencing depths. BCR network analysis therefore allows the direct and quantifiable comparison of BCR repertoires between samples and intra-individual population changes between temporal or spatially separated samples and over the course of therapy.

Molecular evolution of broadly neutralizing Llama antibodies to the CD4-binding site of HIV-1.

To date, no immunization of humans or animals has elicited broadly neutralizing sera able to prevent HIV-1 transmission; however, elicitation of broad and potent heavy chain only antibodies (HCAb) has previously been reported in llamas. In this study, the anti-HIV immune responses in immunized llamas were studied via deep sequencing analysis using broadly neutralizing monoclonal HCAbs as a guides. Distinct neutralizing antibody lineages were identified in each animal, including two defined by novel antibodies (as variable regions called VHH) identified by robotic screening of over 6000 clones. The combined application of five VHH against viruses from clades A, B, C and CRF_AG resulted in neutralization as potent as any of the VHH individually and a predicted 100% coverage with a median IC50 of 0.17 µg/ml for the panel of 60 viruses tested. Molecular analysis of the VHH repertoires of two sets of immunized animals showed that each neutralizing lineage was only observed following immunization, demonstrating that they were elicited de novo. Our results show that immunization can induce potent and broadly neutralizing antibodies in llamas with features similar to human antibodies and provide a framework to analyze the effectiveness of immunization protocols.

Capturing needles in haystacks: a comparison of B-cell receptor sequencing methods.

BACKGROUND: Deep-sequencing methods are rapidly developing in the field of B-cell receptor (BCR) and T-cell receptor (TCR) diversity. These promise to revolutionise our understanding of adaptive immune dynamics, identify novel antibodies, and allow monitoring of minimal residual disease. However, different methods for BCR and TCR enrichment and amplification have been proposed. Here we perform the first systematic comparison between different methods of enrichment, amplification and sequencing for generating BCR and TCR repertoires using large sample numbers. RESULTS: Resampling from the same RNA or cDNA pool results in highly correlated and reproducible repertoires, but resampling low frequency clones leads to stochastic variance. Repertoires generated by different sequencing methods (454 Roche and Illumina MiSeq) and amplification methods (multiplex PCR, 5' Rapid amplification of cDNA ends (5'RACE), and RNA-capture) are highly correlated, and resulting IgHV gene frequencies between the different methods were not significantly different. Read length has an impact on captured repertoire structure, and ultimately full-length BCR sequences are most informative for repertoire analysis as diversity outside of the CDR is very useful for phylogenetic analysis. Additionally, we show RNA-based BCR repertoires are more informative than using DNA. CONCLUSIONS: Repertoires generated by different sequencing and amplification methods are consistent, but we show that read lengths, depths and error profiles should be considered in experimental design, and multiple sampling approaches could be employed to minimise stochastic sampling variation. This detailed investigation of immune repertoire sequencing methods is essential for informing basic and clinical research.