Search

Suppressive function of regulatory T cells (Treg) is dependent on signaling of their antigen receptors triggered by cognate self, dietary, or microbial peptides presented on MHC II. However, it remains largely unknown whether distinct or shared repertoires of Treg TCRs are mobilized in response to different challenges in the same tissue or the same challenge in different tissues. Here we use a fixed TCRβ chain FoxP3-GFP mouse model to analyze conventional (eCD4) and regulatory (eTreg) effector TCRα repertoires in response to six distinct antigenic challenges to the lung and skin. This model shows highly ‘digital’ repertoire behavior with easy-to-track challenge-specific TCRα CDR3 clusters. For both eCD4 and eTreg subsets, we observe challenge-specific clonal expansions yielding homologous TCRα clusters within and across animals and exposure sites, which are also reflected in the draining lymph nodes but not systemically. Some CDR3 clusters are shared across cancer challenges, suggesting a response to common tumor-associated antigens. For most challenges, eCD4 and eTreg clonal response does not overlap. Such overlap is exclusively observed at the sites of certain tumor challenges, and not systematically, suggesting transient and local tumor-induced eCD4=>eTreg plasticity. This transition includes a dominant tumor-responding eCD4 CDR3 motif, as well as characteristic iNKT TCRα CDR3. In addition, we examine the homeostatic tissue residency of clonal eTreg populations by excluding the site of challenge from our analysis. We demonstrate that distinct CDR3 motifs are characteristic of eTreg cells residing in particular lymphatic tissues, regardless of the challenge. This observation reveals the tissue-resident, antigen-specific clonal Treg populations.

An individual’s T-cell repertoire constantly changes under the influence of external and internal factors. Cells that do not receive a stimulatory signal die, while those that encounter and recognize a pathogen or receive a co-stimulatory signal divide, resulting in clonal expansions. T-cell clones can be traced by monitoring the presence of their unique T-cell receptor (TCR) sequence, which is assembled de novo through a process known as V(D)J rearrangement. Tracking T cells can provide valuable insights into the survival of cells after hematopoietic stem cell transplantation (HSCT) or cancer treatment response and can indicate the induction of protective immunity by vaccination. In this study, we report a bioinformatic method for quantifying the T-cell repertoire dynamics from TCR sequencing data. We demonstrate its utility by measuring the T-cell repertoire stability in healthy donors, by quantifying the effect of donor lymphocyte infusion (DLI), and by tracking the fate of the different T-cell subsets in HSCT patients and the expansion of pathogen-specific clones in vaccinated individuals.

IntroductionT-cell receptor (TCR) recognition of foreign peptides presented by the major histocompatibility complex (MHC) initiates the adaptive immune response against pathogens. While a large number of TCR sequences specific to different antigenic peptides are known to date, the structural data describing the conformation and contacting residues for TCR-peptide-MHC complexes is relatively limited. In the present study we aim to extend and analyze the set of available structures by performing highly accurate template-based modeling of these complexes using TCR sequences with known specificity. MethodsIdentification of CDR3 sequences and their further clustering, based on available spatial structures, V- and J-genes of corresponding T-cell receptors, and epitopes, was performed using the VDJdb database. Modeling of the selected CDR3 loops was conducted using a stepwise introduction of single amino acid substitutions to the template PDB structures, followed by optimization of the TCR-peptide-MHC contacting interface using the Rosetta package applications. Statistical analysis and recursive feature elimination procedures were carried out on computed energy values and properties of contacting amino acid residues between CDR3 loops and peptides, using R.ResultsUsing the set of 29 complex templates (including a template with SARS-CoV-2 antigen) and 732 specificity records, we built a database of 1585 model structures carrying substitutions in either TCRα or TCRβ chains with some models representing the result of different mutation pathways for the same final structure. This database allowed us to analyze features of amino acid contacts in TCR - peptide interfaces that govern antigen recognition preferences and interpret these interactions in terms of physicochemical properties of interacting residues.ConclusionOur results provide a methodology for creating high-quality TCR-peptide-MHC models for antigens of interest that can be utilized to predict TCR specificity.

