Search

T-cell receptor (TR) diversity of the variable domains is generated by recombination of both the alpha (TRA) and beta (TRB) chains. The textbook process of TRB chain production starts with TRBD and TRBJ gene rearrangement, followed by the rearrangement of a TRBV gene to the partially rearranged D-J gene. Unsuccessful V-D-J TRB rearrangements lead to apoptosis of the cell. Here, we performed deep sequencing of the poorly explored pool of partial TRBD1-TRBD2 rearrangements in T-cell genomic DNA. We reconstructed full repertoires of human partial TRBD1-TRBD2 rearrangements using novel sequencing and validated them by detecting V-D-J recombination-specific byproducts: excision circles containing the recombination signal (RS) joint 5'D2-RS - 3'D1-RS. Identified rearrangements were in compliance with the classical 12/23 rule, common for humans, rats, and mice and contained typical V-D-J recombination footprints. Interestingly, we detected a bimodal distribution of D-D junctions indicating two active recombination sites producing long and short D-D rearrangements. Long TRB D-D rearrangements with two D-regions are coding joints D1-D2 remaining classically on the chromosome. The short TRB D-D rearrangements with no D-region are signal joints, the coding joint D1-D2 being excised from the chromosome. They both contribute to the TRB V-(D)-J combinatorial diversity. Indeed, short D-D rearrangements may be followed by direct V-J2 recombination. Long D-D rearrangements may recombine further with J2 and V genes forming partial D1-D2-J2 and then complete V-D1-D2-J2 rearrangement. Productive TRB V-D1-D2-J2 chains are present and expressed in thousands of clones of human antigen-experienced memory T cells proving their capacity for antigen recognition and actual participation in the immune response., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The sequences of primers specific for D-D and D-J rearrangements used in this study are intellectual property of MiLaboratories LLC., USA., (Copyright © 2023 Smirnova, Miroshnichenkova, Belyaeva, Kelmanson, Lebedev, Mamedov, Chudakov and Komkov.)

High-throughput sequencing of adaptive immune receptor repertoires is a valuable tool for receiving insights in adaptive immunity studies. Several powerful TCR/BCR repertoire reconstruction and analysis methods have been developed in the past decade. However, detecting and correcting the discrepancy between real and experimentally observed lymphocyte clone frequencies are still challenging. Here, we discovered a hallmark anomaly in the ratio between read count and clone count-based frequencies of non-functional clonotypes in multiplex PCR-based immune repertoires. Calculating this anomaly, we formulated a quantitative measure of V- and J-genes frequency bias driven by multiplex PCR during library preparation called Over Amplification Rate (OAR). Based on the OAR concept, we developed an original software for multiplex PCR-specific bias evaluation and correction named iROAR: immune Repertoire Over Amplification Removal (https://github.com/smiranast/iROAR). The iROAR algorithm was successfully tested on previously published TCR repertoires obtained using both 5' RACE (Rapid Amplification of cDNA Ends)-based and multiplex PCR-based approaches and compared with a biological spike-in-based method for PCR bias evaluation. The developed approach can increase the accuracy and consistency of repertoires reconstructed by different methods making them more applicable for comparative analysis., Competing Interests: AS, AM, YO, MM, YL, DC, IM, AK No competing interests declared, (© 2023, Smirnova et al.)

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., Competing Interests: EL is a member of the Speakers bureau of Janssen. TK is a member of Speakers bureau of: Pfizer, MSD, Novartis, AbbVie, Janssen, Lilly, Celgene, JSC BIOCAD, and Novartis-Sandoz. SE is a member of the Speaker’s bureau of KRKKA, MCB and JSC BIOCAD. SL and IZ provide scientific advice for JSC BIOCAD. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Komech, Koltakova, Barinova, Minervina, Salnikova, Shmidt, Korotaeva, Loginova, Erdes, Bogdanova, Shugay, Lukyanov, Lebedev and Zvyagin.)

The development and implementation of vaccines have been growing exponentially, remaining one of the major successes of healthcare over the last century. Nowadays, active regular immunizations prevent epidemics of many viral diseases, including tick-borne encephalitis (TBE). Along with the generation of virus-specific antibodies, a highly effective vaccine should induce T cell responses providing long-term immune defense. In this study, we performed longitudinal high-throughput T cell receptor (TCR) sequencing to characterize changes in individual T cell repertoires of 11 donors immunized with an inactivated TBE vaccine. After two-step immunization, we found significant clonal expansion of both CD4 + and CD8 + T cells, ranging from 302 to 1706 vaccine-associated TCRβ clonotypes in different donors. We detected several waves of T cell clonal expansion generated by distinct groups of vaccine-responding clones. Both CD4 + and CD8 + vaccine-responding T cell clones formed 17 motifs in TCRβ sequences shared by donors with identical HLA alleles. Our results indicate that TBE vaccination leads to a robust T cell response due to the production of a variety of T cell clones with a memory phenotype, which recognize a large set of epitopes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Sycheva, Komech, Pogorelyy, Minervina, Urazbakhtin, Salnikova, Vorovitch, Kopantzev, Zvyagin, Komkov, Mamedov and Lebedev.)

