InterTADs: integration of multi-omics data on topologically associated domains, application to chronic lymphocytic leukemia

ABSTRACT The integration of multi-omics data can greatly facilitate the advancement of research in Life Sciences by highlighting new interactions. However, there is currently no widespread procedure for meaningful multi-omics data integration. Here, we present a robust framework, called InterTADs, for integrating multi-omics data derived from the same sample, and considering the chromatin configuration of the genome, i.e. the topologically associating domains (TADs). Following the integration process, statistical analysis highlights the differences between the groups of interest (normal versus cancer cells) relating to (i) independent and (ii) integrated events through TADs. Finally, enrichment analysis using KEGG database, Gene Ontology and transcription factor binding sites and visualization approaches are available. We applied InterTADs to multi-omics datasets from 135 patients with chronic lymphocytic leukemia (CLL) and found that the integration through TADs resulted in a dramatic reduction of heterogeneity compared to individual events. Significant differences for individual events and on TADs level were identified between patients differing in the somatic hypermutation status of the clonotypic immunoglobulin genes, the core biological stratifier in CLL, attesting to the biomedical relevance of InterTADs. In conclusion, our approach suggests a new perspective towards analyzing multi-omics data, by offering reasonable execution time, biological benchmarking and potentially contributing to pattern discovery through TADs.

Association of mutation signature effectuating processes with mutation hotspots in driver genes and non-coding regions

Cancer driving mutations are difficult to identify especially in the non-coding part of the genome. Here, we present sigDriver, an algorithm dedicated to call driver mutations. Using 3813 whole-genome sequenced tumors from International Cancer Genome Consortium, The Cancer Genome Atlas Program, and a childhood pan-cancer cohort, we employ mutational signatures based on single-base substitution in the context of tri- and penta-nucleotide motifs for hotspot discovery. Knowledge-based annotations on mutational hotspots reveal enrichment in coding regions and regulatory elements for 6 mutational signatures, including APOBEC and somatic hypermutation signatures. APOBEC activity is associated with 32 hotspots of which 11 are known and 11 are putative regulatory drivers. Somatic single nucleotide variants clusters detected at hypermutation-associated hotspots are distinct from translocation or gene amplifications. Patients carrying APOBEC induced PIK3CA driver mutations show lower occurrence of signature SBS39. In summary, sigDriver uncovers mutational processes associated with known and putative tumor drivers and hotspots particularly in the non-coding regions of the genome.

Tandem Substitutions in Somatic Hypermutation

Upon antigen recognition, activation-induced cytosine deaminase initiates affinity maturation of the B-cell receptor by somatic hypermutation (SHM) through error-prone DNA repair pathways. SHM typically creates single nucleotide substitutions, but tandem substitutions may also occur. We investigated incidence and sequence context of tandem substitutions by massive parallel sequencing of V(D)J repertoires in healthy human donors. Mutation patterns were congruent with SHM-derived single nucleotide mutations, delineating initiation of the tandem substitution by AID. Tandem substitutions comprised 5,7% of AID-induced mutations. The majority of tandem substitutions represents single nucleotide juxtalocations of directly adjacent sequences. These observations were confirmed in an independent cohort of healthy donors. We propose a model where tandem substitutions are predominantly generated by translesion synthesis across an apyramidinic site that is typically created by UNG. During replication, apyrimidinic sites transiently adapt an extruded configuration, causing skipping of the extruded base. Consequent strand decontraction leads to the juxtalocation, after which exonucleases repair the apyramidinic site and any directly adjacent mismatched base pairs. The mismatch repair pathway appears to account for the remainder of tandem substitutions. Tandem substitutions may enhance affinity maturation and expedite the adaptive immune response by overcoming amino acid codon degeneracies or mutating two adjacent amino acid residues simultaneously.

Altered Germinal-Center Metabolism in B Cells in Autoimmunity

