Epstein-Barr Virus and immune status imprint the immunogenomics of non-Hodgkin lymphomas occurring in immune-suppressed environments.

Non-Hodgkin lymphomas (NHL) commonly occur in immune-deficient (ID) patients, both HIV-infected and transplanted, and are often EBV-driven with cerebral localization, raising the question of tumor immunogenicity, a critical issue for treatment responses. We investigated the immunogenomics of 68 lymphoproliferative disorders from 51 ID (34 posttransplant, 17 HIV+) and 17 immunocompetent patients. Overall, 72% were Large B Cells Lymphoma (LBCL) and 25% were primary central-nervous-system lymphoma (PCNSL) while 40% were EBV-positive. Tumor whole-exome and RNA sequencing, along with a bioinformatics pipeline allowed analysis of tumor mutational burden (TMB), tumor landscape and microenvironment (TME) and prediction of tumor neoepitopes. Both TMB (2.2 vs 3.4/Mb, p=0.001) and neoepitopes numbers (40 vs 200, p=0.00019) were lower in EBVpositive than in EBV-negative NHL, regardless of the immune status. In contrast both EBV and the immune status influenced the tumor mutational profile, with HNRNPF and STAT3 mutations exclusively observed in EBV-positive and ID NHL, respectively. Peripheral blood T-cell responses against tumor neoepitopes were detected in all EBV-negative cases but in only half EBV-positive ones, including responses against IgH-derived MHC-class-II restricted neoepitopes. The TME analysis showed higher CD8 T cell infiltrates in EBVpositive vs EBV-negative NHL, together with a more tolerogenic profile composed of Tregs, type-M2 macrophages and an increased expression of negative immune-regulators. Our results highlight that the immunogenomics of NHL in patients with immunodeficiency primarily relies on the tumor EBV status, while T cell recognition of tumor- and IgH-specific neoepitopes is conserved in EBV-negative patients, offering potential opportunities for future T cell-based immune therapies.

Reconstructing B cell lineage trees with minimum spanning tree and genotype abundances.

B cell receptor (BCR) genes exposed to an antigen undergo somatic hypermutations and Darwinian antigen selection, generating a large BCR-antibody diversity. This process, known as B cell affinity maturation, increases antibody affinity, forming a specific B cell lineage that includes the unmutated ancestor and mutated variants. In a B cell lineage, cells with a higher antigen affinity will undergo clonal expansion, while those with a lower affinity will not proliferate and probably be eliminated. Therefore, cellular (genotype) abundance provides a valuable perspective on the ongoing evolutionary process. Phylogenetic tree inference is often used to reconstruct B cell lineage trees and represents the evolutionary dynamic of BCR affinity maturation. However, such methods should process B-cell population data derived from experimental sampling that might contain different cellular abundances. There are a few phylogenetic methods for tracing the evolutionary events occurring in B cell lineages; best-performing solutions are time-demanding and restricted to analysing a reduced number of sequences, while time-efficient methods do not consider cellular abundances. We propose ClonalTree, a low-complexity and accurate approach to construct B-cell lineage trees that incorporates genotype abundances into minimum spanning tree (MST) algorithms. Using both simulated and experimental data, we demonstrate that ClonalTree outperforms MST-based algorithms and achieves a comparable performance to a method that explores tree-generating space exhaustively. Furthermore, ClonalTree has a lower running time, being more convenient for building B-cell lineage trees from high-throughput BCR sequencing data, mainly in biomedical applications, where a lower computational time is appreciable. It is hundreds to thousands of times faster than exhaustive approaches, enabling the analysis of a large set of sequences within minutes or seconds and without loss of accuracy. The source code is freely available at github.com/julibinho/ClonalTree.

Higher-order connections between stereotyped subsets: implications for improved patient classification in CLL.

Chronic lymphocytic leukemia (CLL) is characterized by the existence of subsets of patients with (quasi)identical, stereotyped B-cell receptor (BcR) immunoglobulins. Patients in certain major stereotyped subsets often display remarkably consistent clinicobiological profiles, suggesting that the study of BcR immunoglobulin stereotypy in CLL has important implications for understanding disease pathophysiology and refining clinical decision-making. Nevertheless, several issues remain open, especially pertaining to the actual frequency of BcR immunoglobulin stereotypy and major subsets, as well as the existence of higher-order connections between individual subsets. To address these issues, we investigated clonotypic IGHV-IGHD-IGHJ gene rearrangements in a series of 29 856 patients with CLL, by far the largest series worldwide. We report that the stereotyped fraction of CLL peaks at 41% of the entire cohort and that all 19 previously identified major subsets retained their relative size and ranking, while 10 new ones emerged; overall, major stereotyped subsets had a cumulative frequency of 13.5%. Higher-level relationships were evident between subsets, particularly for major stereotyped subsets with unmutated IGHV genes (U-CLL), for which close relations with other subsets, termed "satellites," were identified. Satellite subsets accounted for 3% of the entire cohort. These results confirm our previous notion that major subsets can be robustly identified and are consistent in relative size, hence representing distinct disease variants amenable to compartmentalized research with the potential of overcoming the pronounced heterogeneity of CLL. Furthermore, the existence of satellite subsets reveals a novel aspect of repertoire restriction with implications for refined molecular classification of CLL.

