Single-cell functional genomics reveals determinants of sensitivity and resistance to natural killer cells in blood cancers.

Cancer cells can evade natural killer (NK) cell activity, thereby limiting anti-tumor immunity. To reveal genetic determinants of susceptibility to NK cell activity, we examined interacting NK cells and blood cancer cells using single-cell and genome-scale functional genomics screens. Interaction of NK and cancer cells induced distinct activation and type I interferon (IFN) states in both cell types depending on the cancer cell lineage and molecular phenotype, ranging from more sensitive myeloid to less sensitive B-lymphoid cancers. CRISPR screens in cancer cells uncovered genes regulating sensitivity and resistance to NK cell-mediated killing, including adhesion-related glycoproteins, protein fucosylation genes, and transcriptional regulators, in addition to confirming the importance of antigen presentation and death receptor signaling pathways. CRISPR screens with a single-cell transcriptomic readout provided insight into underlying mechanisms, including regulation of IFN-γ signaling in cancer cells and NK cell activation states. Our findings highlight the diversity of mechanisms influencing NK cell susceptibility across different cancers and provide a resource for NK cell-based therapies.

Variation in Microbiome LPS Immunogenicity Contributes to Autoimmunity in Humans.

According to the hygiene hypothesis, the increasing incidence of autoimmune diseases in western countries may be explained by changes in early microbial exposure, leading to altered immune maturation. We followed gut microbiome development from birth until age three in 222 infants in Northern Europe, where early-onset autoimmune diseases are common in Finland and Estonia but are less prevalent in Russia. We found that Bacteroides species are lowly abundant in Russians but dominate in Finnish and Estonian infants. Therefore, their lipopolysaccharide (LPS) exposures arose primarily from Bacteroides rather than from Escherichia coli, which is a potent innate immune activator. We show that Bacteroides LPS is structurally distinct from E. coli LPS and inhibits innate immune signaling and endotoxin tolerance; furthermore, unlike LPS from E. coli, B. dorei LPS does not decrease incidence of autoimmune diabetes in non-obese diabetic mice. Early colonization by immunologically silencing microbiota may thus preclude aspects of immune education.

TET proteins regulate the lineage specification and TCR-mediated expansion of iNKT cells.

TET proteins oxidize 5-methylcytosine in DNA to 5-hydroxymethylcytosine and other oxidation products. We found that simultaneous deletion of Tet2 and Tet3 in mouse CD4+CD8+ double-positive thymocytes resulted in dysregulated development and proliferation of invariant natural killer T cells (iNKT cells). Tet2-Tet3 double-knockout (DKO) iNKT cells displayed pronounced skewing toward the NKT17 lineage, with increased DNA methylation and impaired expression of genes encoding the key lineage-specifying factors T-bet and ThPOK. Transfer of purified Tet2-Tet3 DKO iNKT cells into immunocompetent recipient mice resulted in an uncontrolled expansion that was dependent on the nonclassical major histocompatibility complex (MHC) protein CD1d, which presents lipid antigens to iNKT cells. Our data indicate that TET proteins regulate iNKT cell fate by ensuring their proper development and maturation and by suppressing aberrant proliferation mediated by the T cell antigen receptor (TCR).

Quantitative proteomics analysis of signalosome dynamics in primary T cells identifies the surface receptor CD6 as a Lat adaptor-independent TCR signaling hub.

T cell antigen receptor (TCR)-mediated activation of T cells requires the interaction of dozens of proteins. Here we used quantitative mass spectrometry and activated primary CD4(+) T cells from mice in which a tag for affinity purification was knocked into several genes to determine the composition and dynamics of multiprotein complexes that formed around the kinase Zap70 and the adaptors Lat and SLP-76. Most of the 112 high-confidence time-resolved protein interactions we observed were previously unknown. The surface receptor CD6 was able to initiate its own signaling pathway by recruiting SLP-76 and the guanine nucleotide-exchange factor Vav1 regardless of the presence of Lat. Our findings provide a more complete model of TCR signaling in which CD6 constitutes a signaling hub that contributes to the diversification of TCR signaling.

EPIC-TRACE: predicting TCR binding to unseen epitopes using attention and contextualized embeddings.

