HMCES protects immunoglobulin genes specifically from deletions during somatic hypermutation.

Somatic hypermutation (SHM) produces point mutations in immunoglobulin (Ig) genes in B cells when uracils created by the activation-induced deaminase are processed in a mutagenic manner by enzymes of the base excision repair (BER) and mismatch repair (MMR) pathways. Such uracil processing creates DNA strand breaks and is susceptible to the generation of deleterious deletions. Here, we demonstrate that the DNA repair factor HMCES strongly suppresses deletions without significantly affecting other parameters of SHM in mouse and human B cells, thereby facilitating the production of antigen-specific antibodies. The deletion-prone repair pathway suppressed by HMCES operates downstream from the uracil glycosylase UNG and is mediated by the combined action of BER factor APE2 and MMR factors MSH2, MSH6, and EXO1. HMCES's ability to shield against deletions during SHM requires its capacity to form covalent cross-links with abasic sites, in sharp contrast to its DNA end-joining role in class switch recombination but analogous to its genome-stabilizing role during DNA replication. Our findings lead to a novel model for the protection of Ig gene integrity during SHM in which abasic site cross-linking by HMCES intercedes at a critical juncture during processing of vulnerable gapped DNA intermediates by BER and MMR enzymes.

Essential role for the transcription factor Bhlhe41 in regulating the development, self-renewal and BCR repertoire of B-1a cells.

Innate-like B-1a cells provide a first line of defense against pathogens, yet little is known about their transcriptional control. Here we identified an essential role for the transcription factor Bhlhe41, with a lesser contribution by Bhlhe40, in controlling B-1a cell differentiation. Bhlhe41-/-Bhlhe40-/- B-1a cells were present at much lower abundance than were their wild-type counterparts. Mutant B-1a cells exhibited an abnormal cell-surface phenotype and altered B cell receptor (BCR) repertoire exemplified by loss of the phosphatidylcholine-specific VH12Vκ4 BCR. Expression of a pre-rearranged VH12Vκ4 BCR failed to 'rescue' the mutant phenotype and revealed enhanced proliferation accompanied by increased cell death. Bhlhe41 directly repressed the expression of cell-cycle regulators and inhibitors of BCR signaling while enabling pro-survival cytokine signaling. Thus, Bhlhe41 controls the development, BCR repertoire and self-renewal of B-1a cells.

Polyreactive Broadly Neutralizing B cells Are Selected to Provide Defense against Pandemic Threat Influenza Viruses.

Polyreactivity is the ability of a single antibody to bind to multiple molecularly distinct antigens and is a common feature of antibodies induced upon pathogen exposure. However, little is known about the role of polyreactivity during anti-influenza virus antibody responses. By analyzing more than 500 monoclonal antibodies (mAbs) derived from B cells induced by numerous influenza virus vaccines and infections, we found mAbs targeting conserved neutralizing influenza virus hemagglutinin epitopes were polyreactive. Polyreactive mAbs were preferentially induced by novel viral exposures due to their broad viral binding breadth. Polyreactivity augmented mAb viral binding strength by increasing antibody flexibility, allowing for adaption to imperfectly conserved epitopes. Lastly, we found affinity-matured polyreactive B cells were typically derived from germline polyreactive B cells that were preferentially selected to participate in B cell responses over time. Together, our data reveal that polyreactivity is a beneficial feature of antibodies targeting conserved epitopes.

DNA flexibility can shape the preferential hypermutation of antibody genes.

Antibody-coding genes accumulate somatic mutations to achieve antibody affinity maturation. Genetic dissection using various mouse models has shown that intrinsic hypermutations occur preferentially and are predisposed in the DNA region encoding antigen-contacting residues. The molecular basis of nonrandom/preferential mutations is a long-sought question in the field. Here, we summarize recent findings on how single-strand (ss)DNA flexibility facilitates activation-induced cytidine deaminase (AID) activity and fine-tunes the mutation rates at a mesoscale within the antibody variable domain exon. We propose that antibody coding sequences are selected based on mutability during the evolution of adaptive immunity and that DNA mechanics play a noncoding role in the genome. The mechanics code may also determine other cellular DNA metabolism processes, which awaits future investigation.

Naive B Cells with High-Avidity Germline-Encoded Antigen Receptors Produce Persistent IgM+ and Transient IgG+ Memory B Cells.

