Tracking SARS-CoV-2 variants of concern in wastewater: an assessment of nine computational tools using simulated genomic data.

Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta-Omicron recombinant and a synthetic 'novel' lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

IκBε deficiency accelerates disease development in chronic lymphocytic leukemia.

The NFKBIE gene, which encodes the NF-κB inhibitor IκBε, is mutated in 3-7% of patients with chronic lymphocytic leukemia (CLL). The most recurrent alteration is a 4-bp frameshift deletion associated with NF-κB activation in leukemic B cells and poor clinical outcome. To study the functional consequences of NFKBIE gene inactivation, both in vitro and in vivo, we engineered CLL B cells and CLL-prone mice to stably down-regulate NFKBIE expression and investigated its role in controlling NF-κB activity and disease expansion. We found that IκBε loss leads to NF-κB pathway activation and promotes both migration and proliferation of CLL cells in a dose-dependent manner. Importantly, NFKBIE inactivation was sufficient to induce a more rapid expansion of the CLL clone in lymphoid organs and contributed to the development of an aggressive disease with a shortened survival in both xenografts and genetically modified mice. IκBε deficiency was associated with an alteration of the MAPK pathway, also confirmed by RNA-sequencing in NFKBIE-mutated patient samples, and resistance to the BTK inhibitor ibrutinib. In summary, our work underscores the multimodal relevance of the NF-κB pathway in CLL and paves the way to translate these findings into novel therapeutic options.

Intranuclear inclusions of polyQ-expanded ATXN1 sequester RNA molecules.

Spinocerebellar ataxia type 1 (SCA1) is an autosomal dominant neurodegenerative disease caused by a trinucleotide (CAG) repeat expansion in the ATXN1 gene. It is characterized by the presence of polyglutamine (polyQ) intranuclear inclusion bodies (IIBs) within affected neurons. In order to investigate the impact of polyQ IIBs in SCA1 pathogenesis, we generated a novel protein aggregation model by inducible overexpression of the mutant ATXN1(Q82) isoform in human neuroblastoma SH-SY5Y cells. Moreover, we developed a simple and reproducible protocol for the efficient isolation of insoluble IIBs. Biophysical characterization showed that polyQ IIBs are enriched in RNA molecules which were further identified by next-generation sequencing. Finally, a protein interaction network analysis indicated that sequestration of essential RNA transcripts within ATXN1(Q82) IIBs may affect the ribosome resulting in error-prone protein synthesis and global proteome instability. These findings provide novel insights into the molecular pathogenesis of SCA1, highlighting the role of polyQ IIBs and their impact on critical cellular processes.

High risk HPV-positive women cervicovaginal microbial profiles in a Greek cohort: a retrospective analysis of the GRECOSELF study.

Increasing evidence supports a role for the vaginal microbiome (VM) in the severity of HPV infection and its potential link to cervical intraepithelial neoplasia. However, a lot remains unclear regarding the precise role of certain bacteria in the context of HPV positivity and persistence of infection. Here, using next generation sequencing (NGS), we comprehensively profiled the VM in a series of 877 women who tested positive for at least one high risk HPV (hrHPV) type with the COBAS® 4,800 assay, after self-collection of a cervico-vaginal sample. Starting from gDNA, we PCR amplified the V3-V4 region of the bacterial 16S rRNA gene and applied a paired-end NGS protocol (Illumina). We report significant differences in the abundance of certain bacteria compared among different HPV-types, more particularly concerning species assigned to Lacticaseibacillus, Megasphaera and Sneathia genera. Especially for Lacticaseibacillus, we observed significant depletion in the case of HPV16, HPV18 versus hrHPVother. Overall, our results suggest that the presence or absence of specific cervicovaginal microbial genera may be linked to the observed severity in hrHPV infection, particularly in the case of HPV16, 18 types.

Software pipelines for RNA-Seq, ChIP-Seq and germline variant calling analyses in common workflow language (CWL).

