Epstein-Barr virus perpetuates B cell germinal center dynamics and generation of autoimmune-associated phenotypes in vitro.

Human B cells encompass functionally diverse lineages and phenotypic states that contribute to protective as well as pathogenic responses. Epstein-Barr virus (EBV) provides a unique lens for studying heterogeneous B cell responses, given its adaptation to manipulate intrinsic cell programming. EBV promotes the activation, proliferation, and eventual outgrowth of host B cells as immortalized lymphoblastoid cell lines (LCLs) in vitro, which provide a foundational model of viral latency and lymphomagenesis. Although cellular responses and outcomes of infection can vary significantly within populations, investigations that capture genome-wide perspectives of this variation at single-cell resolution are in nascent stages. We have recently used single-cell approaches to identify EBV-mediated B cell heterogeneity in de novo infection and within LCLs, underscoring the dynamic and complex qualities of latent infection rather than a singular, static infection state. Here, we expand upon these findings with functional characterizations of EBV-induced dynamic phenotypes that mimic B cell immune responses. We found that distinct subpopulations isolated from LCLs could completely reconstitute the full phenotypic spectrum of their parental lines. In conjunction with conserved patterns of cell state diversity identified within scRNA-seq data, these data support a model in which EBV continuously drives recurrent B cell entry, progression through, and egress from the Germinal Center (GC) reaction. This "perpetual GC" also generates tangent cell fate trajectories including terminal plasmablast differentiation, which constitutes a replicative cul-de-sac for EBV from which lytic reactivation provides escape. Furthermore, we found that both established EBV latency and de novo infection support the development of cells with features of atypical memory B cells, which have been broadly associated with autoimmune disorders. Treatment of LCLs with TLR7 agonist or IL-21 was sufficient to generate an increased frequency of IgD-/CD27-/CD23-/CD38+/CD138+ plasmablasts. Separately, de novo EBV infection led to the development of CXCR3+/CD11c+/FCRL4+ B cells within days, providing evidence for possible T cell-independent origins of a recently described EBV-associated neuroinvasive CXCR3+ B cell subset in patients with multiple sclerosis. Collectively, this work reveals unexpected virus-driven complexity across infected cell populations and highlights potential roles of EBV in mediating or priming foundational aspects of virus-associated immune cell dysfunction in disease.

Establishment of mouse model of inherited PIGO deficiency and therapeutic potential of AAV-based gene therapy.

Inherited glycosylphosphatidylinositol (GPI) deficiency (IGD) is caused by mutations in GPI biosynthesis genes. The mechanisms of its systemic, especially neurological, symptoms are not clarified and fundamental therapy has not been established. Here, we report establishment of mouse models of IGD caused by PIGO mutations as well as development of effective gene therapy. As the clinical manifestations of IGD are systemic and lifelong lasting, we treated the mice with adeno-associated virus for homology-independent knock-in as well as extra-chromosomal expression of Pigo cDNA. Significant amelioration of neuronal phenotypes and growth defect was achieved, opening a new avenue for curing IGDs.

Phenotype and cytokine profile in a TRAPS Syndrome family with TNFRSF1A p.(Thr79Met): Association with sacro-iliitis.

OBJECTIVE: A comparative retrospective and analytic study was performed in 13 members of a family, affected or not by tumor necrosis factor TNF-α receptor-associated periodic syndrome (TRAPS), including one patient with sacro-illitis. METHODS: Clinical features and TNFRSF1A gene analysis were reported for each family member, symptomatic or asymptomatic. Biological features including CRP/SAA, IgD and ex vivo T lymphocytes and myeloid-derived cytokine profile (IL1-ß, IL-1α, IL-1ra, IL-4, IL-10, IL-17, IFN-γ, TNF-α, IL-6) were characterized for all family members. CRP and IgD only were compared with those of 13 axial SA patients with TNF-blockade indication. RESULTS: Symptomatic TRAPS patients presented p.(Thr79Met) variant and generally complained of abdominal pain. They displayed higher SAA/IgD levels and IL-6/IL-1RA/IL-10 production by PBMC compared to asymptomatic family members. IgD serum level was higher in symptomatic members compared to SA patients. The patient with sacro-iliitis displayed the highest cytokine production and IgD serum levels. CONCLUSION: This study describes clinical heterogeneity in a family examined for TRAPS syndrome and reports the first sacro-iliitis in a patient with pathogenic TNFRSF1A variant. Dysregulated PBMC-derived cytokine and Il-10/IgD dysregulation in the patient with sacro-iliitis raises the issue of sacro-iliitis pathophysiology.

