Transcriptome Profile Analyses of Head Kidney in Roach (Rutilus rutilus), Common Bream (Abramis brama) and Their Hybrids: Does Infection by Monogenean Parasites in Freshwater Fish Reveal Differences in Fish Vigour among Parental Species and Their Hybrids?

Hybrid generations usually face either a heterosis advantage or a breakdown, that can be expressed by the level of parasite infection in hybrid hosts. Hybrids are less infected by parasites than parental species (especially F1 generations) or more infected than parental species (especially post-F1 generations). We performed the experiment with blood-feeding gill parasite Paradiplozoon homoion (Monogenea) infecting leuciscid species, Abramis brama and Rutilus rutilus, their F1 generation and two backcross generations. Backcross generations tended to be more parasitized than parental lines and the F1 generation. The number of differentially expressed genes (DEGs) was lower in F1 hybrids and higher in backcross hybrids when compared to each of the parental lines. The main groups of DEGs were shared among lines; however, A. brama and R. rutilus differed in some of the top gene ontology (GO) terms. DEG analyses revealed the role of heme binding and erythrocyte differentiation after infection by blood-feeding P. homoion. Two backcross generations shared some of the top GO terms, representing mostly downregulated genes associated with P. homoion infection. KEGG analysis revealed the importance of disease-associated pathways; the majority of them were shared by two backcross generations. Our study revealed the most pronounced DEGs associated with blood-feeding monogeneans in backcross hybrids, potentially (but not exclusively) explainable by hybrid breakdown. The lower DEGs reported in F1 hybrids being less parasitized than backcross hybrids is in line with the hybrid advantage.

From fish to cells: Establishment of continuous cell lines from embryos of annual killifish Nothobranchius furzeri and N. kadleci.

There is a growing need of alternative experimental models that avoid or minimize the use of animals due to ethical, economical, and scientific reasons. Surprisingly, the stable embryonic cell lines representing Nothobranchius spp., emerging vertebrate models in aging research, regenerative medicine, ecotoxicology, or genomics, have been not derived so far. This paper reports establishment and deep characterization of ten continuous cell lines from annual killifish embryos of N. furzeri and N. kadleci. The established cell lines exhibited mostly fibroblast- and epithelial-like morphology and steady growth rates with cell doubling time ranging from 27 to 40 h. All cell lines retained very similar characteristics even after continuous subcultivation (more than 100 passages) and extended storage in liquid nitrogen (∼3 years). The cytogenetic analysis of the cell lines revealed a diploid chromosome number mostly equal to 38 elements (i.e., the native chromosome count for both killifish species), with minor but diverse line/passage-specific karyotype changes compared to the patterns observed in non-cultured N. furzeri and N. kadleci somatic cells. Based on transcriptional analysis of marker genes, the cell lines displayed features of an undifferentiated state without signs of senescence even in advanced passages. We confirmed that the cell lines are transfectable and can form viable 3-D spheroids. The applicability of the cell lines for (eco)toxicological surveys was confirmed by assessing the effect of cytotoxic and growth inhibitory agents. Properties of established Nothobranchius embryonic cell lines open new possibilities for the application of this model in various fields of life sciences including molecular mechanisms of aging, karyotype (in)stability or differences in lifespan.

Radiogenomic markers enable risk stratification and inference of mutational pathway states in head and neck cancer.

