Chromatin activation profiling of stereotyped chronic lymphocytic leukemias reveals a subset 8-specific signature.

The chromatin activation landscape of chronic lymphocytic leukemia (CLL) with stereotyped B-cell receptor immunoglobulin is currently unknown. In this study, we report the results of a whole-genome chromatin profiling of histone 3 lysine 27 acetylation of 22 CLLs from major subsets, which were compared against nonstereotyped CLLs and normal B-cell subpopulations. Although subsets 1, 2, and 4 did not differ much from their nonstereotyped CLL counterparts, subset 8 displayed a remarkably distinct chromatin activation profile. In particular, we identified 209 de novo active regulatory elements in this subset, which showed similar patterns with U-CLLs undergoing Richter transformation. These regions were enriched for binding sites of 9 overexpressed transcription factors. In 78 of 209 regions, we identified 113 candidate overexpressed target genes, 11 regions being associated with more than 2 adjacent genes. These included blocks of up to 7 genes, suggesting local coupregulation within the same genome compartment. Our findings further underscore the uniqueness of subset 8 CLL, notable for the highest risk of Richter's transformation among all CLLs and provide additional clues to decipher the molecular basis of its clinical behavior.

N-Glycosylation of the Ig Receptors Shapes the Antigen Reactivity in Chronic Lymphocytic Leukemia Subset #201.

Subset #201 is a clinically indolent subgroup of patients with chronic lymphocytic leukemia defined by the expression of stereotyped, mutated IGHV4-34/IGLV1-44 BCR Ig. Subset #201 is characterized by recurrent somatic hypermutations (SHMs) that frequently lead to the creation and/or disruption of N-glycosylation sites within the Ig H and L chain variable domains. To understand the relevance of this observation, using next-generation sequencing, we studied how SHM shapes the subclonal architecture of the BCR Ig repertoire in subset #201, particularly focusing on changes in N-glycosylation sites. Moreover, we profiled the Ag reactivity of the clonotypic BCR Ig expressed as rmAbs. We found that almost all analyzed cases from subset #201 carry SHMs potentially affecting N-glycosylation at the clonal and/or subclonal level and obtained evidence for N-glycan occupancy in SHM-induced novel N-glycosylation sites. These particular SHMs impact (auto)antigen recognition, as indicated by differences in Ag reactivity between the authentic rmAbs and germline revertants of SHMs introducing novel N-glycosylation sites in experiments entailing 1) flow cytometry for binding to viable cells, 2) immunohistochemistry against various human tissues, 3) ELISA against microbial Ags, and 4) protein microarrays testing reactivity against multiple autoantigens. On these grounds, N-glycosylation appears as relevant for the natural history of at least a fraction of Ig-mutated chronic lymphocytic leukemia. Moreover, subset #201 emerges as a paradigmatic case for the role of affinity maturation in the evolution of Ag reactivity of the clonotypic BCR Ig.

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.

High-Throughput immunogenetics for precision medicine in cancer.

Cancer is characterized by an extremely complex biological background, which hinders personalized therapeutic interventions. Precision medicine promises to overcome this obstacle through integrating information from different 'subsystems', including the host, the external environment, the tumor itself and the tumor micro-environment. Immunogenetics is an essential tool that allows dissecting both lymphoid cancer ontogeny at both a cell-intrinsic and a cell-extrinsic level, i.e. through characterizing micro-environmental interactions, with a view to precision medicine. This is particularly thanks to the introduction of powerful, high-throughput approaches i.e. next generation sequencing, which allow the comprehensive characterization of immune repertoires. Indeed, NGS immunogenetic analysis (Immune-seq) has emerged as key to both understanding cancer pathogenesis and improving the accuracy of clinical decision making in oncology. Immune-seq has applications in lymphoid malignancies, assisting in the diagnosis e.g. through differentiating from reactive conditions, as well as in disease monitoring through accurate assessment of minimal residual disease. Moreover, Immune-seq facilitates the study of T cell receptor clonal dynamics in critical clinical contexts, including transplantation as well as innovative immunotherapy for solid cancers. The clinical utility of Immune-seq represents the focus of the present contribution, where we highlight what can be achieved but also what must be addressed in order to maximally realize the promise of Immune-seq in precision medicine in cancer.

Higher-order immunoglobulin repertoire restrictions in CLL: the illustrative case of stereotyped subsets 2 and 169.