Although immune memory often lasts for life, this is not the case for certain vaccines in some individuals. We sought a mechanism for this phenomenon by studying B cell responses to phycoerythrin (PE). PE immunization of mouse strains with Ighb immunoglobulin (Ig) variable heavy chain (VH) genes elicited affinity-matured switched Ig memory B cells that declined with time, while the comparable population from an Igha strain was numerically stable. Ighb strains had larger numbers of PE-specific naive B cells and generated smaller germinal center responses and larger numbers of IgM memory cells than the Igha strain. The properties of PE-specific B cells in Ighb mice correlated with usage of a single VH that afforded high-affinity PE binding in its germline form. These results suggest that some individuals may be genetically predisposed to generate non-canonical memory B cell responses to certain antigens because of avid antigen binding via germline-encoded VH elements.

The evolutionary and functional significance of germline immunoglobulin gene variation.

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Addressing Technical Pitfalls in Pursuit of Molecular Factors That Mediate Immunoglobulin Gene Regulation.

The expressed Ab repertoire is a critical determinant of immune-related phenotypes. Ab-encoding transcripts are distinct from other expressed genes because they are transcribed from somatically rearranged gene segments. Human Abs are composed of two identical H and L chain polypeptides derived from genes in IGH locus and one of two L chain loci. The combinatorial diversity that results from Ab gene rearrangement and the pairing of different H and L chains contributes to the immense diversity of the baseline Ab repertoire. During rearrangement, Ab gene selection is mediated by factors that influence chromatin architecture, promoter/enhancer activity, and V(D)J recombination. Interindividual variation in the composition of the Ab repertoire associates with germline variation in IGH, implicating polymorphism in Ab gene regulation. Determining how IGH variants directly mediate gene regulation will require integration of these variants with other functional genomic datasets. In this study, we argue that standard approaches using short reads have limited utility for characterizing regulatory regions in IGH at haplotype resolution. Using simulated and chromatin immunoprecipitation sequencing reads, we define features of IGH that limit use of short reads and a single reference genome, namely 1) the highly duplicated nature of the DNA sequence in IGH and 2) structural polymorphisms that are frequent in the population. We demonstrate that personalized diploid references enhance performance of short-read data for characterizing mappable portions of the locus, while also showing that long-read profiling tools will ultimately be needed to fully resolve functional impacts of IGH germline variation on expressed Ab repertoires.

Gene prediction in the immunoglobulin loci.

The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin (IG) loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) are critical for studies of the adaptive immune system. However, because the standard gene finding algorithms are not suitable for predicting IG genes, they have been semimanually annotated in very few species. We developed the IGDetective algorithm for predicting IG genes and applied it to species with the assembled IG loci. IGDetective generated the first large collection of IG genes across many species and enabled their evolutionary analysis, including the analysis of the "bat IG diversity" hypothesis. This analysis revealed extremely conserved V genes in evolutionary distant species, indicating that these genes may be subjected to the same selective pressure, for example, pressure driven by common pathogens. IGDetective also revealed extremely diverged V genes and a new family of evolutionary conserved V genes in bats with unusual noncanonical cysteines. Moreover, unlike all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Because cysteines form disulfide bonds, we hypothesize that these cysteine-rich V genes might generate antibodies with noncanonical conformations and could potentially form a unique part of the immune repertoire in bats. We also analyzed the diversity landscape of the recombination signal sequences and revealed their features that trigger the high/low usage of the IG genes.

Mutagenic repair during antibody diversification: emerging insights.

Deoxyuracils (dUs) produced by activation-induced cytidine deaminase (AID) during antibody diversification are processed by base excision repair (BER) and mismatch repair (MMR) pathways that paradoxically expand this lesion within jawed vertebrate immunoglobulin (Ig) genes. We highlight new findings describing mechanisms that allow B cells to carry out mutagenic DNA repair, an essential process for antibody maturation with implications in cancer pathogenesis.

Diversity of B Cell Populations and Ig Repertoire in Human Lungs.