PURPOSE: Head and neck squamous cell carcinomas (HNSCCs) are a molecularly, histologically, and clinically heterogeneous set of tumors originating from the mucosal epithelium of the oral cavity, pharynx, and larynx. This heterogeneous nature of HNSCC is one of the main contributing factors to the lack of prognostic markers for personalized treatment. The aim of this study was to develop and identify multi-omics markers capable of improved risk stratification in this highly heterogeneous patient population. METHODS: In this retrospective study, we approached this issue by establishing radiogenomics markers to identify high-risk individuals in a cohort of 127 HNSCC patients. Hybrid in vivo imaging and whole-exome sequencing were employed to identify quantitative imaging markers as well as genetic markers on pathway-level prognostic in HNSCC. We investigated the deductibility of the prognostic genetic markers using anatomical and metabolic imaging using positron emission tomography combined with computed tomography. Moreover, we used statistical and machine learning modeling to investigate whether a multi-omics approach can be used to derive prognostic markers for HNSCC. RESULTS: Radiogenomic analysis revealed a significant influence of genetic pathway alterations on imaging markers. A highly prognostic radiogenomic marker based on cellular senescence was identified. Furthermore, the radiogenomic biomarkers designed in this study vastly outperformed the prognostic value of markers derived from genetics and imaging alone. CONCLUSION: Using the identified markers, a clinically meaningful stratification of patients is possible, guiding the identification of high-risk patients and potentially aiding in the development of effective targeted therapies.

KMT2C methyltransferase domain regulated INK4A expression suppresses prostate cancer metastasis.

BACKGROUND: Frequent truncation mutations of the histone lysine N-methyltransferase KMT2C have been detected by whole exome sequencing studies in various cancers, including malignancies of the prostate. However, the biological consequences of these alterations in prostate cancer have not yet been elucidated. METHODS: To investigate the functional effects of these mutations, we deleted the C-terminal catalytic core motif of Kmt2c specifically in mouse prostate epithelium. We analysed the effect of Kmt2c SET domain deletion in a Pten-deficient PCa mouse model in vivo and of truncation mutations of KMT2C in a large number of prostate cancer patients. RESULTS: We show here for the first time that impaired KMT2C methyltransferase activity drives proliferation and PIN formation and, when combined with loss of the tumour suppressor PTEN, triggers loss of senescence, metastatic dissemination and dramatically reduces life expectancy. In Kmt2c-mutated tumours we show enrichment of proliferative MYC gene signatures and loss of expression of the cell cycle repressor p16INK4A. In addition, we observe a striking reduction in disease-free survival of patients with KMT2C-mutated prostate cancer. CONCLUSIONS: We identified truncating events of KMT2C as drivers of proliferation and PIN formation. Loss of PTEN and KMT2C in prostate cancer results in loss of senescence, metastatic dissemination and reduced life expectancy. Our data demonstrate the prognostic significance of KMT2C mutation status in prostate cancer patients. Inhibition of the MYC signalling axis may be a viable treatment option for patients with KMT2C truncations and therefore poor prognosis.

Memory B-cell like chronic lymphocytic leukaemia is associated with specific methylation profile of WNT5A promoter and undetectable expression of WNT5A gene.

Genome methylation profiles define naïve-like (n-CLL), memory-like (m-CLL), and intermediate (i-CLL) subsets of chronic lymphocytic leukaemia (CLL). The profiles can be easily determined by the analysis of the five-CpG signature. m-CLL, i-CLL, and n-CLL with the good, intermediate, and poor prognoses, respectively, differ by the somatic hypermutation status of the immunoglobulin heavy chain variable gene (IGHV), a widely used prognostic predictor in CLL. We have previously shown that the expression of WNT5A, encoding a ROR1 ligand, distinguishes patients with the worse outcome within the prognostically favourable IGHV-mutated subgroup. To analyse the mechanisms controlling WNT5A expression, we investigated the methylation status of 54 CpG sites within the WNT5A promoter and its relation to the WNT5A gene expression. In a cohort of 59 CLL patients balanced for combinations of IGHV and WNT5A statuses, we identified three promoter CpG sites whose methylation level correlated with the WNT5A expression within the IGHV-mutated subgroup. Further, we complemented our data with the methylation status of the five-CpG signature. IGHV-mutated/WNT5A-negative and IGHV-mutated/WNT5A-positive cases overlapped with m­CLL and i­CLL methylation subgroups, respectively, while most IGHV­unmutated samples were assigned to n-CLL. Median methylation levels of all the three CpG sites in the WNT5A promoter were lowest in i-CLL. Finally, a detailed analysis of m-CLL and i-CLL showed that undetectable WNT5A expression predicts longer treatment-free survival with higher statistical significance than the classification according to the five-CpG signature. To conclude, a favourable m-CLL subgroup is associated with mutated IGHV and undetectable WNT5A expression due to its promoter methylation.