Abstract Background Coinfection with HIV and Plasmodium parasites is fairly common, but the sequence of infection with these two pathogens and their impact on disease progression are poorly understood. Methods A Chinese rhesus macaque HIV and Plasmodium coinfection model was established to compare the impact of pre-existing and subsequent malaria on the progression of SIV infection. Results We found that a pre-existing malaria caused animals to produce a greater number of CD4+CCR5+ T cells for SIV replication, resulting in higher viral loads. Conversely, subsequent malaria induced a substantially larger proportion of CD4+CD28highCD95high central memory T cells and a stronger SIV-specific T cell response, maintained the repertoire diversity of SIV-specific T cell receptors, and generated new SIV-specific T cell clonotypes to trace SIV antigenic variation, resulting in improved survival of SIV-infected animals. Conclusion The complex outcomes of this study may have important implications for research on human HIV and malaria coinfection. The infection order of the two pathogens (HIV and malaria parasites) should be emphasized. Video abstract

Summary: Common variable immune deficiency (CVID) is a heterogeneous disorder characterized by recurrent infections, low levels of serum immunoglobulins, and impaired vaccine responses. Autoimmune manifestations are common, but B cell central and peripheral selection mechanisms in CVID are incompletely understood. Here, we find that receptor editing, a measure of central tolerance, is increased in transitional B cells from CVID patients and that these cells have a higher immunoglobulin κ:λ ratio in CVID patients with autoimmune manifestations than in those with infection only. Contrariwise, the selection pressure in the germinal center on CD27bright memory B cells is decreased in CVID patients with autoimmune manifestations. Finally, functionally, T cell-dependent activation showed that naive B cells in CVID patients are badly equipped for activation and induction of mismatch repair genes. We conclude that central tolerance is functional whereas peripheral selection is defective in CVID patients with autoimmune manifestations, which could underpin the development of autoimmunity.

T-cell receptors (TCRs) are formed by stochastic gene rearrangements, theoretically generating >1019 sequences. They are selected during thymopoiesis, which releases a repertoire of about 108 unique TCRs per individual. How evolution shaped a process that produces TCRs that can effectively handle a countless and evolving set of infectious agents is a central question of immunology. The paradigm is that a diverse enough repertoire of TCRs should always provide a proper, though rare, specificity for any given need. Expansion of such rare T cells would provide enough fighters for an effective immune response and enough antigen-experienced cells for memory. We show here that human thymopoiesis releases a large population of clustered CD8+ T cells harboring α/β paired TCRs that (i) have high generation probabilities and (ii) a preferential usage of some V and J genes, (iii) which CDR3 are shared between individuals, and (iv) can each bind and be activated by multiple unrelated viral peptides, notably from EBV, CMV, and influenza. These polyspecific T cells may represent a first line of defense that is mobilized in response to infections before a more specific response subsequently ensures viral elimination. Our results support an evolutionary selection of polyspecific α/β TCRs for broad antiviral responses and heterologous immunity.