MOTIVATION: T cells play an essential role in adaptive immune system to fight pathogens and cancer but may also give rise to autoimmune diseases. The recognition of a peptide-MHC (pMHC) complex by a T cell receptor (TCR) is required to elicit an immune response. Many machine learning models have been developed to predict the binding, but generalizing predictions to pMHCs outside the training data remains challenging. RESULTS: We have developed a new machine learning model that utilizes information about the TCR from both α and ß chains, epitope sequence, and MHC. Our method uses ProtBERT embeddings for the amino acid sequences of both chains and the epitope, as well as convolution and multi-head attention architectures. We show the importance of each input feature as well as the benefit of including epitopes with only a few TCRs to the training data. We evaluate our model on existing databases and show that it compares favorably against other state-of-the-art models. AVAILABILITY AND IMPLEMENTATION: https://github.com/DaniTheOrange/EPIC-TRACE.

TSignal: a transformer model for signal peptide prediction.

MOTIVATION: Signal peptides (SPs) are short amino acid segments present at the N-terminus of newly synthesized proteins that facilitate protein translocation into the lumen of the endoplasmic reticulum, after which they are cleaved off. Specific regions of SPs influence the efficiency of protein translocation, and small changes in their primary structure can abolish protein secretion altogether. The lack of conserved motifs across SPs, sensitivity to mutations, and variability in the length of the peptides make SP prediction a challenging task that has been extensively pursued over the years. RESULTS: We introduce TSignal, a deep transformer-based neural network architecture that utilizes BERT language models and dot-product attention techniques. TSignal predicts the presence of SPs and the cleavage site between the SP and the translocated mature protein. We use common benchmark datasets and show competitive accuracy in terms of SP presence prediction and state-of-the-art accuracy in terms of cleavage site prediction for most of the SP types and organism groups. We further illustrate that our fully data-driven trained model identifies useful biological information on heterogeneous test sequences. AVAILABILITY AND IMPLEMENTATION: TSignal is available at: https://github.com/Dumitrescu-Alexandru/TSignal.

TCRconv: predicting recognition between T cell receptors and epitopes using contextualized motifs.

MOTIVATION: T cells use T cell receptors (TCRs) to recognize small parts of antigens, called epitopes, presented by major histocompatibility complexes. Once an epitope is recognized, an immune response is initiated and T cell activation and proliferation by clonal expansion begin. Clonal populations of T cells with identical TCRs can remain in the body for years, thus forming immunological memory and potentially mappable immunological signatures, which could have implications in clinical applications including infectious diseases, autoimmunity and tumor immunology. RESULTS: We introduce TCRconv, a deep learning model for predicting recognition between TCRs and epitopes. TCRconv uses a deep protein language model and convolutions to extract contextualized motifs and provides state-of-the-art TCR-epitope prediction accuracy. Using TCR repertoires from COVID-19 patients, we demonstrate that TCRconv can provide insight into T cell dynamics and phenotypes during the disease. AVAILABILITY AND IMPLEMENTATION: TCRconv is available at https://github.com/emmijokinen/tcrconv. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The transcription factor NFAT promotes exhaustion of activated CD8⁺ T cells.

During persistent antigen stimulation, CD8(+) T cells show a gradual decrease in effector function, referred to as exhaustion, which impairs responses in the setting of tumors and infections. Here we demonstrate that the transcription factor NFAT controls the program of T cell exhaustion. When expressed in cells, an engineered form of NFAT1 unable to interact with AP-1 transcription factors diminished T cell receptor (TCR) signaling, increased the expression of inhibitory cell surface receptors, and interfered with the ability of CD8(+) T cells to protect against Listeria infection and attenuate tumor growth in vivo. We defined the genomic regions occupied by endogenous and engineered NFAT1 in primary CD8(+) T cells and showed that genes directly induced by the engineered NFAT1 overlapped with genes expressed in exhausted CD8(+) T cells in vivo. Our data show that NFAT promotes T cell anergy and exhaustion by binding at sites that do not require cooperation with AP-1.

The human gut microbiome in early-onset type 1 diabetes from the TEDDY study.