Histopathology of laboratory-reared Nothobranchius fishes: Mycobacterial infections versus neoplastic lesions.

Short-lived killifishes of the genus Nothobranchius Peters, 1868 (Cyprinodontiformes) are considered promising model organisms for biomedical research on ageing and tumorigenesis. We conducted histopathological analysis of 411 adult individuals from three Nothobranchius species to study details on spontaneous age-related neoplastic lesions. Light microscopy based on H&E and toluidine blue-stained sections revealed (a) non-proliferative liver changes with pronounced vacuolation of hepatocytes; (b) proliferation of kidney haemopoietic tissue contributing to excretory system damage; (c) proliferation of splenic mononuclear haemoblasts accompanied by reduced erythropoiesis; (d) proliferation of mononuclear cell aggregates in the liver parenchyma; and (e) rare occurrence of hepatocellular adenomas. Ziehl-Neelsen (ZN) staining revealed that the proliferative lesions are a host defence response to mycobacterial infections manifested by activation of the mononuclear phagocytic system and atypical granulomatous inflammatory reaction. 16S rRNA analysis identified three species of Mycobacterium in our samples. Our findings turn attention to lesions which mimic neoplasms by their gross appearance and question the light microscopic interpretation of lesions unless differential ZN staining is included. Beyond the limitations of our morphological approach, the intensity of mycobacterial infections is a challenging opportunity for research into the molecular-genetic background of the mononuclear phagocytic system reaction in Nothobranchius killifish.

A Simple RNA Target Capture NGS Strategy for Fusion Genes Assessment in the Diagnostics of Pediatric B-cell Acute Lymphoblastic Leukemia.

Acute lymphoblastic leukemia (ALL) is the most frequent pediatric cancer. Fusion genes are hallmarks of ALL, and they are used as biomarkers for risk stratification as well as targets for precision medicine. Hence, clinical diagnostics pursues broad and comprehensive strategies for accurate discovery of fusion genes. Currently, the gold standard methodologies for fusion gene detection are fluorescence in situ hybridization and polymerase chain reaction; these, however, lack sensitivity for the identification of new fusion genes and breakpoints. In this study, we implemented a simple operating procedure (OP) for detecting fusion genes. The OP employs RNA CaptureSeq, a versatile and effortless next-generation sequencing assay, and an in-house as well as a purpose-built bioinformatics pipeline for the subsequent data analysis. The OP was evaluated on a cohort of 89 B-cell precursor ALL (BCP-ALL) pediatric samples annotated as negative for fusion genes by the standard techniques. The OP confirmed 51 samples as negative for fusion genes, and, more importantly, it identified known (KMT2A rearrangements) as well as new fusion events (JAK2 rearrangements) in the remaining 38 investigated samples, of which 16 fusion genes had prognostic significance. Herein, we describe the OP and its deployment into routine ALL diagnostics, which will allow substantial improvements in both patient risk stratification and precision medicine.

Standardized next-generation sequencing of immunoglobulin and T-cell receptor gene recombinations for MRD marker identification in acute lymphoblastic leukaemia; a EuroClonality-NGS validation study.