T cells that express CD56 in peripheral blood of healthy humans represent a heterogeneous and poorly studied subset. In this work, we analyzed this subset for NKG2C expression. In both CD56 + and CD56 - subsets most of the NKG2C + T cells had a phenotype of highly differentiated CD8 + TEMRA cells. The CD56 + NKG2C + T cells also expressed a number of NK cell receptors, such as NKG2D, CD16, KIR2DL2/DL3, and maturation marker CD57 more often than the CD56 - NKG2C + CD3 + cells. TCR β-chain repertoire of the CD3 + CD56 + NKG2C + cell fraction was limited by the prevalence of one or several clonotypes which can be found within the most abundant clonotypes in total or CD8 + T cell fraction TCRβ repertoire. Thus, NKG2C expression in highly differentiated CD56 + T cells was associated with the most expanded αβ T cell clones. NKG2C + T cells produced almost no IFN-γ in response to stimulation with HCMV pp65-derived peptides. This may be partially due to the high content of CD45RA + CD57 + cells in the fraction. CD3 + NKG2C + cells showed signs of activation, and the frequency of this T-cell subset in HCMV-positive individuals was positively correlated with the frequency of NKG2C + NK cells that may imply a coordinated in a certain extent development of the NKG2C + T and NK cell subsets under HCMV infection., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Kovalenko, Zvyagin, Streltsova, Mikelov, Erokhina, Telford, Sapozhnikov and Lebedev.)

COVID-19 is a global pandemic caused by the SARS-CoV-2 coronavirus. T cells play a key role in the adaptive antiviral immune response by killing infected cells and facilitating the selection of virus-specific antibodies. However, neither the dynamics and cross-reactivity of the SARS-CoV-2-specific T-cell response nor the diversity of resulting immune memory is well understood. In this study, we use longitudinal high-throughput T-cell receptor (TCR) sequencing to track changes in the T-cell repertoire following two mild cases of COVID-19. In both donors, we identified CD4 + and CD8 + T-cell clones with transient clonal expansion after infection. We describe characteristic motifs in TCR sequences of COVID-19-reactive clones and show preferential occurrence of these motifs in publicly available large dataset of repertoires from COVID-19 patients. We show that in both donors, the majority of infection-reactive clonotypes acquire memory phenotypes. Certain T-cell clones were detected in the memory fraction at the pre-infection time point, suggesting participation of pre-existing cross-reactive memory T cells in the immune response to SARS-CoV-2., Competing Interests: AM, EK, AT, MB, ER, IM, AF, GE, DC, TM, YL, MP No competing interests declared, AW Senior editor, eLife, (© 2021, Minervina et al.)

Background: Retroelements (REs) occupy a significant part of all eukaryotic genomes including humans. The majority of retroelements in the human genome are inactive and unable to retrotranspose. Dozens of active copies are repressed in most normal tissues by various cellular mechanisms. These copies can become active in normal germline and brain tissues or in cancer, leading to new retroposition events. The consequences of such events and their role in normal cell functioning and carcinogenesis are not yet fully understood. If new insertions occur in a small portion of cells they can be found only with the use of specific methods based on RE enrichment and high-throughput sequencing. The downside of the high sensitivity of such methods is the presence of various artifacts imitating real insertions, which in many cases cannot be validated due to lack of the initial template DNA. For this reason, adequate assessment of rare (< 1%) subclonal cancer specific RE insertions is complicated., Results: Here we describe a new copy-capture technique which we implemented in a method called SeqURE for Sequencing Unknown of Retroposition Events that allows for efficient and reliable identification of new genomic RE insertions. The method is based on the capture of copies of target molecules (copy-capture), selective amplification and sequencing of genomic regions adjacent to active RE insertions from both sides. Importantly, the template genomic DNA remains intact and can be used for validation experiments. In addition, we applied a novel system for testing method sensitivity and precisely showed the ability of the developed method to reliably detect insertions present in 1 out of 100 cells and a substantial portion of insertions present in 1 out of 1000 cells. Using advantages of the method we showed the absence of somatic Alu insertions in colorectal cancer samples bearing tumor-specific L1HS insertions., Conclusions: This study presents the first description and implementation of the copy-capture technique and provides the first methodological basis for the quantitative assessment of RE insertions present in a small portion of cells.

Background and Aims: Intestinal inflammation in inflammatory bowel diseases [IBD] is thought to be T cell mediated and therefore dependent on the interaction between the T cell receptor [TCR] and human leukocyte antigen [HLA] proteins expressed on antigen presenting cells. The collection of all TCRs in one individual, known as the TCR repertoire, is characterised by enormous diversity and inter-individual variability. It was shown that healthy monozygotic [MZ] twins are more similar in their TCR repertoire than unrelated individuals. Therefore MZ twins, concordant or discordant for IBD, may be useful to identify disease-related and non-genetic factors in the TCR repertoire which could potentially be used as disease biomarkers., Methods: Employing unique molecular barcoding that can distinguish between polymerase chain reaction [PCR] artefacts and true sequence variation, we performed deep TCRα and TCRβ repertoire profiling of the peripheral blood of 28 MZ twin pairs from Denmark and Germany, 24 of whom were discordant and four concordant for IBD., Results: We observed disease- and smoking-associated traits such as sharing, diversity and abundance of specific clonotypes in the TCR repertoire of IBD patients, and particularly in patients with active disease, compared with their healthy twins., Conclusions: Our findings identified TCR repertoire features specific for smokers and IBD patients, particularly when signs of disease activity were present. These findings are a first step towards the application of TCR repertoire analyses as a valuable tool to characterise inflammatory bowel diseases and to identify potential biomarkers and true disease causes., (© European Crohn’s and Colitis Organisation (ECCO) 2019.)

