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

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.</p

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

Supporting data for "Software pipelines for RNA-Seq, ChIP-Seq and Germline Variant calling analyses in Common Workflow Language (CWL)"

Datasets produced during the validation of CWL-based pipelines, designed for the analysis of data from RNA-Seq, ChIP-Seq and germline variant calling experiments. Specifically, the workflows were tested using publicly available High-throughput (HTS) data from published studies on Chronic Lymphocytic Leukemia (CLL) (accession numbers: E-MTAB-6962, GSE115772) and Genome in a Bottle (GIAB) project samples (accession numbers: SRR6794144, SRR22476789, SRR22476790, SRR22476791). The supporting data include: Differential transcript and gene expression results produced during the analysis with the CWL-based RNA-Seq pipeline Bigwig and narrowPeak files, differential binding results, table of consensus peaks and read counts of EZH2 and H3K27me3, produced during the analysis with the CWL-based ChIP-Seq pipeline VCF files containing the detected and filtered variants, along with the respective hap.py () results regarding comparisons against the GIAB golden standard truth sets for both CWL-based germline variant calling pipeline

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.‘Hellenic Network for Precision Medicine’ in the framework of the Hellenic Republic – Siemens Settlement Agreement; European Regional Development Fund of the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship and Innovation, under the call RESEARCH – CREATE – INNOVATE [GenOptics, project code: T2E1DK-00407

Blood Metabolomics May Discriminate a Sub-Group of Patients with First Demyelinating Episode in the Context of RRMS with Increased Disability and MRI Characteristics Indicative of Poor Prognosis

Biomarker research across the health-to-disease continuum is being increasingly applied. We applied blood-based metabolomics in order to identify patient clusters with a first demyelinating episode, and explored the prognostic potential of the method by thoroughly characterizing each cluster in terms of clinical, laboratory and MRI markers of established prognostic potential for Multiple Sclerosis (MS). Recruitment consisted of 11 patients with Clinically Isolated Syndrome (CIS), 37 patients with a first demyelinating episode in the context of Relapsing-Remitting MS (RRMS) and 11 control participants. Blood-based metabolomics and hierarchical clustering analysis (HCL) were applied. Constructed OPLS-DA models illustrated a discrimination between patients with CIS and the controls (p = 0.0014), as well as between patients with RRMS and the controls (p = 1 × 10−5). Hierarchical clustering analysis (HCL) for patients with RRMS identified three clusters. RRMS-patients-cluster-3 exhibited higher mean cell numbers in the Cerebro-spinal Fluid (CSF) compared to patients with CIS (18.17 ± 6.3 vs. 1.09 ± 0.41, p = 0.004). Mean glucose CSF/serum ratio and infratentorial lesion burden significantly differed across CIS- and HCL-derived RRMS-patient clusters (F = 14.95, p p = 0.002, respectively), mainly due to increased mean values for patients with RRMS-cluster-3. HCL discriminated a cluster of patients with a first demyelinating episode in the context of RRMS with increased disability, laboratory findings linked with increased pathology burden and MRI markers of poor prognosis

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

The chromatin activation landscape of chronic lymphocytic leukemia (CLL) with stereo-typed 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 non-stereotyped CLL counterparts, subset 8 displayed a remarkably distinct chromatin acti-vation 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

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

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.</p

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

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

BackgroundMicroenvironmental 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 designTR 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.ResultsOligoclonal 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.ConclusionsDistinct 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.</p