Type 1 diabetes (T1D) is an autoimmune disease that targets pancreatic islet beta cells and incorporates genetic and environmental factors1, including complex genetic elements2, patient exposures3 and the gut microbiome4. Viral infections5 and broader gut dysbioses6 have been identified as potential causes or contributing factors; however, human studies have not yet identified microbial compositional or functional triggers that are predictive of islet autoimmunity or T1D. Here we analyse 10,913 metagenomes in stool samples from 783 mostly white, non-Hispanic children. The samples were collected monthly from three months of age until the clinical end point (islet autoimmunity or T1D) in the The Environmental Determinants of Diabetes in the Young (TEDDY) study, to characterize the natural history of the early gut microbiome in connection to islet autoimmunity, T1D diagnosis, and other common early life events such as antibiotic treatments and probiotics. The microbiomes of control children contained more genes that were related to fermentation and the biosynthesis of short-chain fatty acids, but these were not consistently associated with particular taxa across geographically diverse clinical centres, suggesting that microbial factors associated with T1D are taxonomically diffuse but functionally more coherent. When we investigated the broader establishment and development of the infant microbiome, both taxonomic and functional profiles were dynamic and highly individualized, and dominated in the first year of life by one of three largely exclusive Bifidobacterium species (B. bifidum, B. breve or B. longum) or by the phylum Proteobacteria. In particular, the strain-specific carriage of genes for the utilization of human milk oligosaccharide within a subset of B. longum was present specifically in breast-fed infants. These analyses of TEDDY gut metagenomes provide, to our knowledge, the largest and most detailed longitudinal functional profile of the developing gut microbiome in relation to islet autoimmunity, T1D and other early childhood events. Together with existing evidence from human cohorts7,8 and a T1D mouse model9, these data support the protective effects of short-chain fatty acids in early-onset human T1D.

LuxHMM: DNA methylation analysis with genome segmentation via hidden Markov model.

BACKGROUND: DNA methylation plays an important role in studying the epigenetics of various biological processes including many diseases. Although differential methylation of individual cytosines can be informative, given that methylation of neighboring CpGs are typically correlated, analysis of differentially methylated regions is often of more interest. RESULTS: We have developed a probabilistic method and software, LuxHMM, that uses hidden Markov model (HMM) to segment the genome into regions and a Bayesian regression model, which allows handling of multiple covariates, to infer differential methylation of regions. Moreover, our model includes experimental parameters that describe the underlying biochemistry in bisulfite sequencing and model inference is done using either variational inference for efficient genome-scale analysis or Hamiltonian Monte Carlo (HMC). CONCLUSIONS: Analyses of real and simulated bisulfite sequencing data demonstrate the competitive performance of LuxHMM compared with other published differential methylation analysis methods.

ChromDMM: a Dirichlet-multinomial mixture model for clustering heterogeneous epigenetic data.

MOTIVATION: Research on epigenetic modifications and other chromatin features at genomic regulatory elements elucidates essential biological mechanisms including the regulation of gene expression. Despite the growing number of epigenetic datasets, new tools are still needed to discover novel distinctive patterns of heterogeneous epigenetic signals at regulatory elements. RESULTS: We introduce ChromDMM, a product Dirichlet-multinomial mixture model for clustering genomic regions that are characterized by multiple chromatin features. ChromDMM extends the mixture model framework by profile shifting and flipping that can probabilistically account for inaccuracies in the position and strand-orientation of the genomic regions. Owing to hyper-parameter optimization, ChromDMM can also regularize the smoothness of the epigenetic profiles across the consecutive genomic regions. With simulated data, we demonstrate that ChromDMM clusters, shifts and strand-orients the profiles more accurately than previous methods. With ENCODE data, we show that the clustering of enhancer regions in the human genome reveals distinct patterns in several chromatin features. We further validate the enhancer clusters by their enrichment for transcriptional regulatory factor binding sites. AVAILABILITY AND IMPLEMENTATION: ChromDMM is implemented as an R package and is available at https://github.com/MariaOsmala/ChromDMM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Probabilistic modeling methods for cell-free DNA methylation based cancer classification.