Human leukocyte antigen (HLA) genes are the most polymorphic loci in the human genome and code for proteins that play a key role in guiding adaptive immune responses by presenting foreign and self peptides (ligands) to T cells. Each person carries up to 6 HLA class I variants (maternal and paternal copies of HLA-A, HLA-B and HLA-C genes) and also multiple HLA class II variants, which cumulatively define the landscape of peptides presented to T cells. Each HLA variant has its own repertoire of presented peptides with a certain sequence motif which is mainly defined by peptide anchor residues (typically the second and the last positions for HLA class I ligands) forming key interactions with the peptide-binding groove of HLA. In this study, we aimed to characterize HLA binding preferences in terms of molecular functions of presented proteins. To focus on the ligand presentation bias introduced specifically by HLA-peptide interaction we performed large-scale in silico predictions of binding of all peptides from human proteome for a wide range of HLA variants and established which functions are characteristic for proteins that are more or less preferentially presented by different HLA variants using statistical calculations and gene ontology (GO) analysis. We demonstrated marked distinctions between HLA variants in molecular functions of preferentially presented proteins (e.g. some HLA variants preferentially present membrane and receptor proteins, while others – ribosomal and DNA-binding proteins) and reduced presentation of extracellular matrix and collagen proteins by the majority of HLA variants. To explain these observations we demonstrated that HLA preferentially presents proteins enriched in amino acids which are required as anchor residues for the particular HLA variant. Our observations can be extrapolated to explain the protective effect of certain HLA alleles in infectious diseases, and we hypothesize that they can also explain susceptibility to certain autoimmune diseases and cancers. We demonstrate that these differences lead to differential presentation of HIV, influenza virus, SARS-CoV-1 and SARS-CoV-2 proteins by various HLA alleles. Taking into consideration that HLA alleles are inherited in haplotypes, we hypothesized that haplotypes composed of a combination of HLA variants with different presentation preferences should be more advantageous as they allow presenting a larger repertoire of peptides and avoiding holes in immunopeptidome. Indeed, we demonstrated that HLA-A/HLA-B and HLA-A/HLA-C haplotypes which have a high frequency in the human population are comprised of HLA variants that are more distinct in terms of functions of preferentially presented proteins than the control pairs.

Spondyloarthritis (SpA) comprises a number of inflammatory rheumatic diseases with overlapping clinical manifestations. Strong association with several HLA-I alleles and T cell infiltration into an inflamed joint suggest involvement of T cells in SpA pathogenesis. In this study, we performed high-throughput T cell repertoire profiling of synovial fluid (SF) and peripheral blood (PB) samples collected from a large cohort of SpA patients. We showed that synovial fluid is enriched with expanded T cell clones that are shared between patients with similar HLA genotypes and persist during recurrent synovitis. Using an algorithm for identification of TCRs involved in immune response we discovered several antigen-driven CD8+ clonal groups associated with risk HLA-B*27 or HLA-B*38 alleles. We further show that these clonal groups were enriched in SF and had higher frequency in PB of SpA patients vs healthy donors, implying their relevance to SpA pathogenesis. Several of the groups were shared among patients with different SpAs that suggests a common immunopathological mechanism of the diseases. In summary, our results provide evidence for the role of specific CD8+ T cell clones in pathogenesis of SpA.

The stability and plasticity of B cell-mediated immune memory ensures the ability to respond to the repeated challenges. We have analyzed the longitudinal dynamics of immunoglobulin heavy chain repertoires from memory B cells, plasmablasts, and plasma cells from the peripheral blood of generally healthy volunteers. We reveal a high degree of clonal persistence in individual memory B cell subsets, with inter-individual convergence in memory and antibody-secreting cells (ASCs). ASC clonotypes demonstrate clonal relatedness to memory B cells, and are transient in peripheral blood. We identify two clusters of expanded clonal lineages with differing prevalence of memory B cells, isotypes, and persistence. Phylogenetic analysis revealed signs of reactivation of persisting memory B cell-enriched clonal lineages, accompanied by new rounds of affinity maturation during proliferation and differentiation into ASCs. Negative selection contributes to both persisting and reactivated lineages, preserving the functionality and specificity of B cell receptors (BCRs) to protect against current and future pathogens.

