Diagnosing Viral Infections Through T-Cell Receptor Sequencing of Activated CD8+ T Cells.

T-cell-based diagnostic tools identify pathogen exposure but lack differentiation between recent and historical exposures in acute infectious diseases. Here, T-cell receptor (TCR) RNA sequencing was performed on HLA-DR+/CD38+CD8+ T-cell subsets of hospitalized coronavirus disease 2019 (COVID-19) patients (n = 30) and healthy controls (n = 30; 10 of whom had previously been exposed to severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]). CDR3α and CDR3ß TCR regions were clustered separately before epitope specificity annotation using a database of SARS-CoV-2-associated CDR3α and CDR3ß sequences corresponding to >1000 SARS-CoV-2 epitopes. The depth of the SARS-CoV-2-associated CDR3α/ß sequences differentiated COVID-19 patients from the healthy controls with a receiver operating characteristic area under the curve of 0.84 ± 0.10. Hence, annotating TCR sequences of activated CD8+ T cells can be used to diagnose an acute viral infection and discriminate it from historical exposure. In essence, this work presents a new paradigm for applying the T-cell repertoire to accomplish TCR-based diagnostics.

Unravelling the immune signature of herpes zoster: Insights into pathophysiology and the HLA risk profile.

The varicella-zoster virus (VZV) infects over 95% of the population. VZV reactivation causes herpes zoster (HZ), known as shingles, primarily affecting the elderly and immunocompromised individuals. However, HZ can also occur in otherwise healthy individuals. We analyzed the immune signature and risk profile in HZ patients using a genome-wide association study across different UK Biobank HZ cohorts. Additionally, we conducted one of the largest HZ HLA association studies to date, coupled with transcriptomic analysis of pathways underlying HZ susceptibility. Our findings highlight the significance of the MHC locus for HZ development, identifying five protective and four risk HLA alleles. This demonstrates that HZ susceptibility is largely governed by variations in the MHC. Furthermore, functional analyses revealed the upregulation of type I interferon and adaptive immune responses. These findings provide fresh molecular insights into the pathophysiology and the activation of innate and adaptive immune responses triggered by symptomatic VZV reactivation.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania.

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Clinical Decision Support for Improved Neonatal Care: The Development of a Machine Learning Model for the Prediction of Late-onset Sepsis and Necrotizing Enterocolitis.

OBJECTIVE: To develop an artificial intelligence-based software system for predicting late-onset sepsis (LOS) and necrotizing enterocolitis (NEC) in infants admitted to the neonatal intensive care unit (NICU). STUDY DESIGN: Single-center, retrospective cohort study, conducted in the NICU of the Antwerp University Hospital. Continuous monitoring data of 865 preterm infants born at <32 weeks gestational age, admitted to the NICU in the first week of life, were used to train an XGBoost machine learning (ML) algorithm for LOS and NEC prediction in a cross-validated setup. Afterward, the model's performance was assessed on an independent test set of 148 patients (internal validation). RESULTS: The ML model delivered hourly risk predictions with an overall sensitivity of 69% (142/206) for all LOS/NEC episodes and 81% (67/83) for severe LOS/NEC episodes. The model showed a median time gain of ≤10 hours (IQR, 3.1-21.0 hours), compared with historical clinical diagnosis. On the complete retrospective dataset, the ML model made 721 069 predictions, of which 9805 (1.3%) depicted a LOS/NEC probability of ≥0.15, resulting in a total alarm rate of <1 patient alarm-day per week. The model reached a similar performance on the internal validation set. CONCLUSIONS: Artificial intelligence technology can assist clinicians in the early detection of LOS and NEC in the NICU, which potentially can result in clinical and socioeconomic benefits. Additional studies are required to quantify further the effect of combining artificial and human intelligence on patient outcomes in the NICU.

Current challenges for unseen-epitope TCR interaction prediction and a new perspective derived from image classification.