BACKGROUND: cfMeDIP-seq is a low-cost method for determining the DNA methylation status of cell-free DNA and it has been successfully combined with statistical methods for accurate cancer diagnostics. We investigate the diagnostic classification aspect by applying statistical tests and dimension reduction techniques for feature selection and probabilistic modeling for the cancer type classification, and we also study the effect of sequencing depth. METHODS: We experiment with a variety of statistical methods that use different feature selection and feature extraction methods as well as probabilistic classifiers for diagnostic decision making. We test the (moderated) t-tests and the Fisher's exact test for feature selection, principal component analysis (PCA) as well as iterative supervised PCA (ISPCA) for feature generation, and GLMnet and logistic regression methods with sparsity promoting priors for classification. Probabilistic programming language Stan is used to implement Bayesian inference for the probabilistic models. RESULTS AND CONCLUSIONS: We compare overlaps of differentially methylated genomic regions as chosen by different feature selection methods, and evaluate probabilistic classifiers by evaluating the area under the receiver operating characteristic scores on discovery and validation cohorts. While we observe that many methods perform equally well as, and occasionally considerably better than, GLMnet that was originally proposed for cfMeDIP-seq based cancer classification, we also observed that performance of different methods vary across sequencing depths, cancer types and study cohorts. Overall, methods that seem robust and promising include Fisher's exact test and ISPCA for feature selection as well as a simple logistic regression model with the number of hyper and hypo-methylated regions as features.

Modeling binding specificities of transcription factor pairs with random forests.

BACKGROUND: Transcription factors (TFs) bind regulatory DNA regions with sequence specificity, form complexes and regulate gene expression. In cooperative TF-TF binding, two transcription factors bind onto a shared DNA binding site as a pair. Previous work has demonstrated pairwise TF-TF-DNA interactions with position weight matrices (PWMs), which may however not sufficiently take into account the complexity and flexibility of pairwise binding. RESULTS: We propose two random forest (RF) methods for joint TF-TF binding site prediction: ComBind and JointRF. We train models with previously published large-scale CAP-SELEX DNA libraries, which comprise DNA sequences enriched for binding of a selected TF pair. JointRF builds a random forest with sub-sequences selected from CAP-SELEX DNA reads with previously proposed pairwise PWM. JointRF outperforms (area under receiver operating characteristics curve, AUROC, 0.75) the current state-of-the-art method i.e. orientation and spacing specific pairwise PWMs (AUROC 0.59). Thus, JointRF may be utilized to improve prediction accuracy for pre-determined binding preferences. However, pairwise TF binding is currently considered flexible; a pair may bind DNA with different orientations and amounts of dinucleotide gaps or overlap between the two motifs. Thus, we developed ComBind, which utilizes random forests by considering simultaneously multiple orientations and spacings of the two factors. Our approach outperforms (AUROC 0.78) PWMs, as well as JointRF (p<0.00195). ComBind provides an approach for predicting TF-TF binding sites without prior knowledge on pairwise binding preferences. However, more research is needed to assess ComBind eligibility for practical applications. CONCLUSIONS: Random forest is well suited for modeling pairwise TF-TF-DNA binding specificities, and ComBind provides an improvement to pairwise binding site prediction accuracy.

LuxRep: a technical replicate-aware method for bisulfite sequencing data analysis.

BACKGROUND: DNA methylation is commonly measured using bisulfite sequencing (BS-seq). The quality of a BS-seq library is measured by its bisulfite conversion efficiency. Libraries with low conversion rates are typically excluded from analysis resulting in reduced coverage and increased costs. RESULTS: We have developed a probabilistic method and software, LuxRep, that implements a general linear model and simultaneously accounts for technical replicates (libraries from the same biological sample) from different bisulfite-converted DNA libraries. Using simulations and actual DNA methylation data, we show that including technical replicates with low bisulfite conversion rates generates more accurate estimates of methylation levels and differentially methylated sites. Moreover, using variational inference speeds up computation time necessary for whole genome analysis. CONCLUSIONS: In this work we show that taking into account technical replicates (i.e. libraries) of BS-seq data of varying bisulfite conversion rates, with their corresponding experimental parameters, improves methylation level estimation and differential methylation detection.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age.

AIMS/HYPOTHESIS: Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies. METHODS: Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4+ T cell, CD8+ T cell and CD4-CD8- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing. RESULTS: We identified 79, 56 and 45 differentially methylated regions in CD4+ T cells, CD8+ T cells and CD4-CD8- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4+ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4+ T cells. CONCLUSIONS/INTERPRETATION: These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

Umbilical cord blood DNA methylation in children who later develop type 1 diabetes.