Background: Automating data analysis pipelines is a key requirement to ensure reproducibility of results, especially when dealing with large volumes of data. Here we assembled automated pipelines for the analysis of High-throughput Sequencing (HTS) data originating from RNA-Seq, ChIP-Seq and Germline variant calling experiments. We implemented these workflows in Common workflow language (CWL) and evaluated their performance by: i) reproducing the results of two previously published studies on Chronic Lymphocytic Leukemia (CLL), and ii) analyzing whole genome sequencing data from four Genome in a Bottle Consortium (GIAB) samples, comparing the detected variants against their respective golden standard truth sets. Findings: We demonstrated that CWL-implemented workflows clearly achieved high accuracy in reproducing previously published results, discovering significant biomarkers and detecting germline SNP and small INDEL variants. Conclusion: CWL pipelines are characterized by reproducibility and reusability; combined with containerization, they provide the ability to overcome issues of software incompatibility and laborious configuration requirements. In addition, they are flexible and can be used immediately or adapted to the specific needs of an experiment or study. The CWL-based workflows developed in this study, along with version information for all software tools, are publicly available on GitHub (https://github.com/BiodataAnalysisGroup/CWL_HTS_pipelines) under the MIT License. They are suitable for the analysis of short-read (such as Illumina-based) data and constitute an open resource that can facilitate automation, reproducibility and cross-platform compatibility for standard bioinformatic analyses.

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.

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.

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.

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.

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.

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

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.

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.

miRkit: R framework analyzing miRNA PCR array data.

OBJECTIVE: The characterization of microRNAs (miRNA) in recent years is an important advance in the field of gene regulation. To this end, several approaches for miRNA expression analysis and various bioinformatics tools have been developed over the last few years. It is a common practice to analyze miRNA PCR Array data using the commercially available software, mostly due to its convenience and ease-of-use. RESULTS: In this work we present miRkit, an open source framework written in R, that allows for the comprehensive analysis of RT-PCR data, from the processing of raw data to a functional analysis of the produced results. The main goal of the proposed tool is to provide an assessment of the samples' quality, perform data normalization by endogenous and exogenous miRNAs, and facilitate differential and functional enrichment analysis. The tool offers fast execution times with low memory usage, and is freely available under a ΜΙΤ license from https://bio.tools/mirkit . Overall, miRkit offers the full analysis from the raw RT-PCR data to functional analysis of targeted genes, and specifically designed to support the popular miScript miRNA PCR Array (Qiagen) technology.

Exploring plant diversity through soil DNA in Thai national parks for influencing land reform and agriculture planning.

BACKGROUND: The severe deforestation, as indicated in national forest data, is a recurring problem in many areas of Northern Thailand, including Doi Suthep-Pui National Park. Agricultural expansion in these areas, is one of the major drivers of deforestation, having adverse consequences on local plant biodiversity. Conserving biodiversity is mainly dependent on the biological monitoring of species distribution and population sizes. However, the existing conventional approaches for monitoring biodiversity are rather limited. METHODS: Here, we explored soil DNA at four forest types in Doi Suthep-Pui National Park in Northern Thailand. Three soil samples, composed of different soil cores mixed together, per sampling location were collected. Soil biodiversity was investigated through eDNA metabarcoding analysis using primers targeting the P6 loop of the plastid DNA trnL (UAA) intron. RESULTS: The distribution of taxa for each sample was found to be similar between replicates. A strong congruence between the conventional morphology- and eDNA-based data of plant diversity in the studied areas was observed. All species recorded by conventional survey with DNA data deposited in the GenBank were detected through the eDNA analysis. Moreover, traces of crops, such as lettuce, maize, wheat and soybean, which were not expected and were not visually detected in the forest area, were identified. It is noteworthy that neighboring land and areas in the studied National Park were once used for crop cultivation, and even to date there is still agricultural land within a 5-10 km radius from the forest sites where the soil samples were collected. The presence of cultivated area near the forest may suggest that we are now facing agricultural intensification leading to deforestation. Land reform for agriculture usage necessitates coordinated planning in order to preserve the forest area. In that context, the eDNA-based data would be useful for influencing policies and management towards this goal.

A Computational Framework for Pattern Detection on Unaligned Sequences: An Application on SARS-CoV-2 Data.

The exponential growth of genome sequences available has spurred research on pattern detection with the aim of extracting evolutionary signal. Traditional approaches, such as multiple sequence alignment, rely on positional homology in order to reconstruct the phylogenetic history of taxa. Yet, mining information from the plethora of biological data and delineating species on a genetic basis, still proves to be an extremely difficult problem to consider. Multiple algorithms and techniques have been developed in order to approach the problem multidimensionally. Here, we propose a computational framework for identifying potentially meaningful features based on k-mers retrieved from unaligned sequence data. Specifically, we have developed a process which makes use of unsupervised learning techniques in order to identify characteristic k-mers of the input dataset across a range of different k-values and within a reasonable time frame. We use these k-mers as features for clustering the input sequences and identifying differences between the distributions of k-mers across the dataset. The developed algorithm is part of an innovative and much promising approach both to the problem of grouping sequence data based on their inherent characteristic features, as well as for the study of changes in the distributions of k-mers, as the k-value is fluctuating within a range of values. Our framework is fully developed in Python language as an open source software licensed under the MIT License, and is freely available at https://github.com/BiodataAnalysisGroup/kmerAnalyzer.

UMIc: A Preprocessing Method for UMI Deduplication and Reads Correction.

A recent refinement in high-throughput sequencing involves the incorporation of unique molecular identifiers (UMIs), which are random oligonucleotide barcodes, on the library preparation steps. A UMI adds a unique identity to different DNA/RNA input molecules through polymerase chain reaction (PCR) amplification, thus reducing bias of this step. Here, we propose an alignment free framework serving as a preprocessing step of fastq files, called UMIc, for deduplication and correction of reads building consensus sequences from each UMI. Our approach takes into account the frequency and the Phred quality of nucleotides and the distances between the UMIs and the actual sequences. We have tested the tool using different scenarios of UMI-tagged library data, having in mind the aspect of a wide application. UMIc is an open-source tool implemented in R and is freely available from https://github.com/BiodataAnalysisGroup/UMIc.

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 .