The human lung carries a unique microbiome adapted to the air-filled, mucous-lined environment, the presence of which requires an immune system capable of recognizing harmful populations while preventing reactions toward commensals. B cells in the lung play a key role in pulmonary immunity, generating Ag-specific Abs, as well as cytokine secretion for immune activation and regulation. In this study, we compared B cell subsets in human lungs versus circulating cells by analyzing patient-paired lung and blood samples. We found a significantly smaller pool of CD19+, CD20+ B cells in the lung relative to the blood. CD27+, IgD-, class-switched memory B cells (Bmems) composed a larger proportion of the pool of pulmonary B cells. The residency marker CD69 was also significantly higher in the lung. We also sequenced the Ig V region genes (IgVRGs) of class-switched Bmems that do, or do not, express CD69. We observed the IgVRGs of pulmonary Bmems to be as heavily mutated from the unmutated common ancestor as those in circulation. Furthermore, we found progenies within a quasi-clone can gain or lose CD69 expression, regardless of whether the parent clone expressed the residency marker. Overall, our results show that despite its vascularized nature, human lungs carry a unique proportion of B cell subsets. The IgVRGs of pulmonary Bmems are as diverse as those in blood, and progenies of Bmems retain the ability to gain or lose residency.

Cell Cycle-Mediated Regulation of Secondary Ig Diversification.

Secondary Ig diversification in B cells requires the deliberate introduction of DNA damage into the Ig genes by the enzyme activation-induced cytidine deaminase (AID) and the error-prone resolution of AID-induced lesions. These processes must be tightly regulated because they may lead to lymphomagenesis if they act on genes other than the Ig genes. Since B cells may limit secondary Ig diversification mechanisms during the cell cycle to minimize genomic instability, we restricted the activity of AID specifically to the G1 or S/G2 phase to investigate the cell cycle contribution to the regulation of somatic hypermutation, class switch recombination, and Ig gene conversion in human, murine, and avian B cells, respectively. The efficient induction of AID in different cell cycle phases allowed us for the first time, to our knowledge, to discriminate G1- from S/G2-specific events of regulation. We show that the processes of Ig gene conversion and C/G mutagenesis during somatic hypermutation can be achieved throughout the cell cycle, whereas A/T mutagenesis and class switch recombination require AID-mediated deamination in G1. Thus, AID activity in G1, but not in S/G2, leads to the efficient accomplishment of all mechanisms of secondary Ig diversification. Our findings refine the current state-of-the-art knowledge in the context of the regulation of secondary Ig diversification.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities.

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Different classes of genomic inserts contribute to human antibody diversity.

Recombination of antibody genes in B cells can involve distant genomic loci and contribute a foreign antigen-binding element to form hybrid antibodies with broad reactivity for Plasmodium falciparum. So far, antibodies containing the extracellular domain of the LAIR1 and LILRB1 receptors represent unique examples of cross-chromosomal antibody diversification. Here, we devise a technique to profile non-VDJ elements from distant genes in antibody transcripts. Independent of the preexposure of donors to malaria parasites, non-VDJ inserts were detected in 80% of individuals at frequencies of 1 in 104 to 105 B cells. We detected insertions in heavy, but not in light chain or T cell receptor transcripts. We classify the insertions into four types depending on the insert origin and destination: 1) mitochondrial and 2) nuclear DNA inserts integrated at VDJ junctions; 3) inserts originating from telomere proximal genes; and 4) fragile sites incorporated between J-to-constant junctions. The latter class of inserts was exclusively found in memory and in in vitro activated B cells, while all other classes were already detected in naïve B cells. More than 10% of inserts preserved the reading frame, including transcripts with signs of antigen-driven affinity maturation. Collectively, our study unravels a mechanism of antibody diversification that is layered on the classical V(D)J and switch recombination.

Systematic characterization of immunoglobulin loci and deep sequencing of the expressed repertoire in the Atlantic cod (Gadus morhua).

BACKGROUND: The Atlantic cod is a prolific species in the Atlantic, despite its inconsistent specific antibody response. It presents a peculiar case within vertebrate immunology due to its distinct immune system, characterized by the absence of MHCII antigen presentation pathway, required for T cell-dependent antibody responses. Thorough characterisation of immunoglobulin loci and analysis of the antibody repertoire is necessary to further our understanding of the Atlantic cod's immune response on a molecular level. RESULTS: A comprehensive search of the cod genome (gadmor3.0) identified the complete set of IgH genes organized into three sequential translocons on chromosome 2, while IgL genes were located on chromosomes 2 and 5. The Atlantic cod displayed a moderate germline V gene diversity, comprising four V gene families for both IgH and IgL, each with distinct chromosomal locations and organizational structures. 5'RACE sequencing revealed a diverse range of heavy chain CDR3 sequences and relatively limited CDR3 diversity in light chains. The analysis highlighted a differential impact of V-gene germline CDR3 length on receptor CDR3 length between heavy and light chains, underlining different recombination processes. CONCLUSIONS: This study reveals that the Atlantic cod, despite its inconsistent antibody response, maintains a level of immunoglobulin diversity comparable to other fish species. The findings suggest that the extensive recent duplications of kappa light chain genes do not result in increased repertoire diversity. This research provides a comprehensive view of the Atlantic cod's immunoglobulin gene organization and repertoire, necessary for future studies of antibody responses at the molecular level.

