MAGERI: Computational pipeline for molecular-barcoded targeted resequencing.

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.

Application of nonsense-mediated primer exclusion (NOPE) for preparation of unique molecular barcoded libraries.

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.

Aberrant adaptive immune response underlies genetic susceptibility to tuberculosis.

Mycobacterium tuberculosis (Mtb) remains a major threat worldwide, although only a fraction of infected individuals develops tuberculosis (TB). TB susceptibility is shaped by multiple genetic factors, and we performed comparative immunological analysis of two mouse strains to uncover relevant mechanisms underlying susceptibility and resistance. C57BL/6 mice are relatively TB-resistant, whereas I/St mice are prone to develop severe TB, partly due to the MHC-II allelic variant that shapes suboptimal CD4+ T cell receptor repertoire. We investigated the repertoires of lung-infiltrating helper T cells and B cells at the progressed stage in both strains. We found that lung CD4+ T cell repertoires of infected C57BL/6 but not I/St mice contained convergent TCR clusters with functionally confirmed Mtb specificity. Transcriptomic analysis revealed a more prominent Th1 signature in C57BL/6, and expression of pro-inflammatory IL-16 in I/St lung-infiltrating helper T cells. The two strains also showed distinct Th2 signatures. Furthermore, the humoral response of I/St mice was delayed, less focused, and dominated by IgG/IgM isotypes, whereas C57BL/6 mice generated more Mtb antigen-focused IgA response. We conclude that the inability of I/St mice to produce a timely and efficient anti-Mtb adaptive immune responses arises from a suboptimal helper T cell landscape that also impacts the humoral response, leading to diffuse inflammation and severe disease.

Local Enrichment with Convergence of Enriched T-Cell Clones Are Hallmarks of Effective Peptide Vaccination against B16 Melanoma.

BACKGROUND: Some peptide anticancer vaccines elicit a strong T-cell memory response but fail to suppress tumor growth. To gain insight into tumor resistance, we compared two peptide vaccines, p20 and p30, against B16 melanoma, with both exhibiting good in vitro T-cell responses but different tumor suppression abilities. METHODS: We compared activation markers and repertoires of T-lymphocytes from tumor-draining (dLN) and non-draining (ndLN) lymph nodes for the two peptide vaccines. RESULTS: We showed that the p30 vaccine had better tumor control as opposed to p20. p20 vaccine induced better in vitro T-cell responsiveness but failed to suppress tumor growth. Efficient antitumor vaccination is associated with a higher clonality of cytotoxic T-cells (CTLs) in dLNs compared with ndLNs and the convergence of most of the enriched clones. With the inefficient p20 vaccine, the most expanded and converged were clones of the bystander T-cells without an LN preference. CONCLUSIONS: Here, we show that the clonality and convergence of the T-cell response are the hallmarks of efficient antitumor vaccination. The high individual and methodological dependencies of these parameters can be avoided by comparing dLNs and ndLNs.

Clonal structure and the specificity of vaccine-induced T cell response to SARS-CoV-2 Spike protein.

Adenovirus vaccines, particularly the COVID-19 Ad5-nCoV adenovirus vaccine, have emerged as promising tools in the fight against infectious diseases. In this study, we investigated the structure of the T cell response to the Spike protein of the SARS-CoV-2 virus used in the COVID-19 Ad5-nCoV adenoviral vaccine in a phase 3 clinical trial (NCT04540419). In 69 participants, we collected peripheral blood samples at four time points after vaccination or placebo injection. Sequencing of T cell receptor repertoires from Spike-stimulated T cell cultures at day 14 from 17 vaccinated revealed a more diverse CD4+ T cell repertoire compared to CD8+. Nevertheless, CD8+ clonotypes accounted for more than half of the Spike-specific repertoire. Our longitudinal analysis showed a peak T cell response at day 14, followed by a decline until month 6. Remarkably, multiple T cell clonotypes persisted for at least 6 months after vaccination, as demonstrated by ex vivo stimulation. Examination of CDR3 regions revealed homologous sequences in both CD4+ and CD8+ clonotypes, with major CD8+ clonotypes sharing high similarity with annotated sequences specific for the NYNYLYRLF peptide, suggesting potential immunodominance. In conclusion, our study demonstrates the immunogenicity of the Ad5-nCoV adenoviral vaccine and highlights its ability to induce robust and durable T cell responses. These findings provide valuable insight into the efficacy of the vaccine against COVID-19 and provide critical information for ongoing efforts to control infectious diseases.