Chronic lymphocytic leukemia (CLL) major stereotyped subset 2 (IGHV3-21/IGLV3-21, ∼2.5% of all cases of CLL) is an aggressive disease variant, irrespective of the somatic hypermutation (SHM) status of the clonotypic IGHV gene. Minor stereotyped subset 169 (IGHV3-48/IGLV3-21, ∼0.2% of all cases of CLL) is related to subset 2, as it displays a highly similar variable antigen-binding site. We further explored this relationship through next-generation sequencing and crystallographic analysis of the clonotypic B-cell receptor immunoglobulin. Branching evolution of the predominant clonotype through intraclonal diversification in the context of ongoing SHM was evident in both heavy and light chain genes of both subsets. Molecular similarities between the 2 subsets were highlighted by the finding of shared SHMs within both the heavy and light chain genes in all analyzed cases at either the clonal or subclonal level. Particularly noteworthy in this respect was a ubiquitous SHM at the linker region between the variable and the constant domain of the IGLV3-21 light chains, previously reported as critical for immunoglobulin homotypic interactions underlying cell-autonomous signaling capacity. Notably, crystallographic analysis revealed that the IGLV3-21-bearing CLL subset 169 immunoglobulin retains the same geometry and contact residues for the homotypic intermolecular interaction observed in subset 2, including the SHM at the linker region, and, from a molecular standpoint, belong to a common structural mode of autologous recognition. Collectively, our findings document that stereotyped subsets 2 and 169 are very closely related, displaying shared immunoglobulin features that can be explained only in the context of shared functional selection.

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.

TRIP - T cell receptor/immunoglobulin profiler.

BACKGROUND: Antigen receptors are characterized by an extreme diversity of specificities, which poses major computational and analytical challenges, particularly in the era of high-throughput immunoprofiling by next generation sequencing (NGS). The T cell Receptor/Immunoglobulin Profiler (TRIP) tool offers the opportunity for an in-depth analysis based on the processing of the output files of the IMGT/HighV-Quest tool, a standard in NGS immunoprofiling, through a number of interoperable modules. These provide detailed information about antigen receptor gene rearrangements, including variable (V), diversity (D) and joining (J) gene usage, CDR3 amino acid and nucleotide composition and clonality of both T cell receptors (TR) and B cell receptor immunoglobulins (BcR IG), and characteristics of the somatic hypermutation within the BcR IG genes. TRIP is a web application implemented in R shiny. RESULTS: Two sets of experiments have been performed in order to evaluate the efficiency and performance of the TRIP tool. The first used a number of synthetic datasets, ranging from 250k to 1M sequences, and established the linear response time of the tool (about 6 h for 1M sequences processed through the entire BcR IG data pipeline). The reproducibility of the tool was tested comparing the results produced by the main TRIP workflow with the results from a previous pipeline used on the Galaxy platform. As expected, no significant differences were noted between the two tools; although the preselection process seems to be stricter within the TRIP pipeline, about 0.1% more rearrangements were filtered out, with no impact on the final results. CONCLUSIONS: TRIP is a software framework that provides analytical services on antigen receptor gene sequence data. It is accurate and contains functions for data wrangling, cleaning, analysis and visualization, enabling the user to build a pipeline tailored to their needs. TRIP is publicly available at https://bio.tools/TRIP_-_T-cell_Receptor_Immunoglobulin_Profiler .

Recent revelations and future directions using single-cell technologies in chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a clinically and biologically heterogeneous disease with varying outcomes. In the last decade, the application of next-generation sequencing technologies has allowed extensive mapping of disease-specific genomic, epigenomic, immunogenetic, and transcriptomic signatures linked to CLL pathogenesis. These technologies have improved our understanding of the impact of tumor heterogeneity and evolution on disease outcome, although they have mostly been performed on bulk preparations of nucleic acids. As a further development, new technologies have emerged in recent years that allow high-resolution mapping at the single-cell level. These include single-cell RNA sequencing for assessment of the transcriptome, both of leukemic and non-malignant cells in the tumor microenvironment; immunogenetic profiling of B and T cell receptor rearrangements; single-cell sequencing methods for investigation of methylation and chromatin accessibility across the genome; and targeted single-cell DNA sequencing for analysis of copy-number alterations and single nucleotide variants. In addition, concomitant profiling of cellular subpopulations, based on protein expression, can also be obtained by various antibody-based approaches. In this review, we discuss different single-cell sequencing technologies and how they have been applied so far to study CLL onset and progression, also in response to treatment. This latter aspect is particularly relevant considering that we are moving away from chemoimmunotherapy to targeted therapies, with a potentially distinct impact on clonal dynamics. We also discuss new possibilities, such as integrative multi-omics analysis, as well as inherent limitations of the different single-cell technologies, from sample preparation to data interpretation using available bioinformatic pipelines. Finally, we discuss future directions in this rapidly evolving field.