Immunoglobulin genes and severity of COVID-19.

There is tremendous interindividual and interracial variability in the outcome of SARS-CoV-2 infection, suggesting the involvement of host genetic factors. Here, we investigated whether IgG allotypes GM (γ marker) 3 and GM 17, genetic markers of IgG1, contributed to the severity of COVID-19. IgG1 plays a pivotal role in response against SARS-CoV-2 infection. We also investigated whether these GM alleles synergistically/epistatically with IGHG3 and FCGR2A alleles-which have been previously implicated in COVID-19-modulated the extent of COVID-19 severity. The study population consisted of 316 COVID-19 patients who needed treatment in the intensive care unit of Hospital Universitario Central de Asturias. All individuals were genotyped for GM 3/17, IGHG3 hinge length, and FCGR2A rs1801274 A/G polymorphisms. Among the 316 critical patients, there were 86 deaths. The risk of death among critical patients was significantly higher in subjects with GM 17 (IgG1) and short hinge length (IgG3). GM 17-carriers were at almost three-fold higher risk of death than non-carriers (p < 0.001; OR = 2.86, CI 1.58-5.16). Subjects with short hinge length of IgG3 had a two-fold higher risk of death than those with medium hinge length (p = 0.01; OR = 2.16, CI 1.19-3.90). GM 3/3 and IGHG3 (MM) genotypes were less frequent among death vs. survivors (9% vs 36%, p < 0.001) and associated with protective effect (OR = 0.18, 95% CI = 0.08-0.39). This is the first report implicating IgG1 allotypes in COVID-19-spurred death. It needs to be replicated in an independent study population.

IgIDivA: immunoglobulin intraclonal diversification analysis.

Intraclonal diversification (ID) within the immunoglobulin (IG) genes expressed by B cell clones arises due to ongoing somatic hypermutation (SHM) in a context of continuous interactions with antigen(s). Defining the nature and order of appearance of SHMs in the IG genes can assist in improved understanding of the ID process, shedding light into the ontogeny and evolution of B cell clones in health and disease. Such endeavor is empowered thanks to the introduction of high-throughput sequencing in the study of IG gene repertoires. However, few existing tools allow the identification, quantification and characterization of SHMs related to ID, all of which have limitations in their analysis, highlighting the need for developing a purpose-built tool for the comprehensive analysis of the ID process. In this work, we present the immunoglobulin intraclonal diversification analysis (IgIDivA) tool, a novel methodology for the in-depth qualitative and quantitative analysis of the ID process from high-throughput sequencing data. IgIDivA identifies and characterizes SHMs that occur within the variable domain of the rearranged IG genes and studies in detail the connections between identified SHMs, establishing mutational pathways. Moreover, it combines established and new graph-based metrics for the objective determination of ID level, combined with statistical analysis for the comparison of ID level features for different groups of samples. Of importance, IgIDivA also provides detailed visualizations of ID through the generation of purpose-built graph networks. Beyond the method design, IgIDivA has been also implemented as an R Shiny web application. IgIDivA is freely available at https://bio.tools/igidiva.

Differences in transcriptional changes in psoriasis and psoriatic arthritis skin with immunoglobulin gene enrichment in psoriatic arthritis.