[Analysis of MVK gene variant in a child with high IgD syndrome caused by mevalonate kinase deficiency].

OBJECTIVE: To analyze the clinical and genetic features of a patient with mevalonate kinase deficiency (MKD). METHODS: Whole exome sequencing was carried out for the proband. Candidate variant was verified by Sanger sequencing. RESULTS: The proband was found to harbor compound heterozygous variants of the MVK gene, including a c.248C>T (p.Phe83Cys) variant derived from his father and a c.971C>T (p.Ala324Val) variant from his mother. Based on the guidelines of the American College of Medical Genetics and Genomics, both variations were predicted to be likely pathogenic (PM1 + PM2 + PM3 + PP3). CONCLUSION: The compound heterozygous variants of the MVK gene probably underlay the MKD in the proband. Above findings have enriched the mutational spectrum of the MVK gene.

Rad52 mediates class-switch DNA recombination to IgD.

In B cells, IgD is expressed together with IgM through alternative splicing of primary VHDJH-Cµ-s-m-Cδ-s-m RNAs, and also through IgD class switch DNA recombination (CSR) via double-strand DNA breaks (DSB) and synapse of Sµ with σδ. How such DSBs are resolved is still unknown, despite our previous report showing that Rad52 effects the 'short-range' microhomology-mediated synapsis of intra-Sµ region DSBs. Here we find that induction of IgD CSR downregulates Zfp318, and promotes Rad52 phosphorylation and recruitment to Sµ and σδ, thereby leading to alternative end-joining (A-EJ)-mediated Sµ-σδ recombination with extensive microhomologies, VHDJH-Cδs transcription and sustained IgD secretion. Rad52 ablation in mouse Rad52-/- B cells aborts IgD CSR in vitro and in vivo and dampens the specific IgD antibody response to OVA. Rad52 knockdown in human B cells also abrogates IgD CSR. Finally, Rad52 phosphorylation is associated with high levels of IgD CSR and anti-nuclear IgD autoantibodies in patients with systemic lupus erythematosus and in lupus-prone mice. Our findings thus show that Rad52 mediates IgD CSR through microhomology-mediated A-EJ in concert with Zfp318 downregulation.

B cell receptor isotypes differentially associate with cell signaling, kinetics, and outcome in chronic lymphocytic leukemia.

In chronic lymphocytic leukemia (CLL), the B cell receptor (BCR) plays a critical role in disease development and progression, as indicated by the therapeutic efficacy of drugs blocking BCR signaling. However, the mechanism(s) underlying BCR responsiveness are not completely defined. Selective engagement of membrane IgM or IgD on CLL cells, each coexpressed by more than 90% of cases, leads to distinct signaling events. Since both IgM and IgD carry the same antigen-binding domains, the divergent actions of the receptors are attributed to differences in immunoglobulin (Ig) structure or the outcome of signal transduction. We showed that IgM, not IgD, level and organization associated with CLL-cell birth rate and the type and consequences of BCR signaling in humans and mice. The latter IgM-driven effects were abrogated when BCR signaling was inhibited. Collectively, these studies demonstrated a critical, selective role for IgM in BCR signaling and B cell fate decisions, possibly opening new avenues for CLL therapy.

Analysis of MVK gene variant in a child with high IgD syndrome caused by mevalonate kinase deficiency / 中华医学遗传学杂志

OBJECTIVE@#To analyze the clinical and genetic features of a patient with mevalonate kinase deficiency (MKD).@*METHODS@#Whole exome sequencing was carried out for the proband. Candidate variant was verified by Sanger sequencing.@*RESULTS@#The proband was found to harbor compound heterozygous variants of the MVK gene, including a c.248C>T (p.Phe83Cys) variant derived from his father and a c.971C>T (p.Ala324Val) variant from his mother. Based on the guidelines of the American College of Medical Genetics and Genomics, both variations were predicted to be likely pathogenic (PM1 + PM2 + PM3 + PP3).@*CONCLUSION@#The compound heterozygous variants of the MVK gene probably underlay the MKD in the proband. Above findings have enriched the mutational spectrum of the MVK gene.