The prediction of epitope recognition by T-cell receptors (TCRs) has seen many advancements in recent years, with several methods now available that can predict recognition for a specific set of epitopes. However, the generic case of evaluating all possible TCR-epitope pairs remains challenging, mainly due to the high diversity of the interacting sequences and the limited amount of currently available training data. In this work, we provide an overview of the current state of this unsolved problem. First, we examine appropriate validation strategies to accurately assess the generalization performance of generic TCR-epitope recognition models when applied to both seen and unseen epitopes. In addition, we present a novel feature representation approach, which we call ImRex (interaction map recognition). This approach is based on the pairwise combination of physicochemical properties of the individual amino acids in the CDR3 and epitope sequences, which provides a convolutional neural network with the combined representation of both sequences. Lastly, we highlight various challenges that are specific to TCR-epitope data and that can adversely affect model performance. These include the issue of selecting negative data, the imbalanced epitope distribution of curated TCR-epitope datasets and the potential exchangeability of TCR alpha and beta chains. Our results indicate that while extrapolation to unseen epitopes remains a difficult challenge, ImRex makes this feasible for a subset of epitopes that are not too dissimilar from the training data. We show that appropriate feature engineering methods and rigorous benchmark standards are required to create and validate TCR-epitope predictive models.

ClusTCR: a python interface for rapid clustering of large sets of CDR3 sequences with unknown antigen specificity.