Insulinomatosis is characterized by monohormonality of multiple macro-tumors and micro-tumors that arise synchronously and metachronously in all regions of the pancreas, and often recurring hypoglycemia. One of the main characteristics of insulinomatosis is the presence of insulin-expressing monohormonal endocrine cell clusters that are exclusively composed of proliferating insulin-positive cells, are less than 1 mm in size, and show solid islet-like structure. It is presumed that insulinomatosis affects the entire population of β-cells. With regards to molecular genetics, this phenomenon is not related to mutation in MEN1 gene and is more similar to sporadic benign insulinomas, however, at the moment molecular genetics of this disease remains poorly investigated. NGS sequencing was performed with a panel of 409 cancer-related genes. Results of sequencing were analyzed by bioinformatic algorithms for detecting point mutations and copy number variations. DNA copy number variations were detected that harbor a large number of genes in insulinoma and fewer genes in micro-tumors. qPCR was used to confirm copy number variations at ATRX, FOXL2, IRS2 and CEBPA genes. Copy number alterations involving FOXL2, IRS2, CEBPA and ATRX genes were observed in insulinoma as well as in micro-tumors samples, suggesting that alterations of these genes may promote malignization in the β-cells population.

Summary: Biology of the response to anti-CTLA-4 involves the dynamics of specific T cell clones. Reasons for clinical success and failure of this treatment are still largely unknown. Here, we quantified the dynamics of the T cell receptor (TCR) repertoire, throughout 4 weeks involving treatment with anti-CTLA-4, in a syngeneic mouse model for colorectal cancer. These dynamics show an initial increase in clonality in tandem with a decrease in diversity, effects which gradually subside. Furthermore, response to treatment is tightly connected to the shared and public parts of the T cell repertoire. We were able to recognize time-dependent behaviors of specific TCR sequences and cell types and to show the response is dominated by specific motifs. We see that a single, specific time point might be useful to inform a physician of the true response to treatmentThe research further highlights the importance of temporal analyses of the immune response.

Summary: T cell recognition of peptides presented by human leukocyte antigens (HLAs) is mediated by the highly variable T cell receptor (TCR). Despite this built-in TCR variability, individuals can mount immune responses against viral epitopes by using identical or highly related TCRs expressed on CD8+ T cells. Characterization of these TCRs has extended our understanding of the molecular mechanisms that govern the recognition of peptide-HLA. However, few examples exist for CD4+ T cells. Here, we investigate CD4+ T cell responses to the internal proteins of the influenza A virus that correlate with protective immunity. We identify five internal epitopes that are commonly recognized by CD4+ T cells in five HLA-DR1+ subjects and show conservation across viral strains and zoonotic reservoirs. TCR repertoire analysis demonstrates several shared gene usage biases underpinned by complementary biochemical features evident in a structural comparison. These epitopes are attractive targets for vaccination and other T cell therapies.

Antigen recognition by T-cells is guided by the T-cell receptor (TCR) heterodimer formed by α and β chains. A huge diversity of TCR sequences should be maintained by the immune system in order to be able to mount an effective response towards foreign pathogens, so, due to cooperative binding of α and β chains to the pathogen, any constraints on chain pairing can have a profound effect on immune repertoire structure, diversity and antigen specificity. By integrating available structural data and paired chain sequencing results we were able to show that there are almost no constraints on pairing in TCRαβ complexes, allowing naive T-cell repertoire to reach the highest possible diversity. Additional analysis reveals that the specific choice of contacting amino acids can still have a profound effect on complex conformation. Moreover, antigen-driven selection can distort the uniform landscape of chain pairing, while small, yet significant, differences in the pairing can be attributed to various specialized T-cell subsets such as MAIT and iNKT T-cells, as well as other TCR sets specific to certain antigens.