The diverse repertoire of T-cell receptors (TCR) plays a key role in the adaptive immune response to infections. Using TCR alpha and beta repertoire sequencing for T-cell subsets, as well as single-cell RNAseq and TCRseq, we track the concentrations and phenotypes of individual T-cell clones in response to primary and secondary yellow fever immunization - the model for acute infection in humans - showing their large diversity. We confirm the secondary response is an order of magnitude weaker, albeit ∼10 days faster than the primary one. Estimating the fraction of the T-cell response directed against the single immunodominant epitope, we identify the sequence features of TCRs that define the high precursor frequency of the two major TCR motifs specific for this particular epitope. We also show the consistency of clonal expansion dynamics between bulk alpha and beta repertoires, using a new methodology to reconstruct alpha-beta pairings from clonal trajectories., Competing Interests: AM, MP, EK, VK, PB, ER, AF, DC, IM, YL, TM No competing interests declared, AW Senior editor, eLife, (© 2020, Minervina et al.)

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., Competing Interests: The authors have declared that no competing interests exist.

Retroelements are considered as one of the important sources of genomic variability in modern humans. It is known that transposition activity of retroelements in germline cells generates new insertions in various genomic loci and sometimes results in genetic diseases. Retroelements activity in somatic cells is restricted by different cellular mechanisms; however, there is an evidence for it in some tissue types. Somatic insertions can trigger tumorigenesis or participate in normal functioning such as generation of neurons' plasticity. In spite of the rapid development of high-throughput sequencing methods a confident detection of somatic insertions is still quite a challenging task. That, in part, is due to the absence of adequate bioinformatic tools for the analysis of sequencing data. Here, we propose an advanced computational pipeline for the identification of somatic insertions in datasets generated by selective amplification and high-throughput sequencing of genomic regions flanking insertions of AluYa5. Particular attention is paid for the identification of various artifacts arising in course of library preparation and the parameters for their filtration. Pipeline sensitivity is confirmed by in silico experiments with artificial datasets. Using the proposed pipeline we remove at least 80% of artifacts and preserve 75% of potentially somatic insertions. The approaches used in this work can be applied for the study of other mobile elements insertion variability.

T cell receptor (TCR) repertoire data contain information about infections that could be used in disease diagnostics and vaccine development, but extracting that information remains a major challenge. Here we developed a statistical framework to detect TCR clone proliferation and contraction from longitudinal repertoire data. We applied this framework to data from three pairs of identical twins immunized with the yellow fever vaccine. We identified 600 to 1,700 responding TCRs in each donor and validated them using three independent assays. While the responding TCRs were mostly private, albeit with higher overlap between twins, they could be well-predicted using a classifier based on sequence similarity. Our method can also be applied to samples obtained postinfection, making it suitable for systematic discovery of new infection-specific TCRs in the clinic., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

Background: There is increasing evidence that the transpositional activity of retroelements (REs) is not limited to germ line cells, but often occurs in tumor and normal somatic cells. Somatic transpositions were found in several human tissues and are especially typical for the brain. Several computational and experimental approaches for detection of somatic retroelement insertions was developed in the past few years. These approaches were successfully applied to detect somatic insertions in clonally expanded tumor cells. At the same time, identification of somatic insertions presented in small proportion of cells, such as neurons, remains a considerable challenge., Results: In this study, we developed a normalization procedure for library enrichment by DNA sequences corresponding to rare somatic RE insertions. Two rounds of normalization increased the number of fragments adjacent to somatic REs in the sequenced sample by more than 26-fold, and the number of identified somatic REs was increased by 8-fold., Conclusions: The developed technique can be used in combination with vast majority of modern RE identification approaches and can dramatically increase their capacity to detect rare somatic RE insertions in different types of cells., Competing Interests: The study was approved by the local ethics committee and conducted in accordance with the Declaration of Helsinki.Not applicable.The authors declare that they have no competing interests.Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Objective: The risk of AS is associated with genomic variants related to antigen presentation and specific cytokine signalling pathways, suggesting the involvement of cellular immunity in disease initiation/progression. The aim of the present study was to explore the repertoire of TCR sequences in healthy donors and AS patients to uncover AS-linked TCR variants., Methods: Using quantitative molecular-barcoded 5'-RACE, we performed deep TCR β repertoire profiling of peripheral blood (PB) and SF samples for 25 AS patients and 108 healthy donors. AS-linked TCR variants were identified using a new computational approach that relies on a probabilistic model of the VDJ rearrangement process., Results: Using the donor-agnostic probabilistic model, we reveal a TCR β motif characteristic for PB of AS patients, represented by eight highly homologous amino acid sequence variants. Some of these variants were previously reported in SF and PB of patients with ReA and in PB of AS patients. We demonstrate that identified AS-linked clones have a CD8+ phenotype, present at relatively low frequencies in PB, and are significantly enriched in matched SF samples of AS patients., Conclusion: Our results suggest the involvement of a particular antigen-specific subset of CD8+ T cells in AS pathogenesis, confirming and expanding earlier findings. The high similarity of the clonotypes with the ones found in ReA implies common mechanisms for the initiation of the diseases.

Autologous hematopoietic stem cell transplantation (HSCT), a safer type of HSCT than allogeneic HSCT, is a promising therapy for patients with severe autoimmune diseases (ADs). Despite the long history of medical practice, structural changes in the adaptive immune system as a result of autologous HSCT in patients with various types of ADs remain poorly understood. In this study, we used high-throughput sequencing to investigate the structural changes in the peripheral blood T-cell repertoire in adult patients with ankylosing spondylitis (AS) during two years after autologous HSCT. The implementation of unique molecular identifiers allowed us to substantially reduce the impact of the biases occurring during the preparation of libraries, to carry out a comparative analysis of the various properties of the T-cell repertoire between different time points, and to track the dynamics of both distinct T-cell clonotypes and T-cell subpopulations. In the first year of the reconstitution, clonal diversity of the T-cell repertoire remained lower than the initial one in both patients. During the second year after HSCT, clonal diversity continued to increase and reached a normal value in one of the patients. The increase in the diversity was associated with the emergence of a large number of low-frequency clonotypes, which were not identified before HSCT. Efficiency of clonotypes detection after HSCT was dependent on their abundance in the initial repertoire. Almost all of the 100 most abundant clonotypes observed before HSCT were detected 2 years after transplantation and remained highly abundant irrespective of their CD4+ or CD8+ phenotype. A total of up to 25% of peripheral blood T cells 2 years after HSCT were represented by clonotypes from the initial repertoire.