Unique features of the TCR repertoire of reactivated memory T cells in the experimental mouse tumor model.

T cell engineering with T cell receptors (TCR) specific to tumor antigens has become a breakthrough towards personalized cancer adoptive cell immunotherapy. However, the search for therapeutic TCRs is often challenging, and effective strategies are strongly required for the identification and enrichment of tumor-specific T cells that express TCRs with superior functional characteristics. Using an experimental mouse tumor model, we studied sequential changes in TCR repertoire features of T cells involved in the primary and secondary immune responses to allogeneic tumor antigens. In-depth bioinformatics analysis of TCR repertoires showed differences in reactivated memory T cells compared to primarily activated effectors. After cognate antigen re-encounter, memory cells were enriched with clonotypes that express α-chain TCR with high potential cross-reactivity and enhanced strength of interaction with both MHC and docked peptides. Our findings suggest that functionally true memory T cells could be a better source of therapeutic TCRs for adoptive cell therapy. No marked changes were observed in the physicochemical characteristics of TCRß in reactivated memory clonotypes, indicative of the dominant role of TCRα in the secondary allogeneic immune response. The results of this study could further contribute to the development of TCR-modified T cell products based on the phenomenon of TCR chain centricity.

A Rapid Method for Detection of Antigen-Specific B Cells.

The global SARS-CoV-2 pandemic has united the efforts of many scientists all over the world to develop wet-lab techniques and computational approaches aimed at the identification of antigen-specific T and B cells. The latter provide specific humoral immunity that is essential for the survival of COVID-19 patients, and vaccine development has essentially been based on these cells. Here, we implemented an approach that integrates the sorting of antigen-specific B cells and B-cell receptor mRNA sequencing (BCR-seq), followed by computational analysis. This rapid and cost-efficient method allowed us to identify antigen-specific B cells in the peripheral blood of patients with severe COVID-19 disease. Subsequently, specific BCRs were extracted, cloned, and produced as full antibodies. We confirmed their reactivity toward the spike RBD domain. Such an approach can be effective for the monitoring and identification of B cells participating in an individual immune response.

Targeted depletion of TRBV9+ T cells as immunotherapy in a patient with ankylosing spondylitis.

Autoimmunity is intrinsically driven by memory T and B cell clones inappropriately targeted at self-antigens. Selective depletion or suppression of self-reactive T cells remains a holy grail of autoimmune therapy, but disease-associated T cell receptors (TCRs) and cognate antigenic epitopes remained elusive. A TRBV9-containing CD8+ TCR motif was recently associated with the pathogenesis of ankylosing spondylitis, psoriatic arthritis and acute anterior uveitis, and cognate HLA-B*27-presented epitopes were identified. Following successful testing in nonhuman primate models, here we report human TRBV9+ T cell elimination in ankylosing spondylitis. The patient achieved remission within 3 months and ceased anti-TNF therapy after 5 years of continuous use. Complete remission has now persisted for 4 years, with three doses of anti-TRBV9 administered per year. We also observed a profound improvement in spinal mobility metrics and the Bath Ankylosing Spondylitis Metrology Index (BASMI). This represents a possibly curative therapy of an autoimmune disease via selective depletion of a TRBV-defined group of T cells. The anti-TRBV9 therapy could potentially be applicable to other HLA-B*27-associated spondyloarthropathies. Such targeted elimination of the underlying cause of the disease without systemic immunosuppression could offer a new generation of safe and efficient therapies for autoimmunity.

Four novel ELANE mutations in patients with congenital neutropenia.

Congenital neutropenia is a heterogeneous bone marrow failure syndrome characterized by a maturation arrest of myelopoesis at the promyelocyte/myelocyte stage. Cyclic neutropenia (CyN) and severe congenital neutropenia (SCN) are two main forms of congenital neutropenia. Genetic analysis has shown that heterozygous mutations in the ELANE gene encoding the neutrophil elastase are the major cause of these disorders. We investigated the prevalence of ELANE mutations in a group of 16 patients from 14 families with congenital neutropenia. Five patients had typical manifestations of CyN, and 11 patients had SCN. Seven different heterozygous ELANE mutations were found, including four novel mutations.

Naïve Regulatory T Cell Subset Is Altered in X-Linked Agammaglobulinemia.