Summary: Memory B cells (MBCs) epitomize the adaptation of the immune system to the environment. We identify two MBC subsets in peripheral blood, CD27dull and CD27bright MBCs, whose frequency changes with age. Heavy chain variable region (VH) usage, somatic mutation frequency replacement-to-silent ratio, and CDR3 property changes, reflecting consecutive selection of highly antigen-specific, low cross-reactive antibody variants, all demonstrate that CD27dull and CD27bright MBCs represent sequential MBC developmental stages, and stringent antigen-driven pressure selects CD27dull into the CD27bright MBC pool. Dynamics of human MBCs are exploited in pregnancy, when 50% of maternal MBCs are lost and CD27dull MBCs transit to the more differentiated CD27bright stage. In the postpartum period, the maternal MBC pool is replenished by the expansion of persistent CD27dull clones. Thus, the stability and flexibility of human B cell memory is ensured by CD27dull MBCs that expand and differentiate in response to change. : Grimsholm et al. show that CD27dull and CD27bright represent sequential MBC developmental stages. T cell- and germinal center (GC)-independent CD27dull MBCs are the plastic source of strongly selected and GC-dependent CD27bright MBCs. CD27dull MBCs, able to expand and differentiate in response to change, ensure stability and flexibility of human B cell memory. Keywords: memory B cells, pregnancy, immunological memory, CD27, VH repertoire, immunodeficiency, aging, spleen, vaccine, germinal center

Abstract Background Adaptive immune responses to newly encountered pathogens depend on the mobilization of antigen-specific clonotypes from a vastly diverse pool of naive T cells. Using recent advances in immune repertoire sequencing technologies, models of the immune receptor rearrangement process, and a database of annotated T cell receptor (TCR) sequences with known specificities, we explored the baseline frequencies of T cells specific for defined human leukocyte antigen (HLA) class I-restricted epitopes in healthy individuals. Methods We used a database of TCR sequences with known antigen specificities and a probabilistic TCR rearrangement model to estimate the baseline frequencies of TCRs specific to distinct antigens epitopespecificT-cells. We verified our estimates using a publicly available collection of TCR repertoires from healthy individuals. We also interrogated a database of immunogenic and non-immunogenic peptides is used to link baseline T-cell frequencies with epitope immunogenicity. Results Our findings revealed a high degree of variability in the prevalence of T cells specific for different antigens that could be explained by the physicochemical properties of the corresponding HLA class I-bound peptides. The occurrence of certain rearrangements was influenced by ancestry and HLA class I restriction, and umbilical cord blood samples contained higher frequencies of common pathogen-specific TCRs. We also identified a quantitative link between specific T cell frequencies and the immunogenicity of cognate epitopes presented by defined HLA class I molecules. Conclusions Our results suggest that the population frequencies of specific T cells are strikingly non-uniform across epitopes that are known to elicit immune responses. This inference leads to a new definition of epitope immunogenicity based on specific TCR frequencies, which can be estimated with a high degree of accuracy in silico, thereby providing a novel framework to integrate computational and experimental genomics with basic and translational research efforts in the field of T cell immunology.

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires.

Abstract The accuracy with which DNA polymerase can replicate a template DNA sequence is an extremely important property that can vary by an order of magnitude from one enzyme to another. The rate of nucleotide misincorporation is shaped by multiple factors, including PCR conditions and proofreading capabilities, and proper assessment of polymerase error rate is essential for a wide range of sensitive PCR-based assays. In this paper, we describe a method for studying polymerase errors with exceptional resolution, which combines unique molecular identifier tagging and high-throughput sequencing. Our protocol is less laborious than commonly-used methods, and is also scalable, robust and accurate. In a series of nine PCR assays, we have measured a range of polymerase accuracies that is in line with previous observations. However, we were also able to comprehensively describe individual errors introduced by each polymerase after either 20 PCR cycles or a linear amplification, revealing specific substitution preferences and the diversity of PCR error frequency profiles. We also demonstrate that the detected high-frequency PCR errors are highly recurrent and that the position in the template sequence and polymerase-specific substitution preferences are among the major factors influencing the observed PCR error rate.