Amplicon-based next-generation sequencing (NGS) of immunoglobulin (IG) and T-cell receptor (TR) gene rearrangements for clonality assessment, marker identification and quantification of minimal residual disease (MRD) in lymphoid neoplasms has been the focus of intense research, development and application. However, standardization and validation in a scientifically controlled multicentre setting is still lacking. Therefore, IG/TR assay development and design, including bioinformatics, was performed within the EuroClonality-NGS working group and validated for MRD marker identification in acute lymphoblastic leukaemia (ALL). Five EuroMRD ALL reference laboratories performed IG/TR NGS in 50 diagnostic ALL samples, and compared results with those generated through routine IG/TR Sanger sequencing. A central polytarget quality control (cPT-QC) was used to monitor primer performance, and a central in-tube quality control (cIT-QC) was spiked into each sample as a library-specific quality control and calibrator. NGS identified 259 (average 5.2/sample, range 0-14) clonal sequences vs. Sanger-sequencing 248 (average 5.0/sample, range 0-14). NGS primers covered possible IG/TR rearrangement types more completely compared with local multiplex PCR sets and enabled sequencing of bi-allelic rearrangements and weak PCR products. The cPT-QC showed high reproducibility across all laboratories. These validated and reproducible quality-controlled EuroClonality-NGS assays can be used for standardized NGS-based identification of IG/TR markers in lymphoid malignancies.

Bioinformatic pipelines for whole transcriptome sequencing data exploitation in leukemia patients with complex structural variants.

BACKGROUND: Extensive genome rearrangements, known as chromothripsis, have been recently identified in several cancer types. Chromothripsis leads to complex structural variants (cSVs) causing aberrant gene expression and the formation of de novo fusion genes, which can trigger cancer development, or worsen its clinical course. The functional impact of cSVs can be studied at the RNA level using whole transcriptome sequencing (total RNA-Seq). It represents a powerful tool for discovering, profiling, and quantifying changes of gene expression in the overall genomic context. However, bioinformatic analysis of transcriptomic data, especially in cases with cSVs, is a complex and challenging task, and the development of proper bioinformatic tools for transcriptome studies is necessary. METHODS: We designed a bioinformatic workflow for the analysis of total RNA-Seq data consisting of two separate parts (pipelines): The first pipeline incorporates a statistical solution for differential gene expression analysis in a biologically heterogeneous sample set. We utilized results from transcriptomic arrays which were carried out in parallel to increase the precision of the analysis. The second pipeline is used for the identification of de novo fusion genes. Special attention was given to the filtering of false positives (FPs), which was achieved through consensus fusion calling with several fusion gene callers. We applied the workflow to the data obtained from ten patients with chronic lymphocytic leukemia (CLL) to describe the consequences of their cSVs in detail. The fusion genes identified by our pipeline were correlated with genomic break-points detected by genomic arrays. RESULTS: We set up a novel solution for differential gene expression analysis of individual samples and de novo fusion gene detection from total RNA-Seq data. The results of the differential gene expression analysis were concordant with results obtained by transcriptomic arrays, which demonstrates the analytical capabilities of our method. We also showed that the consensus fusion gene detection approach was able to identify true positives (TPs) efficiently. Detected coordinates of fusion gene junctions were in concordance with genomic breakpoints assessed using genomic arrays. DISCUSSION: Byapplying our methods to real clinical samples, we proved that our approach for total RNA-Seq data analysis generates results consistent with other genomic analytical techniques. The data obtained by our analyses provided clues for the study of the biological consequences of cSVs with far-reaching implications for clinical outcome and management of cancer patients. The bioinformatic workflow is also widely applicable for addressing other research questions in different contexts, for which transcriptomic data are generated.

Quality control and quantification in IG/TR next-generation sequencing marker identification: protocols and bioinformatic functionalities by EuroClonality-NGS.