MOTIVATION: The T-cell receptor (TCR) determines the specificity of a T-cell towards an epitope. As of yet, the rules for antigen recognition remain largely undetermined. Current methods for grouping TCRs according to their epitope specificity remain limited in performance and scalability. Multiple methodologies have been developed, but all of them fail to efficiently cluster large datasets exceeding 1 million sequences. To account for this limitation, we developed ClusTCR, a rapid TCR clustering alternative that efficiently scales up to millions of CDR3 amino acid sequences, without knowledge about their antigen specificity. RESULTS: Benchmarking comparisons revealed similar accuracy of ClusTCR as compared to other TCR clustering methods, as measured by cluster retention, purity and consistency. ClusTCR offers a drastic improvement in clustering speed, which allows the clustering of millions of TCR sequences in just a few minutes through ultraefficient similarity searching and sequence hashing. AVAILABILITY AND IMPLEMENTATION: ClusTCR was written in Python 3. It is available as an anaconda package (https://anaconda.org/svalkiers/clustcr) and on github (https://github.com/svalkiers/clusTCR). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sorted B cell transcriptomes point towards actively regulated B cell responses during ongoing chronic hepatitis B infections.

The natural course of chronic hepatitis B virus (HBV) infections follows distinct clinical disease phases, characterized by fluctuating levels of serum HBV DNA and ALT. The immune cells and their features that govern these clinical disease transitions remain unknown. In the current study, we performed RNA sequencing on purified B cells from blood (n = 42) and liver (n = 10) of healthy controls and chronic HBV patients. We found distinct gene expression profiles between healthy controls and chronic HBV patients, as evidenced by 190 differentially expressed genes (DEG), but also between the clinical phenotypes of a chronic HBV infection (17-110 DEG between each phase). Numerous immune pathways, including the B cell receptor pathway were upregulated in liver B cells when compared to peripheral B cells. Further investigation of the detected DEG suggested an activation of B cells during HBeAg seroconversion and an active regulation of B cell signalling in the liver.

On the viability of unsupervised T-cell receptor sequence clustering for epitope preference.

MOTIVATION: The T-cell receptor (TCR) is responsible for recognizing epitopes presented on cell surfaces. Linking TCR sequences to their ability to target specific epitopes is currently an unsolved problem, yet one of great interest. Indeed, it is currently unknown how dissimilar TCR sequences can be before they no longer bind the same epitope. This question is confounded by the fact that there are many ways to define the similarity between two TCR sequences. Here we investigate both issues in the context of TCR sequence unsupervised clustering. RESULTS: We provide an overview of the performance of various distance metrics on two large independent datasets with 412 and 2835 TCR sequences respectively. Our results confirm the presence of structural distinct TCR groups that target identical epitopes. In addition, we put forward several recommendations to perform unsupervised T-cell receptor sequence clustering. AVAILABILITY AND IMPLEMENTATION: Source code implemented in Python 3 available at https://github.com/pmeysman/TCRclusteringPaper. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Oocyte maturation under lipotoxic conditions induces carryover transcriptomic and functional alterations during post-hatching development of good-quality blastocysts: novel insights from a bovine embryo-transfer model.

STUDY QUESTION: Does oocyte maturation under lipolytic conditions have detrimental carry-over effects on post-hatching embryo development of good-quality blastocysts after transfer? SUMMARY ANSWER: Surviving, morphologically normal blastocysts derived from bovine oocytes that matured under lipotoxic conditions exhibit long-lasting cellular dysfunction at the transcriptomic and metabolic levels, which coincides with retarded post-hatching embryo development. WHAT IS KNOWN ALREADY: There is increasing evidence showing that following maturation in pathophysiologically relevant lipotoxic conditions (as in obesity or metabolic syndrome), surviving blastocysts of good (transferable) morphological quality have persistent transcriptomic and epigenetic alteration even when in vitro embryo culture takes place under standard conditions. However, very little is known about subsequent development in the uterus after transfer. STUDY DESIGN, SIZE, DURATION: Bovine oocytes were matured in vitro in the presence of pathophysiologically relevant, high non-esterified fatty acid (NEFA) concentrations (HIGH PA), or in basal NEFA concentrations (BASAL) as a physiological control. Eight healthy multiparous non-lactating Holstein cows were used for embryo transfers. Good-quality blastocysts (pools of eight) were transferred per cow, and cows were crossed over for treatments in the next replicate. Embryos were recovered 7 days later and assessed for post-hatching development, phenotypic features and gene expression profile. Blastocysts from solvent-free and NEFA-free maturation (CONTROL) were also tested for comparison. PARTICIPANTS/MATERIALS, SETTING, METHODS: Recovered Day 14 embryos were morphologically assessed and dissected into embryonic disk (ED) and extraembryonic tissue (EXT). Samples of EXT were cultured for 24 h to assess cellular metabolic activity (glucose and pyruvate consumption and lactate production) and embryos' ability to signal for maternal recognition of pregnancy (interferon-τ secretion; IFN-τ). ED and EXT samples were subjected to RNA sequencing to evaluate the genome-wide transcriptome patterns. MAIN RESULTS AND THE ROLE OF CHANCE: The embryo recovery rate at Day 14 p.i. was not significantly different among treatment groups (P > 0.1). However, higher proportions of HIGH PA embryos were retarded in growth (in spherical stage) compared to the more elongated tubular stage embryos in the BASAL group (P < 0.05). Focusing on the normally developed tubular embryos in both groups, HIGH PA exposure resulted in altered cellular metabolism and altered transcriptome profile particularly in pathways related to redox-regulating mechanisms, apoptosis, cellular growth, interaction and differentiation, energy metabolism and epigenetic mechanisms, compared to BASAL embryos. Maturation under BASAL conditions did not have any significant effects on post-hatching development and cellular functions compared to CONTROL. LARGE-SCALE DATA: The datasets of RNA sequencing analysis are available in the NCBI's Gene Expression Omnibus (GEO) repository, series accession number GSE127889 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE127889). Datasets of differentially expressed genes and their gene ontology functions are available in the Mendeley datasets at http://dx.doi.org/10.17632/my2z7dvk9j.2. LIMITATIONS, REASONS FOR CAUTION: The bovine model was used here to allow non-invasive embryo transfer and post-hatching recovery on Day 14. There are physiological differences in some characteristics of post-hatching embryo development between human and cows, such as embryo elongation and trophoblastic invasion. However, the main carry-over effects of oocyte maturation under lipolytic conditions described here are evident at the cellular level and therefore may also occur during post-hatching development in other species including humans. In addition, post-hatching development was studied here under a healthy uterine environment to focus on carry-over effects originating from the oocyte, whereas additional detrimental effects may be induced by maternal metabolic disorders due to adverse changes in the uterine microenvironment. RNA sequencing results were not verified by qPCR, and no solvent control was included. WIDER IMPLICATIONS OF THE FINDINGS: Our observations may increase the awareness of the importance of maternal metabolic stress at the level of the preovulatory oocyte in relation to carry-over effects that may persist in the transferrable embryos. It should further stimulate new research about preventive and protective strategies to optimize maternal metabolic health around conception to maximize embryo viability and thus fertility outcome. STUDY FUNDING/COMPETING INTEREST(S): This study was supported by the Flemish Research Fund (FWO grant 11L8716N and FWO project 42/FAO10300/6541). The authors declare there are no conflicts of interest.