Vaccination against influenza is widely used to protect against seasonal flu epidemic although its effectiveness is debated. Here we performed deep quantitative T cell receptor repertoire profiling in peripheral blood of a healthy volunteer in response to trivalent subunit influenza vaccine. We did not observe significant rebuilding of peripheral blood T cell receptors composition in response to vaccination. However, we found several clonotypes in memory T cell fraction that were undetectable before the vaccination and had a maximum concentration at day 45 after vaccine administration. These cells were found in lower concentration in the course of repertoire monitoring for two years period. Our observation suggests a potential for recruitment of only a limited number of new T cells after each seasonal influenza vaccination., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Diverse repertoires of hypervariable immunoglobulin receptors (TCR and BCR) recognize antigens in the adaptive immune system. The development of immunoglobulin receptor repertoire sequencing methods makes it possible to perform repertoire-wide disease association studies of antigen receptor sequences. We developed a statistical framework for associating receptors to disease from only a small cohort of patients, with no need for a control cohort. Our method successfully identifies previously validated Cytomegalovirus and type one diabetes responsive TCR[Formula: see text] sequences ., Competing Interests: MP, AM, DC, IM, YL, TM No competing interests declared, AW Reviewing editor, eLife, (© 2018, Pogorelyy et al.)

The diversity of T-cell receptors recognizing foreign pathogens is generated through a highly stochastic recombination process, making the independent production of the same sequence rare. Yet unrelated individuals do share receptors, which together constitute a "public" repertoire of abundant clonotypes. The TCR repertoire is initially formed prenatally, when the enzyme inserting random nucleotides is downregulated, producing a limited diversity subset. By statistically analyzing deep sequencing T-cell repertoire data from twins, unrelated individuals of various ages, and cord blood, we show that T-cell clones generated before birth persist and maintain high abundances in adult organisms for decades, slowly decaying with age. Our results suggest that large, low-diversity public clones are created during pre-natal life, and survive over long periods, providing the basis of the public repertoire.

αβT-cell-depleted allogeneic hematopoietic cell transplantation holds promise for the safe and accessible therapy of both malignant and non-malignant blood disorders. Here we employed molecular barcoding normalized T-cell receptor (TCR) profiling to quantitatively track T-cell immune reconstitution after TCRαβ-/CD19-depleted transplantation in children. We demonstrate that seemingly early reconstitution of αβT-cell counts 2 months after transplantation is based on only several hundred rapidly expanded clones originating from non-depleted graft cells. In further months, frequency of these hyperexpanded clones declines, and after 1 year the observed T-cell counts and TCRβ diversity are mostly provided by the newly produced T cells. We also demonstrate that high TCRβ diversity at day 60 observed for some of the patients is determined by recipient T cells and intrathymic progenitors that survived conditioning regimen. Our results indicate that further efforts on optimization of TCRαβ-/CD19-depleted transplantation protocols should be directed toward providing more efficient T-cell defense in the first months after transplantation.

T cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignant neoplasm of the lymphocyte precursors that suffered malignant transformation arresting the lymphoid cell differentiation. Clinical studies revealed monoor, more rarely, oligoclonal nature of the disease. A precise identification of malignant clone markers is both the crucial stage of early diagnostics and the essential prognostic factor for therapeutic treatment. Here we present an improved system for unbiased detection of lymphoblastic clones in bone marrow aspirates of T-ALL patients. The system based on multiplex PCR of rearranged T-cell receptor locus (TRB) and straightforward sequencing of the resulted PCR fragments. Testing of the system on genomic DNA from Jurkat cell line and four clinical bone marrow aspirates revealed a set of unique TRB rearrangements that precisely characterize each of tested samples. Therefore, the outcome of the system produces highly informative molecular genetic markers for further monitoring of minimal residual disease in T-ALL patients.

Emerging high-throughput sequencing methods for the analyses of complex structure of TCR and BCR repertoires give a powerful impulse to adaptive immunity studies. However, there are still essential technical obstacles for performing a truly quantitative analysis. Specifically, it remains challenging to obtain comprehensive information on the clonal composition of small lymphocyte populations, such as Ag-specific, functional, or tissue-resident cell subsets isolated by sorting, microdissection, or fine needle aspirates. In this study, we report a robust approach based on unique molecular identifiers that allows profiling Ag receptors for several hundred to thousand lymphocytes while preserving qualitative and quantitative information on clonal composition of the sample. We also describe several general features regarding the data analysis with unique molecular identifiers that are critical for accurate counting of starting molecules in high-throughput sequencing applications., (Copyright © 2015 by The American Association of Immunologists, Inc.)