AIMS/HYPOTHESIS: Distinct DNA methylation patterns have recently been observed to precede type 1 diabetes in whole blood collected from young children. Our aim was to determine whether perinatal DNA methylation is associated with later progression to type 1 diabetes. METHODS: Reduced representation bisulphite sequencing (RRBS) analysis was performed on umbilical cord blood samples collected within the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) Study. Children later diagnosed with type 1 diabetes and/or who tested positive for multiple islet autoantibodies (n = 43) were compared with control individuals (n = 79) who remained autoantibody-negative throughout the DIPP follow-up until 15 years of age. Potential confounding factors related to the pregnancy and the mother were included in the analysis. RESULTS: No differences in the umbilical cord blood methylation patterns were observed between the cases and controls at a false discovery rate <0.05. CONCLUSIONS/INTERPRETATION: Based on our results, differences between children who progress to type 1 diabetes and those who remain healthy throughout childhood are not yet present in the perinatal DNA methylome. However, we cannot exclude the possibility that such differences would be found in a larger dataset.

lgpr: an interpretable non-parametric method for inferring covariate effects from longitudinal data.

MOTIVATION: Longitudinal study designs are indispensable for studying disease progression. Inferring covariate effects from longitudinal data, however, requires interpretable methods that can model complicated covariance structures and detect non-linear effects of both categorical and continuous covariates, as well as their interactions. Detecting disease effects is hindered by the fact that they often occur rapidly near the disease initiation time, and this time point cannot be exactly observed. An additional challenge is that the effect magnitude can be heterogeneous over the subjects. RESULTS: We present lgpr, a widely applicable and interpretable method for non-parametric analysis of longitudinal data using additive Gaussian processes. We demonstrate that it outperforms previous approaches in identifying the relevant categorical and continuous covariates in various settings. Furthermore, it implements important novel features, including the ability to account for the heterogeneity of covariate effects, their temporal uncertainty, and appropriate observation models for different types of biomedical data. The lgpr tool is implemented as a comprehensive and user-friendly R-package. AVAILABILITY AND IMPLEMENTATION: lgpr is available at jtimonen.github.io/lgpr-usage with documentation, tutorials, test data and code for reproducing the experiments of this article. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Global chromatin state analysis reveals lineage-specific enhancers during the initiation of human T helper 1 and T helper 2 cell polarization.

Naive CD4⁺ T cells can differentiate into specific helper and regulatory T cell lineages in order to combat infection and disease. The correct response to cytokines and a controlled balance of these populations is critical for the immune system and the avoidance of autoimmune disorders. To investigate how early cell-fate commitment is regulated, we generated the first human genome-wide maps of histone modifications that reveal enhancer elements after 72 hr of in vitro polarization toward T helper 1 (Th1) and T helper 2 (Th2) cell lineages. Our analysis indicated that even at this very early time point, cell-specific gene regulation and enhancers were at work directing lineage commitment. Further examination of lineage-specific enhancers identified transcription factors (TFs) with known and unknown T cell roles as putative drivers of lineage-specific gene expression. Lastly, an integrative analysis of immunopathogenic-associated SNPs suggests a role for distal regulatory elements in disease etiology.

Predicting recognition between T cell receptors and epitopes with TCRGP.

Adaptive immune system uses T cell receptors (TCRs) to recognize pathogens and to consequently initiate immune responses. TCRs can be sequenced from individuals and methods analyzing the specificity of the TCRs can help us better understand individuals' immune status in different disorders. For this task, we have developed TCRGP, a novel Gaussian process method that predicts if TCRs recognize specified epitopes. TCRGP can utilize the amino acid sequences of the complementarity determining regions (CDRs) from TCRα and TCRß chains and learn which CDRs are important in recognizing different epitopes. Our comprehensive evaluation with epitope-specific TCR sequencing data shows that TCRGP achieves on average higher prediction accuracy in terms of AUROC score than existing state-of-the-art methods in epitope-specificity predictions. We also propose a novel analysis approach for combined single-cell RNA and TCRαß (scRNA+TCRαß) sequencing data by quantifying epitope-specific TCRs with TCRGP and identify HBV-epitope specific T cells and their transcriptomic states in hepatocellular carcinoma patients.