Assessment of clonality, marker identification and measurement of minimal residual disease (MRD) of immunoglobulin (IG) and T cell receptor (TR) gene rearrangements in lymphoid neoplasms using next-generation sequencing (NGS) is currently under intensive development for use in clinical diagnostics. So far, however, there is a lack of suitable quality control (QC) options with regard to standardisation and quality metrics to ensure robust clinical application of such approaches. The EuroClonality-NGS Working Group has therefore established two types of QCs to accompany the NGS-based IG/TR assays. First, a central polytarget QC (cPT-QC) is used to monitor the primer performance of each of the EuroClonality multiplex NGS assays; second, a standardised human cell line-based DNA control is spiked into each patient DNA sample to work as a central in-tube QC and calibrator for MRD quantification (cIT-QC). Having integrated those two reference standards in the ARResT/Interrogate bioinformatic platform, EuroClonality-NGS provides a complete protocol for standardised IG/TR gene rearrangement analysis by NGS with high reproducibility, accuracy and precision for valid marker identification and quantification in diagnostics of lymphoid malignancies.

First evidence of a paediatric patient with Cornelia de Lange syndrome with acute lymphoblastic leukaemia.

Cornelia de Lange syndrome (CdLS) is a rare autosomal-dominant genetic disorder characterised by prenatal and postnatal growth and mental retardation, facial dysmorphism and upper limb abnormalities. Germline mutations of cohesin complex genes SMC1A, SMC3, RAD21 or their regulators NIPBL and HDAC8 have been identified in CdLS as well as somatic mutations in myeloid disorders. We describe the first case of a paediatric patient with CdLS with B-cell precursor Acute Lymphoblastic Leukaemia (ALL). The patient did not show any unusual cytogenetic abnormality, and he was enrolled into the high risk arm of AIEOP-BFM ALL2009 protocol because of slow early response, but 3 years after discontinuation, he experienced an ALL relapse. We identified a heterozygous mutation in exon 46 of NIPBL, causing frameshift and a premature stop codon (RNA-Targeted Next generation Sequencing Analysis). The analysis of the family indicated a de novo origin of this previously not reported deleterious variant. As for somatic cohesin mutations in acute myeloid leukaemia, also this ALL case was not affected by aneuploidy, thus suggesting a major impact of the non-canonical role of NIPBL in gene regulation. A potential biological role of NIPBL in leukaemia has still to be dissected.

Activation-induced deaminase and its splice variants associate with trisomy 12 in chronic lymphocytic leukemia.

Activation-induced cytidine deaminase (AID) is a mutator enzyme essential for somatic hypermutation (SHM) and class switch recombination (CSR) during effective adaptive immune responses. Its aberrant expression and activity have been detected in lymphomas, leukemias, and solid tumors. In chronic lymphocytic leukemia (CLL) increased expression of alternatively spliced AID variants has been documented. We used real-time RT-PCR to quantify the expression of AID and its alternatively spliced transcripts (AIDΔE4a, AIDΔE4, AIDivs3, and AIDΔE3E4) in 149 CLL patients and correlated this expression to prognostic markers including recurrent chromosomal aberrations, the presence of complex karyotype, mutation status of the immunoglobulin heavy chain variable gene, and recurrent mutations. We report a previously unappreciated association between higher AID transcript levels and trisomy of chromosome 12. Functional analysis of AID splice variants revealed loss of their activity with respect to SHM, CSR, and induction of double-strand DNA breaks. In silico modeling provided insight into the molecular interactions and structural dynamics of wild-type AID and a shortened AID variant closely resembling AIDΔE4, confirming its loss-of-function phenotype.

ToTem: a tool for variant calling pipeline optimization.

BACKGROUND: High-throughput bioinformatics analyses of next generation sequencing (NGS) data often require challenging pipeline optimization. The key problem is choosing appropriate tools and selecting the best parameters for optimal precision and recall. RESULTS: Here we introduce ToTem, a tool for automated pipeline optimization. ToTem is a stand-alone web application with a comprehensive graphical user interface (GUI). ToTem is written in Java and PHP with an underlying connection to a MySQL database. Its primary role is to automatically generate, execute and benchmark different variant calling pipeline settings. Our tool allows an analysis to be started from any level of the process and with the possibility of plugging almost any tool or code. To prevent an over-fitting of pipeline parameters, ToTem ensures the reproducibility of these by using cross validation techniques that penalize the final precision, recall and F-measure. The results are interpreted as interactive graphs and tables allowing an optimal pipeline to be selected, based on the user's priorities. Using ToTem, we were able to optimize somatic variant calling from ultra-deep targeted gene sequencing (TGS) data and germline variant detection in whole genome sequencing (WGS) data. CONCLUSIONS: ToTem is a tool for automated pipeline optimization which is freely available as a web application at https://totem.software .