Background: The Immunoglobulins (IG) and the T cell receptors (TR) play the key role in antigen recognition during the adaptive immune response. Recent progress in next-generation sequencing technologies has provided an opportunity for the deep T cell receptor repertoire profiling. However, a specialised software is required for the rational analysis of massive data generated by next-generation sequencing., Results: Here we introduce tcR, a new R package, representing a platform for the advanced analysis of T cell receptor repertoires, which includes diversity measures, shared T cell receptor sequences identification, gene usage statistics computation and other widely used methods. The tool has proven its utility in recent research studies., Conclusions: tcR is an R package for the advanced analysis of T cell receptor repertoires after primary TR sequences extraction from raw sequencing reads. The stable version can be directly installed from The Comprehensive R Archive Network ( http://cran.r-project.org/mirrors.html ). The source code and development version are available at tcR GitHub ( http://imminfo.github.io/tcr/ ) along with the full documentation and typical usage examples.

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.

Adaptive immunity in humans is provided by hypervariable Ig-like molecules on the surface of B and T cells. The final set of these molecules in each organism is formed under the influence of two forces: individual genetic traits and the environment, which includes the diverse spectra of alien and self-antigens. Here we assess the impact of individual genetic factors on the formation of the adaptive immunity by analyzing the T-cell receptor (TCR) repertoires of three pairs of monozygous twins by next-generation sequencing. Surprisingly, we found that an overlap between the TCR repertoires of monozygous twins is similar to an overlap between the TCR repertoires of nonrelated individuals. However, the number of identical complementary determining region 3 sequences in two individuals is significantly increased for twin pairs in the fraction of highly abundant TCR molecules, which is enriched by the antigen-experienced T cells. We found that the initial recruitment of particular TCR V genes for recombination and subsequent selection in the thymus is strictly determined by individual genetic factors. J genes of TCRs are selected randomly for recombination; however, the subsequent selection in the thymus gives preference to some α but not β J segments. These findings provide a deeper insight into the mechanism of TCR repertoire generation.

The decrease of TCR diversity with aging has never been studied by direct methods. In this study, we combined high-throughput Illumina sequencing with unique cDNA molecular identifier technology to achieve deep and precisely normalized profiling of TCR β repertoires in 39 healthy donors aged 6-90 y. We demonstrate that TCR β diversity per 10(6) T cells decreases roughly linearly with age, with significant reduction already apparent by age 40. The percentage of naive T cells showed a strong correlation with measured TCR diversity and decreased linearly up to age 70. Remarkably, the oldest group (average age 82 y) was characterized by a higher percentage of naive CD4(+) T cells, lower abundance of expanded clones, and increased TCR diversity compared with the previous age group (average age 62 y), suggesting the influence of age selection and association of these three related parameters with longevity. Interestingly, cross-analysis of individual TCR β repertoires revealed a set >10,000 of the most representative public TCR β clonotypes, whose abundance among the top 100,000 clones correlated with TCR diversity and decreased with aging.

The relationship between maternal and child immunity has been actively studied in the context of complications during pregnancy, autoimmune diseases, and haploidentical transplantation of hematopoietic stem cells and solid organs. Here, we have for the first time used high-throughput Illumina HiSeq sequencing to perform deep quantitative profiling of T cell receptor (TCR) repertoires for peripheral blood samples of three mothers and their six children. Advanced technology allowed accurate identification of 5 × 10(5) to 2 × 10(6) TCR beta clonotypes per individual. We performed comparative analysis of these TCR repertoires with the aim of revealing characteristic features that distinguish related mother-child pairs, such as relative TCR beta variable segment usage frequency and relative overlap of TCR beta complementarity-determining region 3 (CDR3) repertoires. We show that thymic selection essentially and similarly shapes the initial output of the TCR recombination machinery in both related and unrelated pairs, with minor effect from inherited differences. The achieved depth of TCR profiling also allowed us to test the hypothesis that mature T cells transferred across the placenta during pregnancy can expand and persist as functional microchimeric clones in their new host, using characteristic TCR beta CDR3 variants as clonal identifiers.

High-throughput sequencing has the power to reveal the nature of adaptive immunity as represented by the full complexity of T-cell receptor (TCR) and antibody (IG) repertoires, but is at present severely compromised by the quantitative bias, bottlenecks, and accumulated errors that inevitably occur in the course of library preparation and sequencing. Here we report an optimized protocol for the unbiased preparation of TCR and IG cDNA libraries for high-throughput sequencing, starting from thousands or millions of live cells in an investigated sample. Critical points to control are revealed, along with tips that allow researchers to minimize quantitative bias, accumulated errors, and cross-sample contamination at each stage, and to enhance the subsequent bioinformatic analysis. The protocol is simple, reliable, and can be performed in 1-2 days.

The TCR repertoire is a mirror of the human immune system that reflects processes caused by infections, cancer, autoimmunity, and aging. Next generation sequencing (NGS) is becoming a powerful tool for deep TCR profiling; yet, questions abound regarding the methodological approaches for sample preparation and correct data interpretation. Accumulated PCR and sequencing errors along with library preparation bottlenecks and uneven PCR efficiencies lead to information loss, biased quantification, and generation of huge artificial TCR diversity. Here, we compare Illumina, 454, and Ion Torrent platforms for individual TCR profiling, evaluate the rate and character of errors, and propose advanced platform-specific algorithms to correct massive sequencing data. These developments are applicable to a wide variety of next generation sequencing applications. We demonstrate that advanced correction allows the removal of the majority of artificial TCR diversity with concomitant rescue of most of the sequencing information. Thus, this correction enhances the accuracy of clonotype identification and quantification as well as overall TCR diversity measurements., (© 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Autologous haematopoietic stem cell transplantation is highly efficient for the treatment of systemic autoimmune diseases, but its consequences for the immune system remain poorly understood. Here, we describe an optimized RNA-based technology for unbiased amplification of T cell receptor beta-chain libraries and use it to perform the first detailed, quantitative tracking of T cell clones during 10 months after transplantation. We show that multiple clones survive the procedure, contribute to the immune response to activated infections, and form a new skewed and stable T cell receptor repertoire., (Copyright © 2011 EMBO Molecular Medicine.)

Increasing evidence points to a role for killer immunoglobulin-like receptors (KIRs) in the development of autoimmune diseases. In particular, a positive association of KIR3DS1 (activating receptor) and a negative association of KIR3DL1 (inhibitory receptor) alleles with ankylosing spondylitis (AS) have been reported by several groups. However, none of the studies analyzed these associations in the context of functionality of polymorphic KIR3DL1. To better understand how the KIR3DL1/3DS1 genes determine susceptibility to AS, we analyzed the frequencies of alleles and genotypes encoding functional (KIR3DL1*F) and non-functional (KIR3DL1*004) receptors. We genotyped 83 AS patients and 107 human leukocyte antigen (HLA)-B27-positive healthy controls from the Russian Caucasian population using a two-stage sequence-specific primer PCR, which distinguishes KIR3DS1, KIR3DL1*F and KIR3DL1*004 alleles. For the patients carrying two functional KIR3DL1 alleles, those alleles were additionally genotyped to identify KIR3DL1*005 and KIR3DL1*007 alleles, which are functional but are expressed at low levels. KIR3DL1 was negatively associated with AS at the expense of KIR3DL1*F but not of KIR3DL1*004. This finding indicates that the inhibitory KIR3DL1 receptor protects against the development of AS and is not simply a passive counterpart of the segregating KIR3DS1 allele encoding the activating receptor. However, analysis of genotype frequencies indicates that the presence of KIR3DS1 is a more important factor for AS susceptibility than the absence of KIR3DL1*F. The activation of either natural killer (NK) or T cells via the KIR3DS1 receptor can be one of the critical events in AS development, while the presence of the functional KIR3DL1 receptor has a protective effect. Nevertheless, even individuals with a genotype that carried two inhibitory KIR3DL1 alleles expressed at high levels could develop AS.

A number of genetic systems for human genetic identification based on short tandem repeats or single nucleotide polymorphisms are widely used for crime detection, kinship studies and in analysis of victims of mass disasters. Here, we have developed a new set of 32 molecular genetic markers for human genetic identification based on polymorphic retroelement insertions. Allele frequencies were determined in a group of 90 unrelated individuals from four genetically distant populations of the Russian Federation. The mean match probability and probability of paternal exclusion, calculated based on population data, were 5.53 x 10(-14) and 99.784%, respectively. The developed system is cheap and easy to use as compared to all previously published methods. The application of fluorescence-based methods for allele discrimination allows to use the human genetic identification set in automatic and high-throughput formats.

Ankylosing spondylitis (AS) belongs to a group of autoimmune diseases affecting the axial skeleton. Beside thehla-b*27allele, several other human genes that control the variety processes of immune homeostasis are considered to be associated with AS manifestation in different human populations. Among strong associated non-MHC geneserap1 encodingthe endoplasmic reticulum aminopeptidase 1 isoform was recently identified by single nucleotide polymorphisms (SNPs) meta analysis. In our study we inspected the genetic association of five non-synonymous coding SNPs fromerap1 withAS in Caucasians. We implemented the SSP-PCR system for precise genotyping of 87hla-b*27positive AS patients and 77hla-b*27healthy donors from the Russian population. Considerable differences in allele's frequencies within patients vs control cohort were shown for 3 of 5 SNPs under investigation. Using the EM-algorhitm we reconstructed 3-marker haplotypes that distinguish with high probability two cohorts due to differences in the haplotypes frequencies. In such a way both the sensitive, CCT, haplotype and the protective, TTC, one were predicted. To verify the calculation we determined genuine frequencies of 5-marker haplotypes in AS cohort by haplotyping of individual cDNA samples using improved SSP-PCR primer set. We demonstrated that the frequencies ofin silicareconstucted haplotypes and the frequencies of experimentally detected haplotypes are in a good agreement. Frequency of the risk haplotype CCT (rs17482078/10050860/2287987) detected within AS cohort reaches 88%, as well as the frequency calculated by EM-algorhitm.

Ankylosing spondylitis (AS) is commonly characterized by clonal expansions of T cells. However, these clonal populations are poorly studied and their role in disease initiation and progression remains unclear. Here, we performed mass sequencing of TCR V beta libraries to search for the expanded T cell clones for two AS patients. A number of clones comprising more than 5% of the corresponding TCR V beta family were identified in both patients. For the first time, expanded clones were shown to be stably abundant in blood samples of AS patients for the prolonged period (1.5 and 2.5 years for two patients, correspondingly). These clones were individually characterized in respect to their differentiation status using fluorescent cell sorting with CD27, CD28, and CD45RA markers followed by quantitative identification of each clone within corresponding fraction using real time PCR analysis. Stable clones differed in phenotype and several were shown to belong to the proinflammatory CD27 - /CD28 - population. Their potentially cytotoxic status was confirmed by staining with perforin-specific antibodies. Search for the TCR V beta CRD3 sequences homologous to the identified clones revealed close matches with the previously reported T cell clones from AS and reactive arthritis patients, thus supporting their role in the disease and proposing consensus TCR V beta CDR3 motifs for AS. Interestingly, these motifs were also found to have homology with earlier reported virus-specific CDR3 variants, indicating that viral infections could play role in development of AS.

Expressed sequence tags (ESTs) represent 500-1000-bp-long sequences corresponding to mRNAs derived from different sources (cell lines, tissues, etc.). The human EST database contains over 8,000,000 sequences, with over 4,000,000,000 total nucleotides. RNA molecules are transcribed from a genomic DNA template; therefore, all ESTs should match corresponding genomes. Nevertheless, we have found in the human EST database approximately 11,000 ESTs not matching sequences in the human genome database. The presence of "trash" ESTs (TESTs) in the EST database could result from DNA or RNA contamination of the laboratory equipment, tissues, or cell lines. TESTs could also represent sequences from unidentified human genes or from species inhabiting the human body. Here, we attempt to identify the sources of human EST database contaminations. In particular, we discuss systematic contamination of the mammalian EST databases with sequences of plants.

Being the most effectively transposed primate-specific SINEs, Alu elements are present in more than one million copies in the human genome and include most recently transposed subsets of AluY elements that are polymorphic in humans. Although Alu elements are commonly thought to play an essential role in shaping and functioning of primate genomes, the understanding of the impact of recent Alu insertions on human gene expression is far from being comprehensive. Here we compared hnRNA contents for allele pairs of genes heterozygous for AluY insertions in their introns in human cell lines of various origins. We demonstrated that some AluY insertions correlated with decreased content of the corresponding hnRNAs. The effect observed does not depend on sequences of Alu elements and their orientation but is likely to be cell type specific.

LINE-1 (L1) retrotransposons comprise about 17% of the human genome and include a recently transposed set of Ta-L1 elements that are polymorphic in humans. Although it is widely believed that L1s play an essential role in shaping and functioning of mammalian genomes, the understanding of the impact of L1 insertions on gene expression is far from being comprehensive. Here we compared hnRNA contents for allele pairs of genes heterozygous for Ta-L1 insertions in their introns in human cell lines of various origin. We demonstrated that some Ta-L1 insertions correlated with decreased content of the corresponding hnRNAs. This effect was characteristic of only nearly full-sized L1s and seemed to be tissue specific.

Human cell lines are an indispensable tool for functional studies of living entities in their numerous manifestations starting with integral complex systems such as signal pathways and networks, regulation of gene ensembles, epigenetic factors, and finishing with pathological changes and impact of artificially introduced elements, such as various transgenes, on the behavior of the cell. Therefore, it is highly desirable to have reliable cell line identification techniques to make sure that the cell lines to be used in experiments are exactly what is expected. To this end, we developed a set of informative markers based on insertion polymorphism of human retroelements (REs). The set includes 47 pairs of PCR primers corresponding to introns of the human genes with dimorphic LINE1 (L1) and Alu insertions. Using locus-specific PCR assays, we have genotyped 10 human cell lines of various origins. For each of these cell lines, characteristic fingerprints were obtained. An estimated probability that two different cell lines possess the same marker genotype is about 10-18. Therefore, the proposed set of markers provides a reliable tool for cell line identification.

A new experimental technique for genome-wide detection of integration sites of polymorphic retroelements (REs) is described. The technique allows one to reveal the absence of a retroelement in an individual genome provided that this retroelement is present in at least one of several other genomes under comparison. Since quite a number of genomes are compared simultaneously, the search for polymorphic REs insertions is very efficient. The technique includes two whole-genome selective PCR amplifications of sequences flanking REs: one for a particular genome and another one for a mixture of ten different genomes. A subsequent subtractive hybridization of the obtained amplicons with DNA of a particular genome as driver results in isolation of polymorphic insertions. The technique was successfully applied for identification of 41 new polymorphic human AluYa5/Ya8 insertions. Among them, 18 individual Alu elements first sequenced in this work were not found in the available human genome databases. This result suggests that significant part of polymorphic REs were not identified during genome sequencing and remain to be detected and characterized. The proposed method does not depend on preliminary knowledge of evolutionary history of retroelements and can be applied for identification of insertion/deletion polymorphic markers in genomes of different species.

Human endogenous retroviruses (HERVs) make up a substantial part of the human genome. HERVs and solitary long terminal repeats (solo LTRs) are usually flanked by 4-6 nt short direct repeats through the well-known mechanism of their integration. A number of solo LTRs flanked by unusually long direct repeats were detected in the human genome. These unusual structures might be a product of an alternative virus insertion mechanism.

One of the evolutionary mechanisms for acquisition of novel functional sequences can be domestication of exogenous retroviruses that have been integrated into the germ line. The whole genome mapping of such elements in various species could reveal differences in positions of the retroviral integration and suggest possible roles of these differences in speciation. Here, we describe the number, locations and sequence features of the human endogenous retrovirus HERV-K (HML-2) long terminal repeat (LTR) sequences on human chromosome 21. We show that their distribution along the chromosome is not only non-random but also roughly correlated with the gene density. Amplification of orthologous LTR sites from a number of primate genomes produced patterns of presence and absence for each LTR sequence and allowed determination of the phylogenetic ages and evolutionary order of appearance of individual LTRs. The identity level and phylogenetic age of the LTRs did not correlate with their map locations. Thus, despite the non-random distribution of LTRs, they have apparently been inserted randomly into the chromosome relative to each other. As evidenced in previous studies of chromosomes 19 and 22, this is a characteristic of HERV-K integration.

A locus harboring a human endogenous retroviral LTR (long terminal repeat) was mapped on the short arm of human chromosome 7 (7p22), and its evolutionary history was investigated. Sequences of two human genome fragments that were homologous to the LTR-flanking sequences were found in human genome databases: (1) an LTR-containing DNA fragment from region 3p13 of the human genome, which includes clusters of olfactory receptor genes and pseudogenes; and (2) a fragment of region 21q22.1 lacking LTR sequences. PCR analysis demonstrated that LTRs with highly homologous flanking sequences could be found in the genomes of human, chimp, gorilla, and orangutan, but were absent from the genomes of gibbon and New World monkeys. A PCR assay with a primer set corresponding to the sequence from human Chr 3 allowed us to detect LTR-containing paralogous sequences on human chromosomes 3, 4, 7, and 11. The divergence times for the LTR-flanking sequences on chromosomes 3 and 7, and the paralogous sequence on chromosome 21, were evaluated and used to reconstruct the order of duplication events and retroviral insertions. (1) An initial duplication event that occurred 14-17 Mya and before LTR insertion - produced two loci, one corresponding to that located on Chr 21, while the second was the ancestor of the loci on chromosomes 3 and 7. (2) Insertion of the LTR (most probably as a provirus) into this ancestral locus took place 13 Mya. (3) Duplication of the LTR-containing ancestral locus occurred 11 Mya, forming the paralogous modern loci on Chr 3 and 7.

Reporter gene analysis of HERV-K solitary long terminal repeats (LTRs) showed that they retain detectable activity in human teratocarcinoma cells, and can direct the transcription in both orientations relative to the reporter gene. Deletion analysis demonstrated the possible existence of alternative promoters within the LTR as well as a silencer-like element in the U5 region. Our results indicate also that all-trans-retinoic acid is capable of modulating expression of the reporter gene directed by a HERV-K LTR in NT2/D1 cells.

The classification of the long terminal repeats (LTRs) of the human endogenous retrovirus HERV-K (HML-2) family was refined according to diagnostic differences between the LTR sequences. The mutation rate was estimated to be approximately equal for LTRs belonging to different families and branches of human endogenous retroviruses (HERVs). An average mutation rate value was calculated based on differences between LTRs of the same HERV and was found to be 0.13% per million years (Myr). Using this value, the ages of different LTR groups belonging to the LTR HML-2 subfamily were found to vary from 3 to 50Myr. Orthologous potential LTR-containing loci from different primate species were PCR amplified using primers corresponding to the genomic sequences flanking LTR integration sites. This allowed us to calculate the phylogenetic times of LTR integrations in primate lineages in the course of the evolution and to demonstrate that they are in good agreement with the LTR ages calculated from the mutation rates. Human-specific integrations for some very young LTRs were demonstrated. The possibility of LTRs and HERVs involvement in the evolution of primates is discussed.

The availability of bacterial genome sequences raises an important new problem - how can one move from completely sequenced microorganisms as a reference to the hundreds and thousands of other strains or isolates of the same or related species that will not be sequenced in the near future? An efficient way to approach this task is the comparison of genomes by subtractive hybridization. Recently we developed a sensitive and reproducible subtraction procedure for comparison of bacterial genomes, based on the method of suppression subtractive hybridization (SSH). In this work we demonstrate the applicability of subtractive hybridization to the comparison of the related but markedly divergent bacterial species Escherichia coli and Salmonella typhimurium. Clone libraries representing sequence differences were obtained and, in the case of completely sequenced E. coli genome, the differences were directly placed in the genome map. About 60% of the differential clones identified by SSH were present in one of the genomes under comparison and absent from the other. Additional differences in most cases represent sequences that have diverged considerably in the course of evolution. Such an approach to comparative bacterial genomics can be applied both to studies of interspecies evolution - to elucidate the "strategies" that enable different genomes to fit their ecological niches - and to development of diagnostic probes for the rapid identification of pathogenic bacterial species.

A highly polymorphic microsatellite (CA)n-marker (CAct685) previously isolated from human chromosome 19 cosmid library was localized near GPI in 19q13.1. For the fine localization of this marker, the hybridization with chromosome 19-specific cosmid libraries assembled in contigs was used. Polymorphism analysis of the marker in 12 populations of Russia and neighboring countries showed 14 alleles containing from 16 to 30 repeat units. Populations belonging to Indo-European, Uralic and Altaic linguistic families demonstrated a great similarity in allele frequency profiles. Differences between these populations were lower for CAct685 than for classical markers. Allele distribution of CAct685 in a Chukchi population belonging to the Chukchi-Kamchatkan linguistic family differs from those in all other populations, that may be typical for Mongoloid population or reflect an ethnic history of Chukchi as a small population. Thus use of the CAct685 marker seems to be effective for analysis of distant peoples.

A chimeric long terminal repeat (LTR) containing the whole LTR of a human endogenous retrovirus-like element of the H family (HERV-H) inserted downstream of the core enhancer region of the 5' LTR of a HERV-K retroelement was detected and sequenced in the human 19p12 locus, known to be enriched with genes encoding zinc finger proteins. Similar chimeras were also detected in human chromosomes 10 and Y in human-hamster hybrid cells containing individual human chromosomes. This finding was interpreted as evidence of transpositions of the chimera in the genome. PCR analyses detected the chimera in the genomes of chimpanzee and gorilla, but not in that of orangutan. These data demonstrate that the chimera appeared in the primate germ cells more than 10 million years ago, before divergence of the human/chimpanzee and the gorilla lineages. The combination of the two LTRs forms a new regulatory system that can be involved in nearby gene expression.

R-loops formed by short RNA transcripts have been imaged by atomic force microscopy (AFM) at a constant force in the height mode. The technique was applied to mapping the human endogenous retrovirus K10 family (HERV-K10) long terminal repeats (LTR) within individual plasmids and cosmids. RNA probes specific for the U3 (384 nt) and U5 (375 nt) LTR regions separated by a span of 200 bp were used for R-loop formation with LTRs located within plasmid (3.8 kb) or cosmid ( approximately 40 kb) DNAs. R-loops stabilized by glyoxal treatment and adsorbed onto the mica surface in the presence of magnesium ions looked like looped out segments of RNA:DNA hybrids. The total yield of R-loops was usually approximately 95%. The RNA:DNA hybrids were found to be 12-15% shorter than the corresponding DNA:DNA duplex. The two regions of the LTR could be easily discerned in the AFM images as clearly separated loops. R-loop positions determined on cosmids by AFM were accurate to approximately 0.5% of the cosmid length. This technique might be easily adapted for mapping various sequences such as gene exons or regulatory regions and for detecting insertions, deletions and rearrangements that cause human genetic diseases.