IgD promotes pannus formation by activating Wnt5A-Fzd5-CTHRC1-NF-κB signaling pathway in FLS of CIA rats and the regulation of IgD-Fc-Ig fusion protein.

Rheumatoid arthritis (RA) is a systemic autoimmune disease characterized by joint inflammation, synovial hyperplasia, cartilage degeneration, bone erosion, and pannus. Immunoglobulin D (IgD) plays an important role in autoimmune diseases although the content of it in vivo is low. Increased concentrations of anti-IgD autoantibodies have been detected in many RA patients. IgD-Fc-Ig fusion protein is constructed by connecting human IgD Fc domain and IgG1 Fc domain, which specifically block the IgD/ IgDR pathway and regulate the function of cells expressing IgDR to treat RA. The expression levels of Wnt5A and Frizzled 5 are higher in RA synovial tissue specimens. The complex of Wnt5A-Fzd5-LRP5/6-CTHRC1 promotes the expression of hypoxia inducible factor-1α by activating nuclear factor kappa-B (NF-κB), leading to high expression of VEGF and participating in angiogenesis. VEGF is the strongest angiogenic factor found so far. Here, we aimed to explore whether IgD participates in synovitis by binding to IgDR and regulating the activation of Wnt5A-Fzd5-CTHRC1-NF-κB signaling pathway in fibroblast synovial cells (FLSs), whether IgD-Fc-Ig fusion protein inhibits VEGF production in FLS of CIA and explore mechanism. We found that IgDR is expressed on MH7A and FLS. IgD promotes VEGF expression by activating Wnt5A-Fzd5-CTHRC1-NF-κB signaling pathway in MH7A and FLS. After activation of Fzd5 with Wnt5A, IgD-Fc-Ig reduced VEGF-A level in the culture supernatant of MH7A stimulation by IgD. The expressions of CTHRC1, Fzd5, p-P65 and VEGF in MH7A and FLSs were down-regulated after IgD-Fc-Ig treatment. IgD-Fc-Ig suppressed the combination of CTHRC1 and Fzd5 as well. By using the animal model, we demonstrated that IgD-Fc-Ig suppress ankle CTHRC1 and Fzd5 production resulted in inhibition of index of joint inflammation of CIA rats, which were consistent with vitro results. Conclusively, IgD-Fc-Ig inhibits IgD and Wnt5A-induced angiogenesis and joint inflammation by suppressing the combination of CTHRC1 and Fzd5. Our results show that IgD-Fc-Ig exerts its suppressive effect on IgD and Wnt5A by Wnt5A-Fzd5-CTHRC1-NF-κB signaling pathway.

High-Throughput Detection of Autoantigen-Specific B Cells Among Distinct Functional Subsets in Autoimmune Donors.

Antigen-specific B cells (ASBCs) can drive autoimmune disease by presenting autoantigen to cognate T cells to drive their activation, proliferation, and effector cell differentiation and/or by differentiating into autoantibody-secreting cells. Autoantibodies are frequently used to predict risk and diagnose several autoimmune diseases. ASBCs can drive type 1 diabetes even when immune tolerance mechanisms block their differentiation into antibody-secreting cells. Furthermore, anti-histidyl tRNA synthetase syndrome patients have expanded IgM+ Jo-1-binding B cells, which clinically diagnostic IgG Jo-1 autoantibodies may not fully reflect. Given the potential disconnect between the pathologic function of ASBCs and autoantibody secretion, direct study of ASBCs is a necessary step towards developing better therapies for autoimmune diseases, which often have no available cure. We therefore developed a high-throughput screening pipeline to 1) phenotypically identify specific B cell subsets, 2) expand them in vitro, 3) drive them to secrete BCRs as antibody, and 4) identify wells enriched for ASBCs through ELISA detection of antibody. We tested the capacity of several B cell subset(s) to differentiate into antibody-secreting cells following this robust stimulation. IgM+ and/or IgD+, CD27- memory, memory, switched memory, and BND B cells secreted B cell receptor (BCR) as antibody following in vitro stimulation, whereas few plasmablasts responded. Bimodal responses were observed across autoimmune donors for IgM+ CD21lo and IgM- CD21lo B cells, consistent with documented heterogeneity within the CD21lo subset. Using this approach, we detected insulin-binding B cell bias towards CD27- memory and CD27+ memory subsets in pre-symptomatic type 1 diabetes donors. We took advantage of routine detection of Jo-1-binding B cells in Jo-1+ anti-histidyl tRNA synthetase syndrome patients to show that Jo-1-binding B cells and total B cells expanded 20-30-fold using this culture system. Overall, these studies highlight technology that is amenable to small numbers of cryopreserved peripheral blood mononuclear cells that enables interrogation of phenotypic and repertoire attributes of ASBCs derived from autoimmune patients.