The use of plasma donor-derived, cell-free DNA to monitor acute rejection after kidney transplantation.

BACKGROUND: After transplantation, cell-free deoxyribonucleic acid (DNA) derived from the donor organ (ddcfDNA) can be detected in the recipient's circulation. We aimed to investigate the role of plasma ddcfDNA as biomarker for acute kidney rejection. METHODS: From 107 kidney transplant recipients, plasma samples were collected longitudinally after transplantation (Day 1 to 3 months) within a multicentre set-up. Cell-free DNA from the donor was quantified in plasma as a fraction of the total cell-free DNA by next generation sequencing using a targeted, multiplex polymerase chain reaction-based method for the analysis of single nucleotide polymorphisms. RESULTS: Increases of the ddcfDNA% above a threshold value of 0.88% were significantly associated with the occurrence of episodes of acute rejection (P = 0.017), acute tubular necrosis (P = 0.011) and acute pyelonephritis (P = 0.032). A receiver operating characteristic curve analysis revealed an equal area under the curve of the ddcfDNA% and serum creatinine of 0.64 for the diagnosis of acute rejection. CONCLUSIONS: Although increases in plasma ddcfDNA% are associated with graft injury, plasma ddcfDNA does not outperform the diagnostic capacity of the serum creatinine in the diagnosis of acute rejection.

Memory CD4+ T cell receptor repertoire data mining as a tool for identifying cytomegalovirus serostatus.

Pathogens of past and current infections have been identified directly by means of PCR or indirectly by measuring a specific immune response (e.g., antibody titration). Using a novel approach, Emerson and colleagues showed that the cytomegalovirus serostatus can also be accurately determined by using a T cell receptor repertoire data mining approach. In this study, we have sequenced the CD4+ memory T cell receptor repertoire of a Belgian cohort with known cytomegalovirus serostatus. A random forest classifier was trained on the CMV specific T cell receptor repertoire signature and used to classify individuals in the Belgian cohort. This study shows that the novel approach can be reliably replicated with an equivalent performance as that reported by Emerson and colleagues. Additionally, it provides evidence that the T cell receptor repertoire signature is to a large extent present in the CD4+ memory repertoire.

Diagnosing enterovirus meningitis via blood transcriptomics: an alternative for lumbar puncture?

BACKGROUND: Meningitis can be caused by several viruses and bacteria. Identifying the causative pathogen as quickly as possible is crucial to initiate the most optimal therapy, as acute bacterial meningitis is associated with a significant morbidity and mortality. Bacterial meningitis requires antibiotics, as opposed to enteroviral meningitis, which only requires supportive therapy. Clinical presentation is usually not sufficient to differentiate between viral and bacterial meningitis, thereby necessitating cerebrospinal fluid (CSF) analysis by PCR and/or time-consuming bacterial cultures. However, collecting CSF in children is not always feasible and a rather invasive procedure. METHODS: In 12 Belgian hospitals, we obtained acute blood samples from children with signs of meningitis (49 viral and 7 bacterial cases) (aged between 3 months and 16 years). After pathogen confirmation on CSF, the patient was asked to give a convalescent sample after recovery. 3' mRNA sequencing was performed to determine differentially expressed genes (DEGs) to create a host transcriptomic profile. RESULTS: Enteroviral meningitis cases displayed the largest upregulated fold change enrichment in type I interferon production, response and signaling pathways. Patients with bacterial meningitis showed a significant upregulation of genes related to macrophage and neutrophil activation. We found several significantly DEGs between enteroviral and bacterial meningitis. Random forest classification showed that we were able to differentiate enteroviral from bacterial meningitis with an AUC of 0.982 on held-out samples. CONCLUSIONS: Enteroviral meningitis has an innate immunity signature with type 1 interferons as key players. Our classifier, based on blood host transcriptomic profiles of different meningitis cases, is a possible strong alternative for diagnosing enteroviral meningitis.