A multi-objective based clustering for inferring BCR clonal lineages from high-throughput B cell repertoire data.

The adaptive B cell response is driven by the expansion, somatic hypermutation, and selection of B cell clonal lineages. A high number of clonal lineages in a B cell population indicates a highly diverse repertoire, while clonal size distribution and sequence diversity reflect antigen selective pressure. Identifying clonal lineages is fundamental to many repertoire studies, including repertoire comparisons, clonal tracking, and statistical analysis. Several methods have been developed to group sequences from high-throughput B cell repertoire data. Current methods use clustering algorithms to group clonally-related sequences based on their similarities or distances. Such approaches create groups by optimizing a single objective that typically minimizes intra-clonal distances. However, optimizing several objective functions can be advantageous and boost the algorithm convergence rate. Here we propose MobiLLe, a new method based on multi-objective clustering. Our approach requires V(D)J annotations to obtain the initial groups and iteratively applies two objective functions that optimize cohesion and separation within clonal lineages simultaneously. We show that our method greatly improves clonal lineage grouping on simulated benchmarks with varied mutation rates compared to other tools. When applied to experimental repertoires generated from high-throughput sequencing, its clustering results are comparable to the most performing tools and can reproduce the results of previous publications. The method based on multi-objective clustering can accurately identify clonally-related antibody sequences and presents the lowest running time among state-of-art tools. All these features constitute an attractive option for repertoire analysis, particularly in the clinical context. MobiLLe can potentially help unravel the mechanisms involved in developing and evolving B cell malignancies.

A mitochondrial origin for frontotemporal dementia and amyotrophic lateral sclerosis through CHCHD10 involvement.

Mitochondrial DNA instability disorders are responsible for a large clinical spectrum, among which amyotrophic lateral sclerosis-like symptoms and frontotemporal dementia are extremely rare. We report a large family with a late-onset phenotype including motor neuron disease, cognitive decline resembling frontotemporal dementia, cerebellar ataxia and myopathy. In all patients, muscle biopsy showed ragged-red and cytochrome c oxidase-negative fibres with combined respiratory chain deficiency and abnormal assembly of complex V. The multiple mitochondrial DNA deletions found in skeletal muscle revealed a mitochondrial DNA instability disorder. Patient fibroblasts present with respiratory chain deficiency, mitochondrial ultrastructural alterations and fragmentation of the mitochondrial network. Interestingly, expression of matrix-targeted photoactivatable GFP showed that mitochondrial fusion was not inhibited in patient fibroblasts. Using whole-exome sequencing we identified a missense mutation (c.176C>T; p.Ser59Leu) in the CHCHD10 gene that encodes a coiled-coil helix coiled-coil helix protein, whose function is unknown. We show that CHCHD10 is a mitochondrial protein located in the intermembrane space and enriched at cristae junctions. Overexpression of a CHCHD10 mutant allele in HeLa cells led to fragmentation of the mitochondrial network and ultrastructural major abnormalities including loss, disorganization and dilatation of cristae. The observation of a frontotemporal dementia-amyotrophic lateral sclerosis phenotype in a mitochondrial disease led us to analyse CHCHD10 in a cohort of 21 families with pathologically proven frontotemporal dementia-amyotrophic lateral sclerosis. We identified the same missense p.Ser59Leu mutation in one of these families. This work opens a novel field to explore the pathogenesis of the frontotemporal dementia-amyotrophic lateral sclerosis clinical spectrum by showing that mitochondrial disease may be at the origin of some of these phenotypes.

Replication fork reversal after replication-transcription collision.