Context-dependent T-cell Receptor Gene Repertoire Profiles in Proliferations of T Large Granular Lymphocytes.

T cell large granular lymphocyte (T-LGL) lymphoproliferations constitute a disease spectrum ranging from poly/oligo to monoclonal. Boundaries within this spectrum of proliferations are not well established. T-LGL lymphoproliferations co-occur with a wide variety of other diseases ranging from autoimmune disorders, solid tumors, hematological malignancies, post solid organ, and hematopoietic stem cell transplantation, and can therefore arise as a consequence of a wide variety of antigenic triggers. Persistence of a dominant malignant T-LGL clone is established through continuous STAT3 activation. Using next-generation sequencing, we profiled a cohort of 27 well-established patients with T-LGL lymphoproliferations, aiming to identify the subclonal architecture of the T-cell receptor beta (TRB) chain gene repertoire. Moreover, we searched for associations between TRB gene repertoire patterns and clinical manifestations, with the ultimate objective of discriminating between T-LGL lymphoproliferations developing in different clinical contexts and/or displaying distinct clinical presentation. Altogether, our data demonstrates that the TRB gene repertoire of patients with T-LGL lymphoproliferations is context-dependent, displaying distinct clonal architectures in different settings. Our results also highlight that there are monoclonal T-LGL cells with or without STAT3 mutations that cause symptoms such as neutropenia on one end of a spectrum and reactive oligoclonal T-LGL lymphoproliferations on the other. Longitudinal analysis revealed temporal clonal dynamics and showed that T-LGL cells might arise as an epiphenomenon when co-occurring with other malignancies, possibly reactive toward tumor antigens.

Variant-Related Differences in Laboratory Biomarkers among Patients Affected with Alpha, Delta and Omicron: A Retrospective Whole Viral Genome Sequencing and Hospital-Setting Cohort Study.

During the COVID-19 pandemic, different SARS-CoV-2 variants of concern (VOC) with specific characteristics have emerged and spread worldwide. At the same time, clinicians routinely evaluate the results of certain blood tests upon patient admission as well as during hospitalization to assess disease severity and the overall patient status. In the present study, we searched for significant cell blood count and biomarker differences among patients affected with the Alpha, Delta and Omicron VOCs at admission. Data from 330 patients were retrieved regarding age, gender, VOC, cell blood count results (WBC, Neut%, Lymph%, Ig%, PLT), common biomarkers (D-dimers, urea, creatinine, SGOT, SGPT, CRP, IL-6, suPAR), ICU admission and death. Statistical analyses were performed using ANOVA, the Kruskal-Wallis test, two-way ANOVA, Chi-square, T-test, the Mann-Whitney test and logistic regression was performed where appropriate using SPSS v.28 and STATA 14. Age and VOC were significantly associated with hospitalization, whereas significant differences among VOC groups were found for WBC, PLT, Neut%, IL-6, creatinine, CRP, D-dimers and suPAR. Our analyses showed that throughout the current pandemic, not only the SARS-CoV-2 VOCs but also the laboratory parameters that are used to evaluate the patient's status at admission are subject to changes.

T cell receptor gene repertoire profiles in subgroups of patients with chronic lymphocytic leukemia bearing distinct genomic aberrations.