Abstract Background Recently we proposed efficient method to exclude undesirable primers at any stage of amplification reaction, here termed NOPE (NOnsense-mediated Primer Exclusion). According to this method, added oligonucleotide overlapping with the 3′-end of unwanted amplification primer (NOPE oligo) simultaneously provides a template for its elongation. This elongation disrupts specificity of unwanted primer, preventing its further participation in PCR. The suggested approach allows to rationally manage the course of PCR reactions in order to facilitate analysis of complex DNA mixtures as well as to perform multistage PCR bypassing intermediate purification steps. Results Here we apply NOPE method to DNA library preparation for the high-throughput sequencing (HTS) with the PCR-based introduction of unique molecular identifiers (UMI). We show that NOPE oligo efficiently neutralizes UMI-containing oligonucleotides after introduction of UMI into sample DNA molecules, thus allowing to proceed with further amplification steps without purification and associated loss of starting material. At the same time, NOPE oligo does not affect the efficiency of target PCR amplification. Conclusion We describe a simple, robust and cheap modification of UMI-labeled HTS libraries preparation procedure, that allows to bypass purification step and thus to preserve starting material which may be limited, e.g. circulating tumor DNA, circulating fetal DNA, or small amounts of isolated cells of interest. Furthermore, demonstrated simplicity and robustness of NOPE method should make it popular in various PCR protocols.

Recently developed molecular methods allow large-scale profiling of T-cell receptor (TCR) sequences that encode for antigen specificity and immunological memory of these cells. However, it is well-known that the even unperturbed TCR repertoire structure is extremely complex due to the high diversity of TCR rearrangements and multiple biases imprinted by VDJ rearrangement process. The latter gives rise to the phenomenon of “public” TCR clonotypes that can be shared across multiple individuals and non-trivial structure of the TCR similarity network. Here, we outline a framework for TCR sequencing data analysis that can control for these biases in order to infer TCRs that are involved in response to antigens of interest. We apply two previously published methods, ALICE and TCRNET, to detect groups of homologous TCRs that are enriched in samples of interest. Using an example dataset of donors with known HLA haplotype and CMV status, we demonstrate that by applying HLA restriction rules and matching against a database of TCRs with known antigen specificity, it is possible to robustly detect motifs of epitope-specific responses in individual repertoires. We also highlight potential shortcomings of TCR clustering methods and demonstrate that highly expanded TCRs should be individually assessed to get the full picture of antigen-specific response.

Hypervariable T cell receptors (TCRs) play a key role in adaptive immunity, recognizing a vast diversity of pathogen-derived antigens. Our ability to extract clinically relevant information from large high-throughput sequencing of TCR repertoires (RepSeq) data is limited, because little is known about TCR-disease associations. We present Antigen-specific Lymphocyte Identification by Clustering of Expanded sequences (ALICE), a statistical approach that identifies TCR sequences actively involved in current immune responses from a single RepSeq sample and apply it to repertoires of patients with a variety of disorders - patients with autoimmune disease (ankylosing spondylitis [AS]), under cancer immunotherapy, or subject to an acute infection (live yellow fever [YF] vaccine). We validate the method with independent assays. ALICE requires no longitudinal data collection nor large cohorts, and it is directly applicable to most RepSeq datasets. Its results facilitate the identification of TCR variants associated with diseases and conditions, which can be used for diagnostics and rational vaccine design.

Age-related changes can significantly alter the state of adaptive immune system and often lead to attenuated response to novel pathogens and vaccination. In present study we employed 5′RACE UMI-based full length and nearly error-free immunoglobulin profiling to compare plasma cell antibody repertoires in young (19–26 years) and middle-age (45–58 years) individuals vaccinated with a live yellow fever vaccine, modeling a newly encountered pathogen. Our analysis has revealed age-related differences in the responding antibody repertoire ranging from distinct IGH CDR3 repertoire properties to differences in somatic hypermutation intensity and efficiency and antibody lineage tree structure. Overall, our findings suggest that younger individuals respond with a more diverse antibody repertoire and employ a more efficient somatic hypermutation process than elder individuals in response to a newly encountered pathogen.