Replication fork arrest is a recognized source of genetic instability, and transcription is one of the most prominent causes of replication impediment. We analyze here the requirement for recombination proteins in Escherichia coli when replication-transcription head-on collisions are induced at a specific site by the inversion of a highly expressed ribosomal operon (rrn). RecBC is the only recombination protein required for cell viability under these conditions of increased replication-transcription collisions. In its absence, fork breakage occurs at the site of collision, and the resulting linear DNA is not repaired and is slowly degraded by the RecJ exonuclease. Lethal fork breakage is also observed in cells that lack RecA and RecD, i.e. when both homologous recombination and the potent exonuclease V activity of the RecBCD complex are inactivated, with a slow degradation of the resulting linear DNA by the combined action of the RecBC helicase and the RecJ exonuclease. The sizes of the major linear fragments indicate that DNA degradation is slowed down by the encounter with another rrn operon. The amount of linear DNA decreases nearly two-fold when the Holliday junction resolvase RuvABC is inactivated in recB, as well as in recA recD mutants, indicating that part of the linear DNA is formed by resolution of a Holliday junction. Our results suggest that replication fork reversal occurs after replication-transcription head-on collision, and we propose that it promotes the action of the accessory replicative helicases that dislodge the obstacle.

Partial deletions of the GRN gene are a cause of frontotemporal lobar degeneration.

Mutations in the progranulin gene (GRN) are an important cause of frontotemporal lobar degeneration (FTLD). Most known GRN mutations are null mutations, such as nonsense and frameshift mutations, which create a premature stop codon resulting in loss of function of the progranulin protein. Complete or near-complete genomic GRN deletions have also been found in three families, but heterozygous partial deletions that remove only one or two exons have not been reported to date. In this study, we analysed three unrelated FTLD patients with low plasma progranulin levels but no point GRN mutations by multiplex ligation-dependent probe amplification (MLPA) and quantitative multiplex polymerase chain reaction of short fluorescent fragments (QMPSF). We detected two heterozygous partial GRN deletions in two patients. One deletion removed exon 1 and part of intron 1. The second deletion was complex: it removed 1,410 bp extending from the part of intron 1 to the part of exon 3, with a small 5-bp insertion at the breakpoint junction (c.-7-1121_159delinsGATCA). Our findings illustrate the usefulness of a quantitative analysis in addition to GRN gene sequencing for a comprehensive genetic diagnosis of FTLD, particularly in patients with low plasma progranulin levels.

The helicases DinG, Rep and UvrD cooperate to promote replication across transcription units in vivo.

How living cells deal with head-on collisions of the replication and transcription complexes has been debated for a long time. Even in the widely studied model bacteria Escherichia coli, the enzymes that take care of such collisions are still unknown. We report here that in vivo, the DinG, Rep and UvrD helicases are essential for efficient replication across highly transcribed regions. We show that when rRNA operons (rrn) are inverted to face replication, the viability of the dinG mutant is affected and over-expression of RNase H rescues the growth defect, showing that DinG acts in vivo to remove R-loops. In addition, DinG, Rep and UvrD exert a common function, which requires the presence of two of these three helicases. After replication blockage by an inverted rrn, Rep in conjunction with DinG or UvrD removes RNA polymerase, a task that is fulfilled in its absence by the SOS-induced DinG and UvrD helicases. Finally, Rep and UvrD also act at inverted sequences other than rrn, and promote replication through highly transcribed regions in wild-type E. coli.

The immunopathological landscape of human pre-TCRα deficiency: From rare to common variants.

We describe humans with rare biallelic loss-of-function PTCRA variants impairing pre-α T cell receptor (pre-TCRα) expression. Low circulating naive αß T cell counts at birth persisted over time, with normal memory αß and high γδ T cell counts. Their TCRα repertoire was biased, which suggests that noncanonical thymic differentiation pathways can rescue αß T cell development. Only a minority of these individuals were sick, with infection, lymphoproliferation, and/or autoimmunity. We also report that 1 in 4000 individuals from the Middle East and South Asia are homozygous for a common hypomorphic PTCRA variant. They had normal circulating naive αß T cell counts but high γδ T cell counts. Although residual pre-TCRα expression drove the differentiation of more αß T cells, autoimmune conditions were more frequent in these patients compared with the general population.

ViCloD, an interactive web tool for visualizing B cell repertoires and analyzing intraclonal diversities: application to human B-cell tumors.