Background: Microenvironmental interactions of the malignant clone with T cells are critical throughout the natural history of chronic lymphocytic leukemia (CLL). Indeed, clonal expansions of T cells and shared clonotypes exist between different CLL patients, strongly implying clonal selection by antigens. Moreover, immunogenic neoepitopes have been isolated from the clonotypic B cell receptor immunoglobulin sequences, offering a rationale for immunotherapeutic approaches. Here, we interrogated the T cell receptor (TR) gene repertoire of CLL patients with different genomic aberration profiles aiming to identify unique signatures that would point towards an additional source of immunogenic neoepitopes for T cells. Experimental design: TR gene repertoire profiling using next generation sequencing in groups of patients with CLL carrying one of the following copy-number aberrations (CNAs): del(11q), del(17p), del(13q), trisomy 12, or gene mutations in TP53 or NOTCH1. Results: Oligoclonal expansions were found in all patients with distinct recurrent genomic aberrations; these were more pronounced in cases bearing CNAs, particularly trisomy 12, rather than gene mutations. Shared clonotypes were found both within and across groups, which appeared to be CLL-biased based on extensive comparisons against TR databases from various entities. Moreover, in silico analysis identified TR clonotypes with high binding affinity to neoepitopes predicted to arise from TP53 and NOTCH1 mutations. Conclusions: Distinct TR repertoire profiles were identified in groups of patients with CLL bearing different genomic aberrations, alluding to distinct selection processes. Abnormal protein expression and gene dosage effects associated with recurrent genomic aberrations likely represent a relevant source of CLL-specific selecting antigens.

Purpose-Built Immunoinformatics for BcR IG/TR Repertoire Data Analysis.

The study of antigen receptor gene repertoires using next-generation sequencing (NGS) technologies has disclosed an unprecedented depth of complexity, requiring novel computational and analytical solutions. Several bioinformatics workflows have been developed to this end, including the T-cell receptor/immunoglobulin profiler (TRIP), a web application implemented in R shiny, specifically designed for the purposes of comprehensive repertoire analysis, which is the focus of this chapter. TRIP has the potential to perform robust immunoprofiling analysis through the extraction and processing of the IMGT/HighV-Quest output, via a series of functions, ensuring the analysis of high-quality, biologically relevant data through a multilevel process of data filtering. Subsequently, it provides in-depth analysis of antigen receptor gene rearrangements, including (a) clonality assessment; (b) extraction of variable (V), diversity (D), and joining (J) gene repertoires; (c) CDR3 characterization at both the nucleotide and amino acid level; and (d) somatic hypermutation analysis, in the case of immunoglobulin gene rearrangements. Relevant to mention, TRIP enables a high level of customization through the integration of various options in key aspects of the analysis, such as clonotype definition and computation, hence allowing for flexibility without compromising on accuracy.

Lumbar spine intrathecal transplantation of neural precursor cells promotes oligodendrocyte proliferation in hot spots of chronic demyelination.

Experimental autoimmune encephalomyelitis (EAE) is a basic and reliable model used to study clinical and pathological hallmarks of multiple sclerosis (MS) in rodents. Several studies suggest neural precursor cells (NPCs) as a significant research tool while reporting that transplanted NPCs are a promising therapeutic approach to treating neurological disorders, such as MS. The main objective was to approach a preclinical, in vivo scenario of oligodendrogenesis with NPCs, targeting the main chronic demyelinated lumbosacral milieu of EAE, via the least invasive delivery method which is lumbar puncture. We utilized MOG35-55 peptide to induce EAE in C57BL/6 mice and prior to the acute relapse, we intervened with either the traceable GFP+ cellular therapy or saline solution in the intrathecal space of their lumbar spine. A BrdU injection, which enabled us to monitor endogenous proliferation, marked the endpoint 50 days post-induction (50 dpi). Neuropathology with high-throughput, triple immunofluorescent, and transmission electron microscopy (TEM) data were extracted and analyzed. The experimental treatment attenuated the chronic phase of EAE (50 dpi; score <1) following an acute, clinical relapse. Myelination and axonal integrity were rescued in the NPC-treated animals along with suppressed immune populations. The differentiation profile of the exogenous NPCs and endogenous BrdU+ cells was location-dependent where GFP+ -rich areas drove undifferentiated phenotypes toward the oligodendrocyte lineage. In situ oligodendrocyte enrichment was demonstrated through increased (p < 0.001) gap junction channels of Cx32 and Cx47, reliable markers for proliferative oligodendroglia syncytium. TEM morphometric analysis ultimately manifested an increased g-ratio in lumbosacral fibers of the recovered animals (p < 0.001). Herein, we suggest that a single, lumbar intrathecal administration of NPCs capacitated a viable cellular load and resulted in clinical and pathological amelioration, stimulating resident OPCs to overcome the remyelination failure in EAE demyelinating locale.