speaq 2.0: A complete workflow for high-throughput 1D NMR spectra processing and quantification.

Nuclear Magnetic Resonance (NMR) spectroscopy is, together with liquid chromatography-mass spectrometry (LC-MS), the most established platform to perform metabolomics. In contrast to LC-MS however, NMR data is predominantly being processed with commercial software. Meanwhile its data processing remains tedious and dependent on user interventions. As a follow-up to speaq, a previously released workflow for NMR spectral alignment and quantitation, we present speaq 2.0. This completely revised framework to automatically analyze 1D NMR spectra uses wavelets to efficiently summarize the raw spectra with minimal information loss or user interaction. The tool offers a fast and easy workflow that starts with the common approach of peak-picking, followed by grouping, thus avoiding the binning step. This yields a matrix consisting of features, samples and peak values that can be conveniently processed either by using included multivariate statistical functions or by using many other recently developed methods for NMR data analysis. speaq 2.0 facilitates robust and high-throughput metabolomics based on 1D NMR but is also compatible with other NMR frameworks or complementary LC-MS workflows. The methods are benchmarked using a simulated dataset and two publicly available datasets. speaq 2.0 is distributed through the existing speaq R package to provide a complete solution for NMR data processing. The package and the code for the presented case studies are freely available on CRAN (https://cran.r-project.org/package=speaq) and GitHub (https://github.com/beirnaert/speaq).

Fast Open Modification Spectral Library Searching through Approximate Nearest Neighbor Indexing.

Open modification searching (OMS) is a powerful search strategy that identifies peptides carrying any type of modification by allowing a modified spectrum to match against its unmodified variant by using a very wide precursor mass window. A drawback of this strategy, however, is that it leads to a large increase in search time. Although performing an open search can be done using existing spectral library search engines by simply setting a wide precursor mass window, none of these tools have been optimized for OMS, leading to excessive runtimes and suboptimal identification results. We present the ANN-SoLo tool for fast and accurate open spectral library searching. ANN-SoLo uses approximate nearest neighbor indexing to speed up OMS by selecting only a limited number of the most relevant library spectra to compare to an unknown query spectrum. This approach is combined with a cascade search strategy to maximize the number of identified unmodified and modified spectra while strictly controlling the false discovery rate as well as a shifted dot product score to sensitively match modified spectra to their unmodified counterparts. ANN-SoLo achieves state-of-the-art performance in terms of speed and the number of identifications. On a previously published human cell line data set, ANN-SoLo confidently identifies more spectra than SpectraST or MSFragger and achieves a speedup of an order of magnitude compared with SpectraST. ANN-SoLo is implemented in Python and C++. It is freely available under the Apache 2.0 license at https://github.com/bittremieux/ANN-SoLo .

On the feasibility of mining CD8+ T cell receptor patterns underlying immunogenic peptide recognition.

Current T cell epitope prediction tools are a valuable resource in designing targeted immunogenicity experiments. They typically focus on, and are able to, accurately predict peptide binding and presentation by major histocompatibility complex (MHC) molecules on the surface of antigen-presenting cells. However, recognition of the peptide-MHC complex by a T cell receptor (TCR) is often not included in these tools. We developed a classification approach based on random forest classifiers to predict recognition of a peptide by a T cell receptor and discover patterns that contribute to recognition. We considered two approaches to solve this problem: (1) distinguishing between two sets of TCRs that each bind to a known peptide and (2) retrieving TCRs that bind to a given peptide from a large pool of TCRs. Evaluation of the models on two HIV-1, B*08-restricted epitopes reveals good performance and hints towards structural CDR3 features that can determine peptide immunogenicity. These results are of particular importance as they show that prediction of T cell epitope and T cell epitope recognition based on sequence data is a feasible approach. In addition, the validity of our models not only serves as a proof of concept for the prediction of immunogenic T cell epitopes but also paves the way for more general and high-performing models.

Increased herpes zoster risk associated with poor HLA-A immediate early 62 protein (IE62) affinity.

Around 30% of individuals will develop herpes zoster (HZ), caused by the varicella zoster virus (VZV), during their life. While several risk factors for HZ, such as immunosuppressive therapy, are well known, the genetic and molecular components that determine the risk of otherwise healthy individuals to develop HZ are still poorly understood. We created a computational model for the Human Leukocyte Antigen (HLA-A, -B, and -C) presentation capacity of peptides derived from the VZV Immediate Early 62 (IE62) protein. This model could then be applied to a HZ cohort with known HLA molecules. We found that HLA-A molecules with poor VZV IE62 presentation capabilities were more common in a cohort of 50 individuals with a history of HZ compared to a nationwide control group, which equated to a HZ risk increase of 60%. This tendency was most pronounced for cases of HZ at a young age, where other risk factors are less prevalent. These findings provide new molecular insights into the development of HZ and reveal a genetic predisposition in those individuals most at risk to develop HZ.

sfinx: an R package for the elimination of false positives from affinity purification-mass spectrometry datasets.

SUMMARY: We describe sfinx, an R package providing access to the straightforward filtering index (SFINX) for the separation of true positive from false positive protein interactions in affinity purification - mass spectrometry datasets. This package maintains the reliability and user-friendliness of the SFINX web site interface but is faster, unlimited in input size, and can be run locally within R. AVAILABILITY AND IMPLEMENTATION: The sfinx R package is available for download at the comprehensive R archive network (CRAN) https://cran.r-project.org/web/packages/sfinx/ under the Apache License 2.0. CONTACT: sven.eyckerman@vib-ugent.be or kevin.titeca@gmail.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Protein complex analysis: From raw protein lists to protein interaction networks.

The elucidation of molecular interaction networks is one of the pivotal challenges in the study of biology. Affinity purification-mass spectrometry and other co-complex methods have become widely employed experimental techniques to identify protein complexes. These techniques typically suffer from a high number of false negatives and false positive contaminants due to technical shortcomings and purification biases. To support a diverse range of experimental designs and approaches, a large number of computational methods have been proposed to filter, infer and validate protein interaction networks from experimental pull-down MS data. Nevertheless, this expansion of available methods complicates the selection of the most optimal ones to support systems biology-driven knowledge extraction. In this review, we give an overview of the most commonly used computational methods to process and interpret co-complex results, and we discuss the issues and unsolved problems that still exist within the field. © 2015 Wiley Periodicals, Inc. Mass Spec Rev 36:600-614, 2017.

COLOMBOS v3.0: leveraging gene expression compendia for cross-species analyses.

COLOMBOS is a database that integrates publicly available transcriptomics data for several prokaryotic model organisms. Compared to the previous version it has more than doubled in size, both in terms of species and data available. The manually curated condition annotation has been overhauled as well, giving more complete information about samples' experimental conditions and their differences. Functionality-wise cross-species analyses now enable users to analyse expression data for all species simultaneously, and identify candidate genes with evolutionary conserved expression behaviour. All the expression-based query tools have undergone a substantial improvement, overcoming the limit of enforced co-expression data retrieval and instead enabling the return of more complex patterns of expression behaviour. COLOMBOS is freely available through a web application at http://colombos.net/. The complete database is also accessible via REST API or downloadable as tab-delimited text files.

A primer to frequent itemset mining for bioinformatics.

Over the past two decades, pattern mining techniques have become an integral part of many bioinformatics solutions. Frequent itemset mining is a popular group of pattern mining techniques designed to identify elements that frequently co-occur. An archetypical example is the identification of products that often end up together in the same shopping basket in supermarket transactions. A number of algorithms have been developed to address variations of this computationally non-trivial problem. Frequent itemset mining techniques are able to efficiently capture the characteristics of (complex) data and succinctly summarize it. Owing to these and other interesting properties, these techniques have proven their value in biological data analysis. Nevertheless, information about the bioinformatics applications of these techniques remains scattered. In this primer, we introduce frequent itemset mining and their derived association rules for life scientists. We give an overview of various algorithms, and illustrate how they can be used in several real-life bioinformatics application domains. We end with a discussion of the future potential and open challenges for frequent itemset mining in the life sciences.