High throughput sequencing of adaptive immune receptor repertoire (AIRR-seq) has provided numerous human immunoglobulin (IG) sequences allowing specific B cell receptor (BCR) studies such as the antigen-driven evolution of antibodies (soluble forms of the membrane-bound IG part of the BCR). AIRR-seq data allows researchers to examine intraclonal differences caused primarily by somatic hypermutations in IG genes and affinity maturation. Exploring this essential adaptive immunity process could help elucidate the generation of antibodies with high affinity or broadly neutralizing activities. Retracing their evolutionary history could also clarify how vaccines or pathogen exposition drive the humoral immune response, and unravel the clonal architecture of B cell tumors. Computational methods are necessary for large-scale analysis of AIRR-seq properties. However, there is no efficient and interactive tool for analyzing intraclonal diversity, permitting users to explore adaptive immune receptor repertoires in biological and clinical applications. Here we present ViCloD, a web server for large-scale visual analysis of repertoire clonality and intraclonal diversity. ViCloD uses preprocessed data in the format defined by the Adaptive Immune Receptor Repertoire (AIRR) Community. Then, it performs clonal grouping and evolutionary analyses, producing a collection of useful plots for clonal lineage inspection. The web server presents diverse functionalities, including repertoire navigation, clonal abundance analysis, and intraclonal evolutionary tree reconstruction. Users can download the analyzed data in different table formats and save the generated plots as images. ViCloD is a simple, versatile, and user-friendly tool that can help researchers and clinicians to analyze B cell intraclonal diversity. Moreover, its pipeline is optimized to process hundreds of thousands of sequences within a few minutes, allowing an efficient investigation of large and complex repertoires.

Immunoglobulin Gene Mutational Status Assessment by Next Generation Sequencing in Chronic Lymphocytic Leukemia.

B cell receptor (BcR) immunoglobulins (IG) display a tremendous diversity due to complex DNA rearrangements, the V(D)J recombination, further enhanced by the somatic hypermutation process. In chronic lymphocytic leukemia (CLL), the mutational load of the clonal BcR IG expressed by the leukemic cells constitutes an important prognostic and predictive biomarker. Here, we provide a reliable methodology capable of determining the mutational status of IG genes in CLL using high-throughput sequencing, starting from leukemic cell DNA or RNA.

Evidence of somatic hypermutation in the antigen binding sites of patients with CLL harboring IGHV genes with 100% germline identity.

Classification of patients with chronic lymphocytic leukemia (CLL) based on the somatic hypermutation (SHM) status of the clonotypic immunoglobulin heavy variable (IGHV) gene has established predictive and prognostic relevance. The SHM status is assessed based on the number of mutations within the IG heavy variable domain sequence, albeit only over the rearranged IGHV gene excluding the variable heavy complementarity determining region 3 (VH CDR3). This may lead to an underestimation of the actual impact of SHM, in fact overlooking the most critical region for antigen-antibody interactions, i.e. the VH CDR3. Here we investigated whether SHM may be present within the VH CDR3 of cases bearing 'truly unmutated' IGHV genes (i.e. 100% germline identity across VH FR1-VH FR3) employing Next Generation Sequencing. We studied 16 patients bearing a 'truly unmutated' CLL clone assigned to stereotyped subsets #1 (n=12) and #6 (n=4). We report the existence of SHM within the germline-encoded 3'IGHV, IGHD, 5'IGHJ regions of the VH CDR3 in both the main IGHV-IGHD-IGHJ gene clonotype and its variants. Recurrent somatic mutations were identified between different patients of the same subset, supporting the notion that they represent true mutational events rather than technical artefacts; moreover, they were located adjacent to/within AID hotspots, pointing to SHM as the underlying mechanism. In conclusion, we provide immunogenetic evidence for intra-VH CDR3 variations, attributed to SHM, in CLL patients carrying 'truly unmutated' IGHV genes. Although the clinical implications of this observation remain to be defined, our findings offer a new perspective into the immunobiology of CLL, alluding to the operation of VH CDR3-restricted SHM in U-CLL.

Immunoglobulin gene analysis in chronic lymphocytic leukemia in the era of next generation sequencing.

Twenty years after landmark publications, there is a consensus that the somatic hypermutation (SHM) status of the clonotypic immunoglobulin heavy variable (IGHV) gene is an important cornerstone for accurate risk stratification and therapeutic decision-making in patients with chronic lymphocytic leukemia (CLL). The IGHV SHM status has traditionally been determined by conventional Sanger sequencing. However, NGS has heralded a new era in medical diagnostics and immunogenetic analysis is following this trend. There is indeed a growing demand for shifting practice and using NGS for IGHV gene SHM assessment, although it is debatable whether it is always justifiable, at least taking into account financial considerations for laboratories with limited resources. Nevertheless, as this analysis impacts on treatment decisions, standardization of both technical aspects, and data interpretation becomes essential. Also, the need for establishing new recommendations and providing dedicated education and training on NGS-based immunogenetics is greater than ever before. Here we address potential and challenges of NGS-based immunogenetics in CLL. We are convinced that this perspective helps the hematological community to better understand the pros and cons of this new technological development for CLL patient management.