Human aging is associated with a profound loss of thymus productivity, yet naïve T lymphocytes still maintain their numbers by division in the periphery for many years. The extent of such proliferation may depend on the cytokine environment, including IL-7 and T-cell receptor (TCR) “tonic” signaling mediated by self pMHCs recognition. Additionally, intrinsic properties of distinct subpopulations of naïve T cells could influence the overall dynamics of aging-related changes within the naïve T cell compartment. Here, we investigated the differences in the architecture of TCR beta repertoires for naïve CD4, naïve CD8, naïve CD4+CD25−CD31+ (enriched with recent thymic emigrants, RTE), and mature naïve CD4+CD25−CD31− peripheral blood subsets between young and middle-age/old healthy individuals. In addition to observing the accumulation of clonal expansions (as was shown previously), we reveal several notable changes in the characteristics of T cell repertoire. We observed significant decrease of CDR3 length, NDN insert, and number of non-template added N nucleotides within TCR beta CDR3 with aging, together with a prominent change of physicochemical properties of the central part of CDR3 loop. These changes were similar across CD4, CD8, RTE-enriched, and mature CD4 subsets of naïve T cells, with minimal or no difference observed between the latter two subsets for individuals of the same age group. We also observed an increase in “publicity” (fraction of shared clonotypes) of CD4, but not CD8 naïve T cell repertoires. We propose several explanations for these phenomena built upon previous studies of naïve T-cell homeostasis, and call for further studies of the mechanisms causing the observed changes and of consequences of these changes in respect of the possible holes formed in the landscape of naïve T cell TCR repertoire.

Unique molecular identifiers (UMIs) show outstanding performance in targeted high-throughput resequencing, being the most promising approach for the accurate identification of rare variants in complex DNA samples. This approach has application in multiple areas, including cancer diagnostics, thus demanding dedicated software and algorithms. Here we introduce MAGERI, a computational pipeline that efficiently handles all caveats of UMI-based analysis to obtain high-fidelity mutation profiles and call ultra-rare variants. Using an extensive set of benchmark datasets including gold-standard biological samples with known variant frequencies, cell-free DNA from tumor patient blood samples and publicly available UMI-encoded datasets we demonstrate that our method is both robust and efficient in calling rare variants. The versatility of our software is supported by accurate results obtained for both tumor DNA and viral RNA samples in datasets prepared using three different UMI-based protocols.

Despite the growing number of immune repertoire sequencing studies, the field still lacks software for analysis and comprehension of this high-dimensional data. Here we report VDJtools, a complementary software suite that solves a wide range of T cell receptor (TCR) repertoires post-analysis tasks, provides a detailed tabular output and publication-ready graphics, and is built on top of a flexible API. Using TCR datasets for a large cohort of unrelated healthy donors, twins, and multiple sclerosis patients we demonstrate that VDJtools greatly facilitates the analysis and leads to sound biological conclusions. VDJtools software and documentation are available at https://github.com/mikessh/vdjtools.