The influence of mutational status and biological characteristics of acute myeloid leukemia on xenotransplantation outcomes in NOD SCID gamma mice.

PURPOSE: This study aimed at analyzing the association of gene mutations and other acute myeloid leukemia (AML) characteristics with engraftment outcomes in immunodeficient mice and to select the engraftment outcomes that best reflect patient survival. METHODS: Mutations in 19 genes as well as leukemia- and patient-related characteristics were analyzed for a group of 47 de novo AML samples with respect to three engraftment outcomes: engraftment ability, engraftment intensity (percentage of hCD45+ cells) and engraftment latency. Leukemia-related characteristics were additionally analyzed in an extended group of 68 samples that included the 47 de novo samples, and additional 21 samples from refractory and relapsed cases. Engraftment outcomes were compared with overall and event-free survival of the patients. RESULTS: For the 47 de novo samples, no single mutation influenced engraftment, whereas the NPM1 mut /DNMT3A mut co-mutation was associated with higher engraftment ability. NPM1 mut /FLT3-ITD neg had lower engraftment intensity. Among leukemia-related characteristics, a complex karyotype was associated with higher engraftment intensity. Among patient-related characteristics, higher cytogenetic risk was associated with higher engraftment intensity, and failure to achieve clinical remission was associated with shorter engraftment latency. In the extended group of 68 samples, white blood count was associated with higher engraftment ability, and the presence of a complex karyotype was associated with higher engraftment intensity. Association with patient overall survival was seen only for engraftment intensity. CONCLUSIONS: The engraftment of AML was influenced by mutation-interactions and other AML characteristics, rather than by single mutated genes, and engraftment intensity best reflected clinical penetrance of AML.

GLASS: assisted and standardized assessment of gene variations from Sanger sequence trace data.

MOTIVATION: Sanger sequencing is still being employed for sequence variant detection by many laboratories, especially in a clinical setting. However, chromatogram interpretation often requires manual inspection and in some cases, considerable expertise. RESULTS: We present GLASS, a web-based Sanger sequence trace viewer, editor, aligner and variant caller, built to assist with the assessment of variations in 'curated' or user-provided genes. Critically, it produces a standardized variant output as recommended by the Human Genome Variation Society. AVAILABILITY AND IMPLEMENTATION: GLASS is freely available at http://bat.infspire.org/genomepd/glass/ with source code at https://github.com/infspiredBAT/GLASS. CONTACT: nikos.darzentas@gmail.com or malcikova.jitka@fnbrno.cz. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

High resolution IgH repertoire analysis reveals fetal liver as the likely origin of life-long, innate B lymphopoiesis in humans.

The ontogeny of the natural, public IgM repertoire remains incompletely explored. Here, high-resolution immunogenetic analysis of B cells from (unrelated) fetal, child, and adult samples, shows that although fetal liver (FL) and bone marrow (FBM) IgM repertoires are equally diversified, FL is the main source of IgM natural immunity during the 2nd trimester. Strikingly, 0.25% of all prenatal clonotypes, comprising 18.7% of the expressed repertoire, are shared with the postnatal samples, consistent with persisting fetal IgM+ B cells being a source of natural IgM repertoire in adult life. Further, the origins of specific stereotypic IgM+ B cell receptors associated with chronic lymphocytic leukemia, can be traced back to fetal B cell lymphopoiesis, suggesting that persisting fetal B cells can be subject to malignant transformation late in life. Overall, these novel data provide unique insights into the ontogeny of physiological and malignant B lymphopoiesis that spans the human lifetime.