Immunoglobulin D and its encoding genes: An updated review.

Since the identification of a functional Cδ gene in ostriches, immunoglobulin (Ig) D has been considered to be an extremely evolutionarily conserved Ig isotype besides the IgM found in all classes of jawed vertebrates. However, in contrast to IgM (which remains stable over evolutionary time), IgD shows considerable structural plasticity among vertebrate species and, moreover, its functions are far from elucidated even in humans and mice. Recently, several studies have shown that high expression of the IgD-B-cell receptor (IgD-BCR) may help physiologically autoreactive B cells survive in peripheral lymphoid tissues thanks to unresponsiveness to self-antigens and help their entry into germinal centers to "redeem" autoreactivity via somatic hypermutation. Other studies have demonstrated that secreted IgD may enhance mucosal homeostasis and immunity by linking B cells with basophils to optimize T-helper-2 cell-mediated responses and to constrain IgE-mediated basophil degranulation. Herein, we review the new discoveries on IgD-encoding genes in jawed vertebrates in the past decade. We also highlight advances in the functions of the IgD-BCR and secreted IgD in humans and mice.

Shared B cell memory to coronaviruses and other pathogens varies in human age groups and tissues.

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

Single-cell analysis of human B cell maturation predicts how antibody class switching shapes selection dynamics.

Protective humoral memory forms in secondary lymphoid organs where B cells undergo affinity maturation and differentiation into memory or plasma cells. Here, we provide a comprehensive roadmap of human B cell maturation with single-cell transcriptomics matched with bulk and single-cell antibody repertoires to define gene expression, antibody repertoires, and clonal sharing of B cell states at single-cell resolution, including memory B cell heterogeneity that reflects diverse functional and signaling states. We reconstruct gene expression dynamics during B cell activation to reveal a pre-germinal center state primed to undergo class switch recombination and dissect how antibody class-dependent gene expression in germinal center and memory B cells is linked with a distinct transcriptional wiring with potential to influence their fate and function. Our analyses reveal the dynamic cellular states that shape human B cell-mediated immunity and highlight how antibody isotype may play a role during their antibody-based selection.

Immunoglobulins in teleosts.

Immunoglobulins are glycoproteins which are produced as membrane-bound receptors on B-cells or in a secreted form, known as antibodies. In teleosts, three immunoglobulin isotypes, IgM, IgT, and IgD, are present, each comprising two identical heavy and two identical light polypeptide chains. The basic mechanisms for generation of immunoglobulin diversity are similar in teleosts and higher vertebrates. The B-cell pre-immune repertoire is diversified by VDJ recombination, junctional flexibility, addition of nucleotides, and combinatorial association of light and heavy chains, while the post-immune repertoire undergoes somatic hypermutation during clonal expansion. Typically, the teleost immunoglobulin heavy chain gene complex has a modified translocon arrangement where the Dτ-Jτ-Cτ cluster of IgT is generally located between the variable heavy chain (VH) region and the Dµ/δ-Jµ/δ-Cµ-Cδ gene segments, or within the set of VH gene segments. However, multiple genome duplication and deletion events and loss of some individual genes through evolution has complicated the IgH gene organization. The IgH gene arrangement allows the expression of either IgT or IgM/IgD. Alternative splicing is responsible for the regulation of IgM/IgD expression and the secreted versus transmembrane forms of IgT, IgD, and IgM. The overall structure of IgM and IgT is usually conserved across species, whereas IgD has a large variety of structures. IgM is the main effector molecule in both systemic and mucosal immunity and shows a broad range of concentrations in different teleost species. Although IgM is usually present in higher concentrations under normal conditions, IgT is considered the main mucosal Ig.