The interplay between T- and B-cell compartments during naïve, effector and memory T cell maturation is critical for a balanced immune response. Primary B-cell immunodeficiency arising from X-linked agammaglobulinemia (XLA) offers a model to explore B cell impact on T cell subsets, starting from the thymic selection. Here we investigated characteristics of naïve and effector T cell subsets in XLA patients, revealing prominent alterations in the corresponding T-cell receptor (TCR) repertoires. We observed immunosenescence in terms of decreased diversity of naïve CD4+ and CD8+ TCR repertoires in XLA donors. The most substantial alterations were found within naïve CD4+ subsets, and we have investigated these in greater detail. In particular, increased clonality and convergence, along with shorter CDR3 regions, suggested narrower focused antigen-specific maturation of thymus-derived naïve Treg (CD4+CD45RA+CD27+CD25+) in the absence of B cells - normally presenting diverse self and commensal antigens. The naïve Treg proportion among naïve CD4 T cells was decreased in XLA patients, supporting the concept of impaired thymic naïve Treg selection. Furthermore, the naïve Treg subset showed prominent differences at the transcriptome level, including increased expression of genes specific for antigen-presenting and myeloid cells. Altogether, our findings suggest active B cell involvement in CD4 T cell subsets maturation, including B cell-dependent expansion of the naïve Treg TCR repertoire that enables better control of self-reactive T cells.

Molecular dissection of Penelope transposable element regulatory machinery.

Penelope-like elements (PLEs) represent a new class of retroelements identified in more than 80 species belonging to at least 10 animal phyla. Penelope isolated from Drosophila virilis is the only known transpositionally active representative of this class. Although the size and structure of the Penelope major transcript has been previously described in both D. virilis and D. melanogaster transgenic strains, the architecture of the Penelope regulatory region remains unknown. In order to determine the localization of presumptive Penelope promoter and enhancer-like elements, segments of the putative Penelope regulatory region were linked to a CAT reporter gene and introduced into D. melanogaster by P-element-mediated transformation. The results obtained using ELISA to measure CAT expression levels and RNA studies, including RT-PCR, suggest that the active Penelope transposon contains an internal promoter similar to the TATA-less promoters of LINEs. The results also suggest that some of the Penelope regulatory sequences control the preferential expression in the ovaries of the adult flies by enhancing expression in the ovary and reducing expression in the carcass. The possible significance of the intron within Penelope for the function and evolution of PLEs, and the effect of Penelope insertions on adjacent genes, are discussed.

RNA-Seq-Based TCR Profiling Reveals Persistently Increased Intratumoral Clonality in Responders to Anti-PD-1 Therapy.

Substantial effort is being invested in the search for peripheral or intratumoral T cell receptor (TCR) repertoire features that could predict the response to immunotherapy. Here we demonstrate the utility of MiXCR software for TCR and immunoglobulin repertoire extraction from RNA-Seq data obtained from sorted tumor-infiltrating T and B cells. We use this approach to extract TCR repertoires from RNA-Seq data obtained from sorted tumor-infiltrating CD4+ and CD8+ T cells in an HKP1 (KrasG12Dp53-/-) syngeneic mouse model of lung cancer after anti-PD-1 treatment. For both subsets, we demonstrate decreased TCR diversity in response to therapy. At a later time point, repertoire diversity is restored in progressing disease but remains decreased in responders to therapy in both CD4+ and CD8+ subsets. These observations complement previous studies and suggest that stably increased intratumoral CD4+ and CD8+ T cell clonality after anti-PD-1/PD-L1 therapy could serve as a predictor of long-term response.

Functionally specialized human CD4+ T-cell subsets express physicochemically distinct TCRs.

The organizational integrity of the adaptive immune system is determined by functionally discrete subsets of CD4+ T cells, but it has remained unclear to what extent lineage choice is influenced by clonotypically expressed T-cell receptors (TCRs). To address this issue, we used a high-throughput approach to profile the αß TCR repertoires of human naive and effector/memory CD4+ T-cell subsets, irrespective of antigen specificity. Highly conserved physicochemical and recombinatorial features were encoded on a subset-specific basis in the effector/memory compartment. Clonal tracking further identified forbidden and permitted transition pathways, mapping effector/memory subsets related by interconversion or ontogeny. Public sequences were largely confined to particular effector/memory subsets, including regulatory T cells (Tregs), which also displayed hardwired repertoire features in the naive compartment. Accordingly, these cumulative repertoire portraits establish a link between clonotype fate decisions in the complex world of CD4+ T cells and the intrinsic properties of somatically rearranged TCRs.