B lymphocytes play an important role in the pathophysiology of many autoimmune disorders by producing autoantibodies, secreting cytokines, and presenting antigens. B cells undergo extreme physiological changes as they develop and differentiate. Aberrant function in tolerogenic checkpoints and the metabolic state of B cells might be the contributing factors to the dysfunctionality of autoimmune B cells. Understanding B-cell metabolism in autoimmunity is important as it can give rise to new treatments. Recent investigations have revealed that alterations in metabolism occur in the activation of B cells. Several reports have suggested that germinal center (GC) B cells of individuals with systemic lupus erythematosus (SLE) have altered metabolic function. GCs are unique microenvironments in which the delicate and complex process of B-cell affinity maturation occurs through somatic hypermutation (SHM) and class switching recombination (CSR) and where Bcl6 tightly regulates B-cell differentiation into memory B-cells or plasma cells. GC B cells rely heavily on glucose, fatty acids, and oxidative phosphorylation (OXPHOS) for their energy requirements. However, the complicated association between GC B cells and their metabolism is still not clearly understood. Here, we review several studies of B-cell metabolism, highlighting the significant transformations that occur in GC progression, and suggest possible approaches that may be investigated to more precisely target aberrant B-cell metabolism in SLE.

USP10 regulates B cell response to SARS-CoV-2 or HIV-1 nanoparticle vaccines through deubiquitinating AID

Activation-induced cytidine deaminase (AID) initiates class-switch recombination and somatic hypermutation (SHM) in antibody genes. Protein expression and activity are tightly controlled by various mechanisms. However, it remains unknown whether a signal from the extracellular environment directly affects the AID activity in the nucleus where it works. Here, we demonstrated that a deubiquitinase USP10, which specifically stabilizes nuclear AID protein, can translocate into the nucleus after AKT-mediated phosphorylation at its T674 within the NLS domain. Interestingly, the signals from BCR and TLR1/2 synergistically promoted this phosphorylation. The deficiency of USP10 in B cells significantly decreased AID protein levels, subsequently reducing neutralizing antibody production after immunization with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or human immunodeficiency virus type 1 (HIV-1) nanoparticle vaccines. Collectively, we demonstrated that USP10 functions as an integrator for both BCR and TLR signals and directly regulates nuclear AID activity. Its manipulation could be used for the development of vaccines and adjuvants.

Structural Basis of Antibody Conformation and Stability Modulation by Framework Somatic Hypermutation

Accumulation of somatic hypermutation (SHM) is the primary mechanism to enhance the binding affinity of antibodies to antigens in vivo. However, the structural basis of the effects of many SHMs remains elusive. Here, we integrated atomistic molecular dynamics (MD) simulation and data mining to build a high-throughput structural bioinformatics pipeline to study the effects of individual and combination SHMs on antibody conformation, flexibility, stability, and affinity. By applying this pipeline, we characterized a common mechanism of modulation of heavy-light pairing orientation by frequent SHMs at framework positions 39H, 91H, 38L, and 87L through disruption of a conserved hydrogen-bond network. Q39LH alone and in combination with light chain framework 4 (FWR4L) insertions further modulated the elbow angle between variable and constant domains of many antibodies, resulting in improved binding affinity for a subset of anti-HIV-1 antibodies. Q39LH also alleviated aggregation induced by FWR4L insertion, suggesting remote epistasis between these SHMs. Altogether, this study provides tools and insights for understanding antibody affinity maturation and for engineering functionally improved antibodies.

Intrinsic Strand-Incision Activity of Human UNG: Implications for Nick Generation in Immunoglobulin Gene Diversification