BackgroundDuchenne muscular dystrophy (DMD) is the most common inherited muscle disease, leading to severe disability and death in young men. Death is caused by the progressive degeneration of striated muscles aggravated by sterile inflammation. The pleiotropic effects of the mutant gene also include cognitive and behavioral impairments and low bone density. Current interventions in DMD are palliative only as no treatment improves the long-term outcome. Therefore, approaches with a translational potential should be investigated, and key abnormalities downstream from the absence of the DMD product, dystrophin, appear to be strong therapeutic targets. We and others have demonstrated that DMD mutations alter ATP signaling and have identified P2RX7 purinoceptor up-regulation as being responsible for the death of muscles in the mdx mouse model of DMD and human DMD lymphoblasts. Moreover, the ATP-P2RX7 axis, being a crucial activator of innate immune responses, can contribute to DMD pathology by stimulating chronic inflammation. We investigated whether ablation of P2RX7 attenuates the DMD model mouse phenotype to assess receptor suitability as a therapeutic target.Methods and findingsUsing a combination of molecular, histological, and biochemical methods and behavioral analyses in vivo we demonstrate, to our knowledge for the first time, that genetic ablation of P2RX7 in the DMD model mouse produces a widespread functional attenuation of both muscle and non-muscle symptoms. In dystrophic muscles at 4 wk there was an evident recovery in key functional and molecular parameters such as improved muscle structure (minimum Feret diameter, p < 0.001), increased muscle strength in vitro (p < 0.001) and in vivo (p = 0.012), and pro-fibrotic molecular signatures. Serum creatine kinase (CK) levels were lower (p = 0.025), and reduced cognitive impairment (p = 0.006) and bone structure alterations (p < 0.001) were also apparent. Reduction of inflammation and fibrosis persisted at 20 mo in leg (p = 0.038), diaphragm (p = 0.042), and heart muscles (p < 0.001). We show that the amelioration of symptoms was proportional to the extent of receptor depletion and that improvements were observed following administration of two P2RX7 antagonists (CK, p = 0.030 and p = 0.050) without any detectable side effects. However, approaches successful in animal models still need to be proved effective in clinical practice.ConclusionsThese results are, to our knowledge, the first to establish that a single treatment can improve muscle function both short and long term and also correct cognitive impairment and bone loss in DMD model mice. The wide-ranging improvements reflect the convergence of P2RX7 ablation on multiple disease mechanisms affecting skeletal and cardiac muscles, inflammatory cells, brain, and bone. Given the impact of P2RX7 blockade in the DMD mouse model, this receptor is an attractive target for translational research: existing drugs with established safety records could potentially be repurposed for treatment of this lethal disease.

Retroelement activity is a common source of polymorphisms in human genome. The mechanism whereby retroelements contribute to the intraindividual genetic heterogeneity by inserting into the DNA of somatic cells is gaining increasing attention. Brain tissues are suspected to accumulate genetic heterogeneity as a result of the retroelements somatic activity. This study aims to expand our understanding of the role retroelements play in generating somatic mosaicism of neural tissues. Whole-genome Alu and L1 profiling of genomic DNA extracted from the cerebellum, frontal cortex, subventricular zone, dentate gyrus, and the myocardium revealed hundreds of somatic insertions in each of the analyzed tissues. Interestingly, the highest concentration of such insertions was detected in the dentate gyrus-the hotspot of adult neurogenesis. Insertions of retroelements and their activity could produce genetically diverse neuronal subsets, which can be involved in hippocampal-dependent learning and memory.

Reciprocal chromosomal translocations (RCTs) leading to the formation of fusion genes are important drivers of hematological cancers. Although the general requirements for breakage and fusion are fairly well understood, quantitative support for a general mechanism of RCT formation is still lacking. The aim of this paper is to analyze available high-throughput datasets with computational and robust statistical methods, in order to identify genomic hallmarks of translocation partner genes (TPGs). Our results show that fusion genes are generally overexpressed due to increased promoter activity of 5' TPGs and to more stable 3'-UTR regions of 3' TPGs. Furthermore, expression profiling of 5' TPGs and of interaction partners of 3' TPGs indicates that these features can help to explain tissue specificity of hematological translocations. Analysis of protein domains retained in fusion proteins shows that the co-occurrence of specific domain combinations is non-random and that distinct functional classes of fusion proteins tend to be associated with different components of the gene fusion network. This indicates that the configuration of fusion proteins plays an important role in determining which 5' and 3' TPGs will combine in specific fusion genes. It is generally accepted that chromosomal proximity in the nucleus can explain the specific pairing of 5' and 3' TPGS and the recurrence of hematological translocations. Using recently available data for chromosomal contact probabilities (Hi-C) we show that TPGs are preferentially located in early replicated regions and occupy distinct clusters in the nucleus. However, our data suggest that, in general, nuclear position of TPGs in hematological cancers explains neither TPG pairing nor clinical frequency. Taken together, our results support a model in which genomic features related to regulation of expression and replication timing determine the set of candidate genes more likely to be translocated in hematological tissues, with functional constraints being responsible for specific gene combinations.