Measuring Intratumoral Heterogeneity of Immune Repertoires.

There is considerable clinical and fundamental value in measuring the clonal heterogeneity of T and B cell expansions in tumors and tumor-associated lymphoid structures-along with the associated heterogeneity of the tumor neoantigen landscape-but such analyses remain challenging to perform. Here, we propose a straightforward approach to analyze the heterogeneity of immune repertoires between different tissue sections in a quantitative and controlled way, based on a beta-binomial noise model trained on control replicates obtained at the level of single-cell suspensions. This approach allows to identify local clonal expansions with high accuracy. We reveal in situ proliferation of clonal T cells in a mouse model of melanoma, and analyze heterogeneity of immunoglobulin repertoires between sections of a metastatically-infiltrated lymph node in human melanoma and primary human colon tumor. On the latter example, we demonstrate the importance of training the noise model on datasets with depth and content that is comparable to the samples being studied. Altogether, we describe here the crucial basic instrumentarium needed to facilitate proper experimental setup planning in the rapidly evolving field of intratumoral immune repertoires, from the wet lab to bioinformatics analysis.

DSN depletion is a simple method to remove selected transcripts from cDNA populations.

A novel DSN-depletion method allows elimination of selected sequences from full-length-enriched cDNA libraries. Depleted cDNA can be applied for subsequent EST sequencing, expression cloning, and functional screening approaches. The method employs specific features of the kamchatka crab duplex-specific nuclease (DSN). This thermostable enzyme is specific for double-stranded (ds) DNA, and is thus used for selective degradation of ds DNA in complex nucleic acids. DSN depletion is performed prior to library cloning, and includes the following steps: target cDNA is mixed with excess driver DNA (representing fragments of the genes to be eliminated), denatured, and allowed to hybridize. During hybridization, driver molecules form hybrids with the target sequences, leading to their removal from the ss DNA fraction. Next, the ds DNA fraction is hydrolyzed by DSN, and the ss fraction is amplified using long-distance PCR. DSN depletion has been tested in model experiments.

Identification of a novel mutation in DKC1 in dyskeratosis congenita.

Dyskeratosis congenita (DC) is a rare congenital syndrome characterized by the triad of reticular skin pigmentation, nail dystrophy and mucosal leukoplakia, and the predisposition to bone marrow failure and malignancy. DC is genetically heterogeneous and X-linked and autosomal forms of the disease exist. Here, we report the clinical description and mutation analysis of a Russian family with X-linked DC. A novel mutation in DKC1 raised de novo in the maternal grandmother's gamete was found; this mutation is a 2 bp inversion in exon 3: NM_001363:c.166_167invCT (Leu56Ser).

A high-throughput assay for quantitative measurement of PCR errors.

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.

A digestive prolyl carboxypeptidase in Tenebrio molitor larvae.

Prolyl carboxypeptidase (PRCP) is a lysosomal proline specific serine peptidase that also plays a vital role in the regulation of physiological processes in mammals. In this report, we isolate and characterize the first PRCP in an insect. PRCP was purified from the anterior midgut of larvae of a stored product pest, Tenebrio molitor, using a three-step chromatography strategy, and it was determined that the purified enzyme was a dimer. The cDNA of PRCP was cloned and sequenced, and the predicted protein was identical to the proteomic sequences of the purified enzyme. The substrate specificity and kinetic parameters of the enzyme were determined. The T. molitor PRCP participates in the hydrolysis of the insect's major dietary proteins, gliadins, and is the first PRCP to be ascribed a digestive function. Our collective data suggest that the evolutionary enrichment of the digestive peptidase complex in insects with an area of acidic to neutral pH in the midgut is a result of the incorporation of lysosomal peptidases, including PRCP.

Normalization of full-length-enriched cDNA.

A well-recognized obstacle to efficient high-throughput analysis of cDNA libraries is the differential abundance of various transcripts in any particular cell type. Decreasing the prevalence of clones representing abundant transcripts before sequencing, using cDNA normalization, may significantly increase the efficacy of random sequencing and is essential for rare gene discovery. Duplex-specific nuclease (DSN) normalization allows the generation of normalized full-length-enriched cDNA libraries to permit a high gene discovery rate. The method is based on the unique properties of DSN from the Kamchatka crab and involves denaturation-reassociation of cDNA, degradation of the ds-fraction formed by abundant transcripts by DSN, and PCR amplification of the remaining ss-DNA fraction. The method has been evaluated in various plant and animal models.