Immunoglobulin gene sequence analysis in chronic lymphocytic leukemia: the 2022 update of the recommendations by ERIC, the European Research Initiative on CLL.

The somatic hypermutation (SHM) status of the clonotypic immunoglobulin heavy variable (IGHV) gene is a critical biomarker for assessing the prognosis of patients with chronic lymphocytic leukemia (CLL). Importantly, independent studies have documented that IGHV SHM status is also a predictor of responses to therapy, including both chemoimmunotherapy (CIT) and novel, targeted agents. Moreover, immunogenetic analysis in CLL has revealed that different patients may express (quasi)identical, stereotyped B cell receptor immunoglobulin (BcR IG) and are classified into subsets based on this common feature. Patients in certain stereotyped subsets display consistent biology, clinical presentation, and outcome that are distinct from other patients, even with concordant IGHV gene SHM status. All of the above highlights the relevance of immunogenetic analysis in CLL, which is considered a cornerstone for accurate risk stratification and clinical decision making. Recommendations for robust immunogenetic analysis exist thanks to dedicated efforts by ERIC, the European Research Initiative on CLL, covering all test phases, from the pre-analytical and analytical to the post-analytical, pertaining to the analysis, interpretation, and reporting of the findings. That said, these recommendations apply to Sanger sequencing, which is increasingly being superseded by next generation sequencing (NGS), further underscoring the need for an update. Here, we present an overview of the clinical utility of immunogenetics in CLL and update our analytical recommendations with the aim to assist in the refined management of patients with CLL.

The TΑp63/BCL2 axis represents a novel mechanism of clinical aggressiveness in chronic lymphocytic leukemia.

The TA-isoform of the p63 transcription factor (TAp63) has been reported to contribute to clinical aggressiveness in chronic lymphocytic leukemia (CLL) in a hitherto elusive way. Here, we sought to further understand and define the role of TAp63 in the pathophysiology of CLL. First, we found that elevated TAp63 expression levels are linked with adverse clinical outcomes, including disease relapse and shorter time-to-first treatment and overall survival. Next, prompted by the fact that TAp63 participates in an NF-κB/TAp63/BCL2 antiapoptotic axis in activated mature, normal B cells, we explored molecular links between TAp63 and BCL2 also in CLL. We documented a strong correlation at both the protein and the messenger RNA (mRNA) levels, alluding to the potential prosurvival role of TAp63. This claim was supported by inducible downregulation of TAp63 expression in the MEC1 CLL cell line using clustered regularly interspaced short palindromic repeats (CRISPR) system, which resulted in downregulation of BCL2 expression. Next, using chromatin immunoprecipitation (ChIP) sequencing, we examined whether BCL2 might constitute a transcriptional target of TAp63 and identified a significant binding profile of TAp63 in the BCL2 gene locus, across a genomic region previously characterized as a super enhancer in CLL. Moreover, we identified high-confidence TAp63 binding regions in genes mainly implicated in immune response and DNA-damage procedures. Finally, we found that upregulated TAp63 expression levels render CLL cells less responsive to apoptosis induction with the BCL2 inhibitor venetoclax. On these grounds, TAp63 appears to act as a positive modulator of BCL2, hence contributing to the antiapoptotic phenotype that underlies clinical aggressiveness and treatment resistance in CLL.

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.

Detecting SARS-CoV-2 lineages and mutational load in municipal wastewater and a use-case in the metropolitan area of Thessaloniki, Greece.

The COVID-19 pandemic represents an unprecedented global crisis necessitating novel approaches for, amongst others, early detection of emerging variants relating to the evolution and spread of the virus. Recently, the detection of SARS-CoV-2 RNA in wastewater has emerged as a useful tool to monitor the prevalence of the virus in the community. Here, we propose a novel methodology, called lineagespot, for the monitoring of mutations and the detection of SARS-CoV-2 lineages in wastewater samples using next-generation sequencing (NGS). Our proposed method was tested and evaluated using NGS data produced by the sequencing of 14 wastewater samples from the municipality of Thessaloniki, Greece, covering a 6-month period. The results showed the presence of SARS-CoV-2 variants in wastewater data. lineagespot was able to record the evolution and rapid domination of the Alpha variant (B.1.1.7) in the community, and allowed the correlation between the mutations evident through our approach and the mutations observed in patients from the same area and time periods. lineagespot is an open-source tool, implemented in R, and is freely available on GitHub and registered on bio.tools.