ARResT/Interrogate: an interactive immunoprofiler for IG/TR NGS data.

Motivation: The study of immunoglobulins and T cell receptors using next-generation sequencing has finally allowed exploring immune repertoires and responses in their immense variability and complexity. Unsurprisingly, their analysis and interpretation is a highly convoluted task. Results: We thus implemented ARResT/Interrogate, a web-based, interactive application. It can organize and filter large amounts of immunogenetic data by numerous criteria, calculate several relevant statistics, and present results in the form of multiple interconnected visualizations. Availability and Implementation: ARResT/Interrogate is implemented primarily in R, and is freely available at http://bat.infspire.org/arrest/interrogate/ Contact: nikos.darzentas@gmail.com Supplementary Information: Supplementary data are available at Bioinformatics online.

A role for palindromic structures in the cis-region of maize Sirevirus LTRs in transposable element evolution and host epigenetic response.

Transposable elements (TEs) proliferate within the genome of their host, which responds by silencing them epigenetically. Much is known about the mechanisms of silencing in plants, particularly the role of siRNAs in guiding DNA methylation. In contrast, little is known about siRNA targeting patterns along the length of TEs, yet this information may provide crucial insights into the dynamics between hosts and TEs. By focusing on 6456 carefully annotated, full-length Sirevirus LTR retrotransposons in maize, we show that their silencing associates with underlying characteristics of the TE sequence and also uncover three features of the host-TE interaction. First, siRNA mapping varies among families and among elements, but particularly along the length of elements. Within the cis-regulatory portion of the LTRs, a complex palindrome-rich region acts as a hotspot of both siRNA matching and sequence evolution. These patterns are consistent across leaf, tassel, and immature ear libraries, but particularly emphasized for floral tissues and 21- to 22-nt siRNAs. Second, this region has the ability to form hairpins, making it a potential template for the production of miRNA-like, hairpin-derived small RNAs. Third, Sireviruses are targeted by siRNAs as a decreasing function of their age, but the oldest elements remain highly targeted, partially by siRNAs that cross-map to the youngest elements. We show that the targeting of older Sireviruses reflects their conserved palindromes. Altogether, we hypothesize that the palindromes aid the silencing of active elements and influence transposition potential, siRNA targeting levels, and ultimately the fate of an element within the genome.

ARResT/AssignSubsets: a novel application for robust subclassification of chronic lymphocytic leukemia based on B cell receptor IG stereotypy.

MOTIVATION: An ever-increasing body of evidence supports the importance of B cell receptor immunoglobulin (BcR IG) sequence restriction, alias stereotypy, in chronic lymphocytic leukemia (CLL). This phenomenon accounts for ∼30% of studied cases, one in eight of which belong to major subsets, and extends beyond restricted sequence patterns to shared biologic and clinical characteristics and, generally, outcome. Thus, the robust assignment of new cases to major CLL subsets is a critical, and yet unmet, requirement. RESULTS: We introduce a novel application, ARResT/AssignSubsets, which enables the robust assignment of BcR IG sequences from CLL patients to major stereotyped subsets. ARResT/AssignSubsets uniquely combines expert immunogenetic sequence annotation from IMGT/V-QUEST with curation to safeguard quality, statistical modeling of sequence features from more than 7500 CLL patients, and results from multiple perspectives to allow for both objective and subjective assessment. We validated our approach on the learning set, and evaluated its real-world applicability on a new representative dataset comprising 459 sequences from a single institution. AVAILABILITY AND IMPLEMENTATION: ARResT/AssignSubsets is freely available on the web at http://bat.infspire.org/arrest/assignsubsets/ CONTACT: nikos.darzentas@gmail.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.