Uracil arises in cellular DNA by cytosine (C) deamination and erroneous replicative incorporation of deoxyuridine monophosphate opposite adenine. The former generates C → thymine transition mutations if uracil is not removed by uracil-DNA glycosylase (UDG) and replaced by C by the base excision repair (BER) pathway. The primary human UDG is hUNG. During immunoglobulin gene diversification in activated B cells, targeted cytosine deamination by activation-induced cytidine deaminase followed by uracil excision by hUNG is important for class switch recombination (CSR) and somatic hypermutation by providing the substrate for DNA double-strand breaks and mutagenesis, respectively. However, considerable uncertainty remains regarding the mechanisms leading to DNA incision following uracil excision: based on the general BER scheme, apurinic/apyrimidinic (AP) endonuclease (APE1 and/or APE2) is believed to generate the strand break by incising the AP site generated by hUNG. We report here that hUNG may incise the DNA backbone subsequent to uracil excision resulting in a 3´-α,β-unsaturated aldehyde designated uracil-DNA incision product (UIP), and a 5´-phosphate. The formation of UIP accords with an elimination (E2) reaction where deprotonation of C2´ occurs via the formation of a C1´ enolate intermediate. UIP is removed from the 3´-end by hAPE1. This shows that the first two steps in uracil BER can be performed by hUNG, which might explain the significant residual CSR activity in cells deficient in APE1 and APE2.

Long-lasting germinal center responses to a priming immunization with continuous proliferation and somatic mutation

Germinal centers (GCs) are the engines of antibody evolution. Using HIV Env protein immunogen priming in rhesus monkeys (RM) followed by a long period without further immunization, we demonstrate GC B cells (BGC) lasted at least 6 months (29 weeks), all the while maintaining rapid proliferation. A 186-fold BGC cell increase was present by week 10 compared to a conventional immunization. Single cell transcriptional profiling revealed that both light zone and dark zone GC states were sustained throughout the 6 months. Antibody somatic hypermutation (SHM) of BGC cells continued to accumulate throughout the 29 week priming period, with evidence of selective pressure. Additionally, Env-binding BGC cells were still 49-fold above baseline 29 weeks after immunization, suggesting that they could be active for significantly longer periods of time. High titers of HIV neutralizing antibodies were generated after a single booster immunization. Fully glycosylated HIV trimer protein is a complex antigen, posing significant immunodominance challenges for B cells, among other difficulties. Memory B cells (BMem) generated under these long priming conditions had higher levels of SHM, and both BMem cells and antibodies were more likely to recognize non-immunodominant epitopes. Numerous BGC cell lineage phylogenies spanning the >6-month GC period were identified, demonstrating continuous GC activity and selection for at least 191 days, with no additional antigen exposure. A long prime, adjuvanted, slow delivery (12-day) immunization approach holds promise for difficult vaccine targets, and suggests that patience can have great value for tuning GCs to maximize antibody responses.

B Cells in Primary Membranous Nephropathy: Escape from Immune Tolerance and Implications for Patient Management

Membranous nephropathy (MN) is an important cause of nephrotic syndrome and chronic kidney disease (CKD) in adults. The pathogenic significance of B cells in MN is increasingly recognized, especially following the discovery of various autoantibodies that target specific podocytic antigens and the promising treatment responses seen with B cell depleting therapies. The presence of autoreactive B cells and autoantibodies that bind to antigens on podocyte surfaces are characteristic features of MN, and are the result of breaches in central and peripheral tolerance of B lymphocytes. These perturbations in B cell tolerance include altered B lymphocyte subsets, dysregulation of genes that govern immunoglobulin production, aberrant somatic hypermutation and co-stimulatory signalling, abnormal expression of B cell-related cytokines, and increased B cell infiltrates and organized tertiary lymphoid structures within the kidneys. An understanding of the role of B cell tolerance and homeostasis may have important implications for patient management in MN, as conventional immunosuppressive treatments and novel B cell-targeted therapies show distinct effects on proliferation, differentiation and reconstitution in different B cell subsets. Circulating B lymphocytes and related cytokines may serve as potential biomarkers for treatment selection, monitoring of therapeutic response and prediction of disease relapse. These recent advances in the understanding of B cell tolerance in MN have provided greater insight into its immunopathogenesis and potential novel strategies for disease monitoring and treatment.