Four immunoglobulin isotypes and IgD splice variants in urodele amphibians.

Until recently, different families of urodele amphibians were thought to express distinct subsets of immunoglobulin (Ig) isotypes. In this study, we explored cDNAs encoding Ig heavy-chains (H-chains) in three species of urodele amphibians. We found that Cynops pyrrhogaster, Pleurodeles waltl, and Ambystoma mexicanum each carry genes encoding four Ig H-chain isotypes, including IgM, IgY, IgD, and IgX, similar to those found in anuran amphibians. We also found that urodele IgDs have a long constant region similar to those found in anuran, reptiles, and bony fishes. We also found several putative IgD splice variants. Our findings indicated that P. waltl IgP is not a novel isotype but an IgD splice variant. Altogether, our findings indicate that IgD splice variants may be universally expressed among amphibian species.

Preimmune Control of the Variance of TCR CDR-B3: Insights Gained From Germline Replacement of a TCR Dß Gene Segment With an Ig D H Gene Segment.

We have previously shown that the sequence of the immunoglobulin diversity gene segment (D H ) helps dictate the structure and composition of complementarity determining region 3 of the immunoglobulin heavy chain (CDR-H3). In order to test the role of germline D sequence on the diversity of the preimmune TCRß repertoire of T cells, we generated a mouse with a mutant TCRß DJC locus wherein the Dß2-Jß2 gene segment cluster was deleted and the remaining diversity gene segment, Dß1 (IMGT:TRDB1), was replaced with DSP2.3 (IMGT:IGHD2-02), a commonly used B cell immunoglobulin D H gene segment. Crystallographic studies have shown that the length and thus structure of TCR CDR-B3 places amino acids at the tip of CDR-B3 in a position to directly interact with peptide bound to an MHC molecule. The length distribution of complementarity determining region 3 of the T cell receptor beta chain (CDR-B3) has been proposed to be restricted largely by MHC-specific selection, disfavoring CDR-B3 that are too long or too short. Here we show that the mechanism of control of CDR-B3 length depends on the Dß sequence, which in turn dictates exonucleolytic nibbling. By contrast, the extent of N addition and the variance of created CDR3 lengths are regulated by the cell of origin, the thymocyte. We found that the sequence of the D and control of N addition collaborate to bias the distribution of CDR-B3 lengths in the pre-immune TCR repertoire and to focus the diversity provided by N addition and the sequence of the D on that portion of CDR-B3 that is most likely to interact with the peptide that is bound to the presenting MHC.

A novel transgenic mouse strain expressing PKCßII demonstrates expansion of B1 and marginal zone B cell populations.

Protein kinase Cß (PKCß) expressed in mammalian cells as two splice variants, PKCßI and PKCßII, functions in the B cell receptor (BCR) signaling pathway and contributes to B cell development. We investigated the relative role of PKCßII in B cells by generating transgenic mice where expression of the transgene is directed to these cells using the Eµ promoter (Eµ-PKCßIItg). Our findings demonstrate that homozygous Eµ-PKCßIItg mice displayed a shift from IgD+IgMdim toward IgDdimIgM+ B cell populations in spleen, peritoneum and peripheral blood. Closer examination of these tissues revealed respective expansion of marginal zone (MZ)-like B cells (IgD+IgM+CD43negCD21+CD24+), increased populations of B-1 cells (B220+IgDdimIgM+CD43+CD24+CD5+), and higher numbers of immature B cells (IgDdimIgMdimCD21neg) at the expense of mature B cells (IgD+IgM+CD21+). Therefore, the overexpression of PKCßII, which is a phenotypic feature of chronic lymphocytic leukaemia cells, can skew B cell development in mice, most likely as a result of a regulatory influence on BCR signaling.

Identification of an IgD/IgT chimera in the European sea bass (Dicentrarchus labrax L.).