OBJECTIVES: Approximately 20% of people with psoriasis develop PsA. Although genetic, clinical and environmental risk factors have been identified, it is not known why some people with psoriasis develop PsA. The skin disease is traditionally considered the same in both. This study compares transcriptional changes in psoriasis and PsA skin for the first time. METHODS: Skin biopsies were collected from healthy controls (HC), and uninvolved and lesional skin from patients with PsA. Bulk tissue sequencing was performed and analysed using the pipeline Searchlight 2.0. Transcriptional changes in PsA skin were compared with existing sequencing data from participants with psoriasis without PsA (GSE121212). Psoriasis and PsA datasets could not be directly compared as different analysis methods were used. Data from participants with PsA in the GSE121212 dataset were used for validation. RESULTS: Skin samples from 9 participants with PsA and 9 HC were sequenced, analysed and compared with available transcriptomic data for 16 participants with psoriasis compared with 16 HC. Uninvolved skin in psoriasis shared transcriptional changes with lesional skin in psoriasis, but uninvolved skin in PsA did not. Most transcriptional changes in psoriasis and PsA lesional skin were shared, but immunoglobulin genes were upregulated in PsA lesional skin specifically. The transcription factor POU2F1, which regulates immunoglobulin gene expression, was enriched in PsA lesional skin. This was confirmed in the validation cohort. CONCLUSIONS: Immunoglobulin genes are upregulated in PsA but not in psoriasis skin lesions. This may have implications for the spread from the cutaneous compartment to other tissues.

Immunoglobulin genes expressed in lymphoblastoid cell lines discern and predict lithium response in bipolar disorder patients.

Bipolar disorder (BD) is a neuropsychiatric mood disorder manifested by recurrent episodes of mania and depression. More than half of BD patients are non-responsive to lithium, the first-line treatment drug, complicating BD clinical management. Given its unknown etiology, it is pertinent to understand the genetic signatures that lead to variability in lithium response. We discovered a set of differentially expressed genes (DEGs) from the lymphoblastoid cell lines (LCLs) of 10 controls and 19 BD patients belonging mainly to the immunoglobulin gene family that can be used as potential biomarkers to diagnose and treat BD. Importantly, we trained machine learning algorithms on our datasets that predicted the lithium response of BD subtypes with minimal errors, even when used on a different cohort of 24 BD patients acquired by a different laboratory. This proves the scalability of our methodology for predicting lithium response in BD and for a prompt and suitable decision on therapeutic interventions.

Activation-Induced Cytidine Deaminase Expression in Human B Cell Precursors Is Essential for Central B Cell Tolerance.

Activation-induced cytidine deaminase (AID), the enzyme-mediating class-switch recombination (CSR) and somatic hypermutation (SHM) of immunoglobulin genes, is essential for the removal of developing autoreactive B cells. How AID mediates central B cell tolerance remains unknown. We report that AID enzymes were produced in a discrete population of immature B cells that expressed recombination-activating gene 2 (RAG2), suggesting that they undergo secondary recombination to edit autoreactive antibodies. However, most AID+ immature B cells lacked anti-apoptotic MCL-1 and were deleted by apoptosis. AID inhibition using lentiviral-encoded short hairpin (sh)RNA in B cells developing in humanized mice resulted in a failure to remove autoreactive clones. Hence, B cell intrinsic AID expression mediates central B cell tolerance potentially through its RAG-coupled genotoxic activity in self-reactive immature B cells.

Somatic Diversification of Rearranged Antibody Gene Segments by Intra- and Interchromosomal Templated Mutagenesis.

The ability of the humoral immune system to generate Abs capable of specifically binding a myriad of Ags is critically dependent on the somatic hypermutation program. This program induces both templated mutations (i.e., gene conversion) and untemplated mutations. In humans, somatic hypermutation is widely believed to result in untemplated point mutations. In this study, we demonstrate detection of large-scale templated events that occur in human memory B cells and circulating plasmablasts. We find that such mutations are templated intrachromosomally from IGHV genes and interchromosomally from IGHV pseudogenes as well as other homologous regions unrelated to IGHV genes. These same donor regions are used in multiple individuals, and they predominantly originate from chromosomes 14, 15, and 16. In addition, we find that exogenous sequences placed at the IgH locus, such as LAIR1, undergo templated mutagenesis and that homology appears to be the major determinant for donor choice. Furthermore, we find that donor tracts originate from areas in proximity with open chromatin, which are transcriptionally active, and are found in spatial proximity with the IgH locus during the germinal center reaction. These donor sequences are inserted into the Ig gene segment in association with overlapping activation-induced cytidine deaminase hotspots. Taken together, these studies suggest that diversity generated during the germinal center response is driven by untemplated point mutations as well as templated mutagenesis using local and distant regions of the genome.