Three classes of immunoglobulins have been identified in Teleosts: IgM, IgT/Z and IgD. They are fundamental for fish immune responses and, therefore, their functional activities are heavily investigated. In this paper, we describe the identification of a new IgD/IgT chimera in sea bass (Dicentrarchus labrax) from a gills transcriptome. This transcript joined the first six constant domains of the IgD chain with the two terminal constant domains of IgT, generating a long in-frame coding sequence with a junction between the canonical δ6 exon splicing donor site and the τ3 exon splicing acceptor site. Studies performed on genomic DNA confirmed the presence of the sequence and identifies and intronic region of 656 bp within this joining region. The basal expression of the IgD/IgT chimera was investigated both in silico and in vivo: high level of expression was found in gills, gut and head kidney. Moreover, IgD/IgT transcripts were up-regulated after in vitro stimulation of sea bass HK leukocytes with LPS. The IgD/IgT chimera was found also in two congener species, Morone saxatilis and Morone chrysops. It is not possible to have a precise idea on the evolutionary scenario that lead to the appearance of this sequence due to the lack of genomic information, but we could speculate that an ancestral duplication of the entire IgH locus was followed by the chimerization of Cδ/Cτ in one of the two loci. Finally, the IgD/IgT high basal expression in tissues and organs fundamental for sea bass immune response and its modulation after LPS stimulation provide a very preliminary indication that this unusual Ig variant could have a functional activity.

Differential Expression of IgM and IgD Discriminates Two Subpopulations of Human Circulating IgM+IgD+CD27+ B Cells That Differ Phenotypically, Functionally, and Genetically.

The origin and function of blood IgM+IgD+CD27+ B cells is controversial, and they are considered a heterogeneous population. Previous staining of circulating B cells of healthy donors with rotavirus fluorescent virus-like particles allowed us to differentiate two subsets of IgM+IgD+CD27+: IgMhi and IgMlo B cells. Here, we confirmed this finding and compared the phenotype, transcriptome, in vitro function, and Ig gene repertoire of these two subsets. Eleven markers phenotypically discriminated both subsets (CD1c, CD69, IL21R, CD27, MTG, CD45RB, CD5, CD184, CD23, BAFFR, and CD38) with the IgMhi phenotypically resembling previously reported marginal zone B cells and the IgMlo resembling both naïve and memory B cells. Transcriptomic analysis showed that both subpopulations clustered close to germinal center-experienced IgM only B cells with a Principal Component Analysis, but differed in expression of 78 genes. Moreover, IgMhi B cells expressed genes characteristic of previously reported marginal zone B cells. After stimulation with CpG and cytokines, significantly (p < 0.05) higher frequencies (62.5%) of IgMhi B cells proliferated, compared with IgMlo B cells (35.37%), and differentiated to antibody secreting cells (14.22% for IgMhi and 7.19% for IgMlo). IgMhi B cells had significantly (p < 0.0007) higher frequencies of mutations in IGHV and IGKV regions, IgMlo B cells had higher usage of IGHJ6 genes (p < 0.0001), and both subsets differed in their HCDR3 properties. IgMhi B cells shared most of their shared IGH clonotypes with IgM only memory B cells, and IgMlo B cells with IgMhi B cells. These results support the notion that differential expression of IgM and IgD discriminates two subpopulations of human circulating IgM+IgD+CD27+ B cells, with the IgMhi B cells having similarities with previously described marginal zone B cells that passed through germinal centers, and the IgMlo B cells being the least differentiated amongst the IgM+CD27+ subsets.

Automated analysis of immunosequencing datasets reveals novel immunoglobulin D genes across diverse species.

Immunoglobulin genes are formed through V(D)J recombination, which joins the variable (V), diversity (D), and joining (J) germline genes. Since variations in germline genes have been linked to various diseases, personalized immunogenomics focuses on finding alleles of germline genes across various patients. Although reconstruction of V and J genes is a well-studied problem, the more challenging task of reconstructing D genes remained open until the IgScout algorithm was developed in 2019. In this work, we address limitations of IgScout by developing a probabilistic MINING-D algorithm for D gene reconstruction, apply it to hundreds of immunosequencing datasets from multiple species, and validate the newly inferred D genes by analyzing diverse whole genome sequencing datasets and haplotyping heterozygous V genes.